o 
o 

(N 



o 



Oh 



Quantum Electrodynamics 
on Background External Fields 



> ■ Dissertation 

zur Erlangung des Doktorgrades 

(N 



des Fachbereichs Physik 
der Universitat Hamburg 



Piotr Marecki 



Hamburg 
2003 



Abstract 

The quantum electrodynamics in the presence of background external fields is de- 
veloped. Modern methods of local quantum physics allow to formulate the theory 
on arbitrarily strong possibly time- dependent external fields. Non-linear observables 
which depend only locally on the external field are constructed. The tools necessary 
for this formulation, the parametrices of the Dirac operator, are investigated. 



Zusammenfassung 

In dieser Arbeit wird die Quantenelektrodynamik in auBeren elektromagnetischen Fel- 
dern entwickelt. Die modernen Methoden der lokalen Quantenphysik ermoglichen es, 
die Theorie so zu formulieren, dass die auBeren Felder weder statisch noch schwach 
sein mussen. Es werden nicht-lineare Observable konstruiert, die nur lokal von den 
Hintergrundfeldern abhangen. Die dazu benotigten Werkzeuge, die Parametrizes des 
Diracoperators, werden untersucht. 
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CHAPTER I 



Introduction 

1.1. Formulation of the problem 

In this work Quantum Electrodynamics will be developed in which the Dirac field 
propagates on an external field background. Perhaps the best way to explain precisely 
what theory we have in mind is to look at its action. Suppose 1 : 

1 



S 



d x < i tp^f a d a il) — mtpij) + ei/j"f a ijj A a 



16/T 



-F ab F ab 



+ J a A c 



Here J a (t,x) denotes some external electromagnetic current which is a fixed function 
of time and space; ip denotes the Dirac field and A a the electromagnetic field. We 
divide A a into two parts, 

Aa = Af ua + A a , (1.1) 
where A class is a solution of the inhomogeneous Maxwell equations 

gbpclass = Alx j^ (L2) 

When substituted into the action S, the splitting (jl.lj) leads to an action of which the 
only dynamical variables are ip and A a : 

1 



S 



x 



</> {ii a d a + e~? a A c J ass -m)1> + e^ a ^A a 



167T 



-F ab F ab 



+ J a A a 



The variation with respect to ip an d A a leads to the Euler-Lagrange equations: 

(i^da + eAi ass - m) V = ~eA a ip, 

d h F* + <?FZr = ^ Ghat + J*)- 

Taking into account fjl.2jl . we get the following system 

(i^da + eAf ss - m) V = -eA a ip, 

That was the classical field theory. Quantum electrodynamics on external field back- 
grounds is the quantum field theory of the interacting Dirac and Maxwell fields. We 



4n the units h = 1 = c; the cgs-Gauss units are restored in appendix 1X1 
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first quantize the free fields, which obey the differential equations 



class 




(1.3) 
(1.4) 



a 



and then investigate their interaction following the steps of the causal perturbation 
theory 2 . We note that the division (jl.lj) is unique only up to the solutions of the 
homogeneous Maxwell equations, which thus can be included either as A class or as 
A. The classical, external current J a (t,x) is produced by some external sources (for 
instance by a heavy nucleus or by charged electrodes) and, by assumption, is not 
influenced by the (charged) quantized Dirac field ip. 



There are good reasons to investigate external field QED. The most important of 



difficult theory of quantum electrodynamics on a background curved spacetime (i.e. in 
the presence of gravitation). The problems posed by the latter theory are tremendous, 
yet nobody doubts it touches the central problem of theoretical physics which is to 
understand the relation between gravitation and quantum phenomena. Perhaps the 
most striking similarity between external field QED and QED on a curved space-time is 
the lack of a preferred vacuum state for the Dirac field. In the absence of a distinguished 
state many traditional concepts require (at least) a redefinition; to name some of them: 
the normal ordering of the field quantities or the concept of particles. Normal ordering 
is crucial if anything else than operators linear in the fields are to be considered 4 . 
The presence of particles in general causes certain characteristic responses of various 
detector arrangements. Particles are quasi-local excitations. However, if no vacuum 
is distinguished, it is impossible to say which configuration describes "excitations". 
Different basis states (the analogues of the vacuum) will give rise to different detector 
responses none of which can be distinguished as "preferred". There is no way to 
calibrate our detectors. 

The definition of non-linear quantities and the understanding of the association 
between detector responses and the presence of particles are not the only important 
issues which, when resolved in the external field QED, may help in the development 

2 The interaction Lagrangian of the perturbation theory is Ci — e : Tjj(x)^ a, ip(x) : A a - 

3 Apart from the fact that the external field QED provides the best currently accepted explanation of 

such a fundamental phenomenon as the spectrum of the hydrogen atom. 

For instance, one would like to investigate the currents, the definition of which requires however the 
normal ordering. 



1.2. General motivation 



which 3 , 



in our opinion, is the fact that this theory has much in common with the more 



of QED on a curved spacetime. After all, the external field theories are by no means 
fundamental theories. It is natural to expect the external field approximation to break 
down in certain regimes. The expectation is that the back reaction effects are to be 
regarded as a test if a given external field theory is a reliable approximation or not. 
The back-reaction in the context of external field QED means the additional (apart 
from J a which is the source of A class ) electromagnetic field produced dynamically by 
the quantum Dirac field if). To say that the external field approximation is justified 
means to regard the quantum fields propagating in it as test fields. Sometimes the 
back-reaction effects are naturally small as is for instance the reaction of an electron 
on the field produced by a macroscopic magnet. In other cases the back reaction 
is essential as for example in the free electron laser (FEL), where the synchrotron 
radiation emitted by a bunch of electrons interacts with this bunch and alters its 
dynamics 5 . In the external field approximation it is possible that every state produces 
some back- reaction effects, even "the vacuum" 6 . More importantly, in QED on a 
curved spacetime it would be interesting to know what is the energy-momentum content 
{T^ u (x)}ci of a certain "vacuum" state Q in the process of a collapse of a heavy star or, 
equally dramatically, does the black hole evaporate due to Hawking radiation. None 
of the above fundamental questions can reliably be addressed at the moment, partially 
because the evolution equations for the gravitational fields are highly complicated. We 
write partially, because there is another fundamental problem: what exactly is the 
back- reaction current /energy-momentum tensor, if no vacuum is distinguished 7 ? Thus 
- partially - the back-reaction question can be investigated more easily in external field 
QED, as the effect would add up to the given external field (Maxwell equations are 
linear) . 



1.3. Relation to other formulations of external-field QED 

The development of external field QED commenced almost simultaneously with the 
development of QED, in part due to the urge to describe atomic systems. The early 
investigations consisted almost exclusively of a double expansion: in 

A dass and in 



This and other main phenomena which occur at the FEL are reported eg. in the paper by S.V.Milton 
et al. |Mi01j . 

6 In quotation marks because there rarely exists a privileged state. 

This question is not trivial, even if a certain vacuum is distinguished - as in the no-external field 
case. Just that there is a unique quantity to subtract from the infinite expectation value does not 
mean that what remains is indeed the source of gravitation/electromagnetism. 



More precisely, the free fields were supposed to fulfill the equations 

{i^da - m) V = 0, (1.5) 
d b F£ = 0, (1.6) 

and the perturbation theory was developed with the external field as well as the quan- 
tum electromagnetic field on the same footing: 

d = e : #r) 7 >(£) :A a + e: ^(x) 7 >(x) : A*™ 3 . 

In such a way many processes of great physical importance have been explained, among 
others bremsstrahlung and e + e~ pair production in the field of a nucleus IBLP82| 
AB65]. Although physically one has learned a lot from those investigations, they 
implicitly assume that the external field is weak. Indeed a more profound theory has 
also been developed called the Furry picture or strong field QED MPS98, BLP82 . 
This theory is very similar to the one developed in this paper. The quantized free fields 
are supposed to fulfill the system of equations 

(tYd a + eAi ass - m) $ = 0, 
d b F a A b = 0, 

which is the same as ours, and the interaction is formally the same, 

though the Wick product in the Furry picture QED means the normal ordering which 
can be written as 

: ipi/j := ipip — (Q, tpifj Q), 

where Q is the vacuum (defined in a certain way). 

We aim at a better understanding of the quantum electrodynamics than the Furry 
picture QED gives. It is therefore necessary to put forward the weaknesses of the latter. 
In our opinion the main unsatisfactory features of this theory which are common to 
all of its formulations are: 

(i) In the definition of quantities nonlinear in the Dirac field (such as, for in- 
stance, the normal ordering required in the first order interaction processes) 
non-local objects are employed. This non-locality (elaborated upon in chap- 
ter El} manifests itself in a delicate way, namely, the observables defined as 



This approximation can also be recognized by the usage of free Dirac field propagators in the 
calculations. 



they are in the Furry picture QED do depend on the external field not only in 
the region of their support. For instance, a detector sensitive to the electric 
charge placed in a region C, 

£>(/) = J d 3 x : V*(x)V>(x) : /(x) with supp / = C, 

would be local if as an operator it depended at most on the external field in 
C. However, if :: means what it does in the Furry picture of QED, then 

even if the support of the variation 5A^ ass does not intersect with C. We 
emphasize the need for local observables. The states of the quantum field 
carry non local information, and that is a characteristic feature of relativistic 
quantum field theory. Locality means that at least observables should be free 
of acausal influences 9 , 
(ii) Almost all of the literature on external field QED assumes the external fields 
to be static. This unnecessary assumption carries with itself a false feeling 
of uniqueness of the vacuum representation which is employed. While it is 
true that the ground state on a static background is privileged as the state of 
lowest energy, we stress that not all external fields are eternally static. Some 
external fields are 10 turned on in the distant past of the experiment. It is 
highly likely that in such situations the state of the Dirac field at later times 
is not the ground state of the static potential. Also concepts like " adiabatic 
switching" of the external field require time dependence of the external field. 



We regard the drawbacks named above as very important, and we will not follow the 
Furry picture of QED any further. On the other hand, these drawbacks do not preclude 
the authors from deriving physically observable properties of matter, which are later 
compared with experimental results and yield a reasonable agreement. It is one of 
the remaining dilemmas whether the same or similar results can be derived from the 
improved foundations which we develop in this thesis. 

In a separate development the theory of quantum fields on curved spacetime has 
recently acquired a very satisfactory status. Indeed the works of many authors over the 



The precise formulation of this new type of locality has been given in |BFV01j . see also chapter ED 
°For instance, the trapping potentials in the ion traps. 



past decade resulted in an almost complete picture of the (interacting) electrodynam- 
ics on curved spacetime 11 |Wa94|. IBFOO|. IHW1|. IHW2|. IBFVOIj . A very modern 
approach allowed to remedy all the drawbacks similar to those named above. The 
renormalization theory in that scheme uses the language of distribution theory. One 
speaks of distributions, their extension to coinciding points and of the uniqueness of 
this procedure. This contrasts sharply with the language of divergent integrals and 
tricky extractions of the finite parts from them which are so common in the litera- 
ture on quantum electrodynamics. Although in the no-external-field context all these 
formulations of the renormalization lead to the same results the mathematical trans- 
parency of causal perturbation theory is encouraging [Sch96 . It seems that certain 
problems of uniqueness of the renormalization of the causal perturbation theory on 
external field backgrounds have not even been realized in the Furry picture QED. 
Our work thus attempts to achieve the following: 

(i) To formulate Quantum Electrodynamics on external field backgrounds in a 
modern way, using the methods of QFT on curved spacetimes together with 
the causal approach to the (perturbative) construction of interacting field 
theories. 

(ii) To construct the theory with a local dependence on the external background. 

(iii) To construct the theory on all possible external field backgrounds, even time- 
dependent ones. 

1.4. Structure of the paper 

The thesis contains seven chapters and five appendices. Here we shall briefly sum- 
marize their content. 

The second chapter is where our investigations begin. It deals with the quantiza- 
tion of the Dirac field in the presence of external field backgrounds. The first section of 
this chapter recalls standard properties of the classical Dirac field on external, possibly 
time-dependent potentials. Results on the selfadjointness and the type of the essen- 
tial spectrum are gathered there. In the second section we attempt to remove one of 
the main unsatisfactory features of the current formulations of the external field QED 
[Sha02, MPS98 . This feature is the restriction to one particular representation 12 of 
the free Dirac field algebra. We remove this unnecessary restriction with the standard 

11 To our great regret the various results have never been gathered together in a single reference. 
The physical (Dirac, Maxwell) fields are investigated by some authors, but the interacting theory (a 
version of causal perturbation theory) is only done for scalar fields. 
12 In the static case this is the ground state based representation. 



methods and results of the algebraic approach to quantum field theory Ha96 . In this 
apparently new application of these methods we rely upon quantum field theory on 
curved spacetimes, where such an application already proved to be useful. It is enlight- 
ening to realize that the global equivalence of states at all times, previously insisted on 
by many authors, is not necessary for the development of quantum electrodynamics. 
Although some observables, for instance the number operator or the total-energy ope- 
rator, are lost in this way, we are still able to describe the response of localized detectors 
which in our opinion link the theoretical description with experimental setups. 

We formulate the theory for a class of locally equivalent states - the Hadamard 
states. We allow all possible, non-singular external fields 13 . The concrete predictions 
can be obtained in any representation based upon an arbitrary quasi-free Hadamard 
state. Such states can be found on time-dependent environments. In particular it is 
relatively easy to construct Hadamard states, if the external field is static for some 
(possibly short) time interval. In the third section of the second chapter we recall the 
standard construction of the ground state representation. Mostly known results are 
gathered there. 

The third chapter deals with the quantization of the free electromagnetic field, 
which is the other basic field of quantum electrodynamics. In our theory, the free 
electromagnetic field A fulfills the standard Maxwell equations, and so the quantization 
procedure is standard (the Gupta-Bleuler method). 

In the fourth chapter, which is rather technical, we develop tools which enable us 
(in later chapters) to remove the other main unsatisfactory feature of the standard 
approaches to QED. This feature is the non-local dependence on the external field 
of these theories. The tools we develop are parametrices of the Dirac operator. To 
our knowledge they have not been extensively studied in the literature. Although 
the coefficients of those parametrices are written down in [DM75 , we have found it 
valuable to present our own derivation of them. It helps us later to study directly their 
short-distance limit, their scaling, uniqueness, dependence on the external field and 
their gauge covariance. Additionally, we expand the parametrix (which is a distribution 
of two variables) in a power series in the distance of its arguments. This straightforward 
computation allows us to see important things. For instance, we can foretell that 



On external gravitational backgrounds the Hadamard property as a spectrum condition rules out 
spacetimes with closed time-like curves - see (KRW97) . The case of non-smooth external fields 
requires a separate investigation. 



the instantaneous ground states (employed by some authors in the context of time- 
dependent external fields) are not Hadamard states which is a drawback of such states. 



The fifth chapter deals with the very important concept of the Hadamard prop- 
erty It describes the short- distance singularity structure of the allowed class of states. 
In this chapter we gather important theorems which assure that a broad class of 
states shares this property. We also recall the connection between two possible ways 
to define Hadamard states, namely, in terms of their short- distance singularity ex- 
pansion (the Hadamard series) of the two-point function and in terms of the wave 
front set of this two-point function. The equivalence of both definitions, first real- 
ized by M.Radzikowski |Ra96j for scalar fields and proven by S.Hollands Hol99j and 



K.Kratzert [KaOO for the Dirac field, is also reported here as it joins together various 
important parts of this thesis. 

The sixth chapter is in many ways the central one. It deals with the construction 
of non-linear field observables. These are the pointwise products of field operators 
smeared with test functions. There are at least two contexts for which non-linear 
observables are of fundamental importance. The first is the investigation of the current 
density and the energy-momentum density of the free quantum Dirac field. The other is 
the perturbative construction of interacting quantum electrodynamics. Our intention 
is to address both of these contexts. 

In the first section we recall the inductive construction of perturbative quantum 
electrodynamics. We use the framework of causal perturbation theory, which on the 
one hand is one of many formulations of the no-external-field quantum electrodynamics 
|Sch96j . and on the other hand is flexible enough to be applied to the construction of 
interacting quantum field theories on background spacetime manifolds [BFOO . The 
purpose of our investigations is to construct the building blocks of causal perturbation 
theory (the time-ordered products) in the lowest orders. In the second section we do a 
step in this direction by defining the algebra W of Wick polynomials of fermionic field 
operators. This algebra will also contain the time-ordered products which describe the 
interacting evolution in a finite order of the perturbation. 

The third section defines the most important concept of this thesis which is the 
local dependence of the observables on the external field. All of our important results 
are consequences of it. We motivate this requirement physically by showing it to be 
closely related to one of the foundations of general relativity. This foundation, the 
local position invariance, is well-tested experimentally and intuitively clear in content. 



Much of our subsequent work is a deduction from this very natural assumption . In 
the later sections we show, by means of simple examples, that both the normal ordering 
prescription and the renormalization subtraction scheme employed in known formula- 
tions of the external-field quantum electrodynamics are not local. Having established 
this, we proceed constructively and build the local Wick and time-ordered products in 
the lowest two orders of perturbation theory. 

In another development in the seventh section we discuss the definition of the 
current operator for the free quantized Dirac field. This is of prime importance for the 
investigations of the back-reaction process. The requirement of locality allows for the 
first time to reduce the huge ambiguity of its definition to a finite number of constants. 
Previously only differences of current densities of two states could be defined; here, we 
can define the absolute charge density of a given state. 

In the seventh chapter we begin to analyze the consequences of the local construc- 
tion of quantum electrodynamics. Here, we only show how various ingredients are 
combined together in calculations of the probability amplitudes of physically impor- 
tant processes in the presence of static external fields. A fair amount of work still 
has to be invested in order to derive concrete predictions of the theory. Specifically 
the construction of states in concrete situations is particularly cumbersome. The pur- 
pose of the seventh chapter is to outline the way in which concrete predictions can be 
obtained. 

The five appendices vary in importance and content. The first one deals with 
the electromechanical units which are employed in this paper. This issue can have 
important consequences, as for instance the dimensional analysis alone (together with 
the postulate of locality) reduces the ambiguity in the definition of the current density 
to three arbitrary numbers. 

The second appendix contains a brief exposition of the main theorems of micro 
local analysis which find their application in the chapter on the Hadamard form. 

The third appendix presents the vacuum representation of the Dirac field in the 
absence of external potentials which may help the reader not familiar with the external- 
field QED to recognize the familiar expressions in their generalization presented in 
chapter [TTJ 



The investigation of the dependence of quantum processes on the background fields on which they 
take place clarifies to a certain extent the meaning of local position invariance - cf. the remark in the 
section lVL3.ll 



The fourth appendix discusses our model of the atomic spontaneous emission of 
light. Although this model is only partially related to the main theme of this thesis, we 
believe it to give an important insight into the dynamics of the interacting theory. Here, 
on the basis of the Weisskopf-Wigner approach to the spontaneous-emission problem, 
we construct a model in which a system of two (non-relativistic) bound states of the 
electron interacts with the radiation field restricted to the vacuum and the one-photon 
sector. Instead of using the perturbation theory, we derive an equation (which is an 
integral equation) for the excited state's amplitude. In contrast to the perturbation 
theory, this equation is reliable also for large times. More importantly the calculation 
shows that the spontaneous emission is influenced directly by the two-point function 
of the radiation field. Due to the nature of this equation, it is straightforward to 
investigate various two-point functions of the radiation field, not only the vacuum one. 
We can, for instance, investigate the modification of the emission process on physical 
spacetimes (eg. Robertson- Walker) or in the presence of boundaries (eg. Casimir-like 
geometry) . 

The fifth appendix discusses the peculiarities of the construction of states of fermionic 
systems (i.e. representations of the CAR algebra). In a simple example we show what 
happens if the GNS construction is performed with a mixed "basis" state. The phe- 
nomena which occur are symptomatic of the problems which might occur in the general 
construction of representations of the CAR algebra in the presence of external back- 
grounds. 



CHAPTER II 



Quantization of the free Dirac field 

This chapter deals with the quantization of the Dirac field in the presence of external 
field backgrounds. It begins with a section on properties of the classical Dirac field and 
of the Dirac operator on various external field backgrounds. Most of the results are 
standard; we recall those which are particularly important for the further development 
of the theory. 

II. 1. Classical Dirac field 

In the following, let the Hilbert space be 

H = L 2 (R 3 ,d 3 x) 4 . 
The Dirac equation governs the time evolution of the vectors ip G H: 

id t m = H(t)ip(t), 

where 

H{t) = afl-idi + eAi(t, x)] - eA (t, x) + 

is the Hamiltonian, e denotes the electronic charge and m is the mass of the electron. 
The symbols a 1 , 7 stand for the Dirac gamma matrices 1 . The classical external 
electromagnetic field A^t, x) is assumed to be such that the Hamiltonian H(t) at each 
instant of time is self-adjoint on a suitable domain V(H) C H. 

II. 1.1. Theorems on properties of the Dirac operator. The Dirac operator 
in the presence of external fields can be split according to 

H = H + V(A), 

where 

H = -itfdi + m-f (II. 1) 

1 In the spinor (Wcyl) representation the gamma matrices are expressed in a simple manner by the 
2x2 Pauli matrices: 

-;■), / = ;°). 
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is the free part (independent of the external field A^it, x)); and 

V(A) = e[a%(t,x)-A)(t,x)] 

is the potential matrix, which is a multiplication operator. The matrix elements of 
V will be denoted by V^, i, j — 1, ... ,4. In the sequel we shall specify the domain 
of definition of H, V(H), and recall some results on its properties depending on the 
potential V(A). 

If V(A) is time-independent and Hermitian, then the following theorems hold true: 

Theorem II. 1 (Theorem 4.3 of |Tha91| ) 

If each matrix element of V is a smooth function of x, 

Vij G C°°(M 3 ), 
then H is essentially self-adjoint on Q°(M 3 ) 4 . 

Remark. The above theorem covers quite a substantial area of physical situations - 
the non-differentiable potentials are often only convenient approximations which are 
physically smoothed out at short distances. Even the Coulomb potential of an atomic 
nucleus is typically smoothed out inside the nucleus; a notable counterexample, where 
the singularity is not smoothed out, is the Coulomb field of an electron which appears 
not to be modified at any distance at all. 

The situation of external fields which possess a Coulomb-like singularity is covered 
by the following 

Theorem II. 2 (Theorem 4.2 of jThafllp 

If all the elements of the matrix are majorized by Coulomb-like terms, 
\Vij\ < ^- + b with x G K 3 \{0}, b > 0, a < 1 

then H is essentially self-adjoint on C^°(1R 3 \{0}); moreover, H is self-adjoint on the 
Sobolev space 2 H^R 3 ) 4 . 

Remark. The above theorem is sensitive to the constant which multiplies the po- 
tential. In the proof the theorem of Kato-Rellich [RS75J is utilized. In case of the 



2 The first Sobolev space H 1 is the space of L 2 -functions whose first derivatives are also square- 
integrable. 



Coulomb field 

A°(x) = eZ/\x\, 

after restoration of units, the theorem guarantees essential self-adjointness up to Z = 
68. For the Coulomb potential this is still not the maximal charge for which the 
essential self-adjointness property holds because of 

Theorem IL3 (Theorem 4.4 of jThafllp 

If the external field is electrostatic, i.e. A = 0, A = A°(x), and singular with the 
singularity not stronger than a/|x|, more precisely 

sup | |x| A°(x)| < a, 

xGR3\{0} 

then the corresponding Dirac Hamiltonian is essentially self-adjoint on Qj°(IR 3 \{0}), 
if a < ca/3/2 = 118.6 in atomic units (see appendix ID]). For larger values of a not 
greater than c = 137 there exists a unique selfadjoint extension of H whose domain is 
contained in the domain of V(A). 

If the potential is time-dependent and bounded, the evolution is described in terms 
of the Dyson series. In the interaction representation the evolution propagator is given 
by 



where T denotes the chronological order of the V operators and 

U(t, s) = e iHot U{t, s) e~ iH ° s , V(t) = e iHot V(t) e~ iHot . 
The unitary propagator U (t, s) fulfills the strong operator equations, 

U(t,s)iP = H(t)U(t,s)iP, 

-i-^ U{t,s)i/> = -U(t,s)H{t)i{>, 

for all ip in the domain of H, only if the family of interaction operators V(t) is strongly 
continuous in time 3 , for otherwise only weak, distributional solutions can be expected. 

In the quantization of the Dirac field on static external fields the question of the type 
of the spectrum of the Dirac operator is of interest. It therefore appears appropriate 



U(t,s) = Y,—r- / dn...dr n T 



This is the case, if the commutator [Hq, V(t)] is strongly continuous, cf. |RS75j chapter X.12. 



to recall general results which settle the question of the essential spectrum 4 a ess of the 
Dirac operator H . 

The essential spectrum of the free Dirac operator H is 



This property is stable under the addition of static potentials decaying at infinity. 
Theorem II. 4 

Let H = H + V(x.) be self-adjoint on a certain domain D(H), and let V(x) be decaying 
at infinity, 



then the Dirac operator H possesses the same essential spectrum as the free Dirac 
operator: 



Remark. As a consequence, it is appropriate to have in mind a picture of the spectrum 
of H consisting of two continua (free Dirac operator) and a ladder of bound states 
which can have {—mc 2 ,mc 2 } as the only accumulation points, as this will be the case 
in interesting applications. More on the essential spectrum of the Dirac field can be 
found in section 4.3.4 of |Tha 91j. 

II. 2. Construction of states on general external field backgrounds 

The purpose of this section is to remove one of the most unsatisfactory features 
of the recent constructions of the external field QED, which is the fact that they are 
founded upon a certain "vacuum" representation of the canonical anticommutation 
relations. This not only introduces an unnecessary assumption that the external field 
is static but also carries an unjustified claim that such a construction is unique and 
necessary. On the other hand, in the constructive approach to the quantum field theory 
the existence of many, even unitarily inequivalent representations of the canonical 
anticommutation relations is well-known. The temptation of applying the modern 
methods of local quantum physics ,Ha96; to the external field problem of quantum 
electrodynamics has resulted in the section that follows. 

This section discusses in detail how to find representations of the CAR algebra on a 
Cauchy surface. This is the algebra of fields ip(f), i>*(f), smeared on a Cauchy surface 

4 The essential spectrum a ess is the set of all accumulation points and infinite-degenerate eigenvalues. 





(II.2) 




with complex functions / G Q°(IR 3 ), together with their polynomials. The fields fulfill 
the Canonical Anticommutation Relations (CAR): 



If the external fields are static, the construction of a (vacuum) representation poses 
no particular difficulty and is described in many textbooks on QED. All that is needed 
in order to define such a vacuum state is a projection operator P + which describes the 
splitting of the underlying Hilbert space TL into the electron/positron subspaces TC±. 
Such a projection on a static background is provided by 5 

P+ = (l + sgnfO/2 

and is distinguished as it leads to the vacuum which is a ground state. 

Static external fields comprise, however, a narrow family of allowed environments. 
After all, hardly any field available in experiments is static for all times. With an 
important exception of the Coulomb field of an eternal charge, all external fields that 
are static during an experiment are rather generated earlier from the no-external-field 
environment, stabilized for the duration of the experiment and later turned off. It is 
important to realize that the ground state of the Dirac field in the static external field 
configuration is different from the state which was a ground state before the experiment 
and evolved in time while the fields were being turned on. A byproduct of this fact 
is the observation of G.Scharf |Sch96j that the adiabatically turned on Coulomb field 



should be modified on short distances by vacuum-back-reaction currents. Such an 
effect does not, however, occur, if the field is strictly static. 

What follows is the adaptation to the external field problem of the methods pre- 
sented in [Hol99j. We shall also make use of various results of [P S70j . 

II. 2.1. Introduction. Suppose the external field A^(x) is time-dependent. The 
classical Dirac equation 

[id t -H(t)) 1>(t,x) = 

with 

H(t) = -i 7 °7 < di + m 7 ° + e[a%(t, x) - A (t, x)] 



The representation produced by such a projection describes what is usually called "the Dirac sea". 
Here we assume that the Hamiltonian has an empty kernel. 



will be investigated. Suppose we restrict ourselves to the Cauchy surface t = t . The 
Hamiltonian at that instant of time, H(to), is an essentially self-adjoint operator on 
some dense domain in Ti. In what follows we shall make use of the spectral properties 
of H(t ); they are 6 

• the spectral measure d/x(k), 

• the (possibly generalized) eigenfunctions of positive ip£ (x) and negative ijj^ (x) 
frequency. 

The smeared two-point function of a state 

restricted to the Cauchy surface under consideration may be parameterized with the 
help of a positive, bounded operator < B < 1. Let us prescribe the action of this 
operator with the help of an integral kernel -B(x, y) 

(B gf(x) = J Sy B(x,yf A g A (y), 

for g eTi. Using the generalized eigenfunctions of H(t ), we define 

B nm (k,p) = J d 3 x d 3 y (^)i(x) 5(x,y) AC (i^%(y), (II.3) 

where the dagger denotes the conjugation in Ti (i.e. complex adjunction supplemented 
by a transposition) and the indices m, n were introduced and summed over in order to 
keep track of the positive /negative frequencies. The action of B can now be expressed 
in terms of B nm (k, p): 

(Bg) c (x) = J dfi(k)d^p) ($|) c (x)£U(k,p) «,s), (II.4) 

where (., .) on the RHS denotes the scalar product in Ti. As the field operator is 
smeared with the test functions which are elements of the underlying Hilbert space Ti: 

= f d 3 x /(x)^(x), 

so that 

uMf)P(9)) = J d 3 x d 3 y /t(x) c AC (x,y) g c (y), 



"In most cases the index k will be continuous if the eigenvalue corresponding to belongs to the 
essential spectrum of -ff(^o) and discrete if Vk i s a bound state of H(t ). The upper index which 
denotes the sign of the eigenvalue will in the following be a variable. The summation rule with 
respect to this index will assumed. 



we infer that the integral kernel of the operator B, B AC (x., y), coincides with the two- 
point distribution cj ac, (x, y). The two-point function of the free Dirac field is thus 
defined as 

uj^(f)r(g)) = (f,B g). 

It is the goal of the whole section to find a representation of the CAR algebra OL on 
certain Hilbert space T . In such a representation there exists typically a certain "base 
state" Q G T from which all other states (excitations) are generated with the help 
of (if they exist) creation operators. If we assume this "base state" to be quasi- free 7 
then a construction named after Gelfand, Naimark and Segal tells us how to find the 
representation n, the Hilbert space T and the "base vector" Q: 

II. 2. 2. GNS construction. The GNS construction starts from an observation 
that the elements of the CAR algebra OL (that is sums of products of field operators and 
their adjoints) may be regarded as vectors in a Hilbert space. Indeed, the quantized 
Dirac field operator if>(f), f G Ti, on static backgrounds is an element of OL. This 
operator (cf. section III.3J) corresponds to a vector in the Fock space T via 

V>(/) = a(p+f) + 6*(p_/) — v>(/)n = b*(P-m, 

which is a one-positron state with wave packet P-f. The idea of GNS is to construct 
the whole Fock space by applying products of field operators to the vacuum fl. At the 
beginning, however, we have only OL and the state a>; what is needed is: 

• the Hilbert space structure (the scalar product), 

• the creation operators which are elements of OL with appropriate smearing 
functions (in the example above b*(P-f) = tp(P-f), f G TC). 

The scalar product of A, C G O is provided by the state uj: 

(AQ, CQ) =u{A*C). 
It is semi-definite, because for some elements of OL it may happen that 

u(A*A) = 0. 



The quasi-free state has the property that all n-point functions can be expressed in terms of the 
two-point function. Only those n-point functions do not vanish which are the expectation values of 
an equal number of field operators and their adjoints, that is: 

w(V>(/i) ■ ■■i>(fn)ip*(9m) ■ ■■tp(9i))= Sn,m det[(/i, Bgj)]. 

What appears on the right hand side is the determinant of a matrix whose i, j-th elements are (/j, Bgj). 
The quasi-free property is thought in the literature to be a modest one, because many known states 
of quantum fields can be expressed by density operators in the representation based upon a quasi-free 
state. 



Those elements are precisely the annihilation operators which are the adjoints of what 
is sought in order to construct the whole representation space T . The search for the 
null space of u is by far the most non-trivial step of the GNS construction. This null 
space forms a left ideal 8 (named after Gelfand) of Ci in agreement with the fact that 
CAQ is always a null vector for all C G Ci, if A is an annihilation operator. 

If B which describes to is a projection operator on 7i, then it is relatively easy to 
find the annihilators in the one-operator subspace 9 of Ci. If A = ip{f) + V^CflO £ C^j 
then 

u{A*A)=0^ \P-f\ 2 + \P+g\ 2 = 0. 

Therefore ip(f+) = a(f+) and ip*(f_) = b{f), f± = P±f G TC, are the annihilation 
operators present in the algebra Ci. Products of their adjoints applied to Q generate 
the electron-positron Fock space. From the commutation relations of the field operators 
it follows 

{a,a} = 0, {b,b} = 0, 

{a,b} = 0, { a ,b*} = 0, 

{a(f), a*(g)} = (f, P+g), {b*(f), b(g)} = (/, P_g). 

The representation obtained in such a way (the standard vacuum representation) is 
irreducible. What happens, if B is not a projection, is studied in appendix |El as it 
leads away from the case of highest interest. 

II. 2.3. Time evolution, local and global quasi-equivalence. There are two 
ways to incorporate the time evolution into the theory. One of them is to view the time 
evolution as an automorphism of the CAR algebra Ci, which was the algebra of field 
operators on a Cauchy surface. The field operators were operator-valued distributions 
smeared with elements of the classical Hilbert space TC. The time evolution can be 
described by requiring the test functions to evolve in time. More precisely, the map 

iJ>{f)^1>(U*(t,to)f) 

preserves the algebraic relations and is reversible. Thus, it describes an automorphism 
of the CAR (which will be denoted by atjj). 

8 Algebraically this follows from the Schwarz inequality 

< \u({CA)*CA)\ 2 < uj{A*A)lu(A*C*CC*CA) = 

9 The annihilators in higher subspaces are a rather trivial addition, for they only enforce Pauli's 
principle (the appropriate statistics) in T . 



The other way to look at the time evolution is to investigate the two-point function 
of the state to. Quite frequently it can be decomposed, as in ()H.4(1 . with the help 
of some generalized eigenf unctions of some selfadjoint operator. The eigenfunctions 
V'k(x) will evolve in time just as classical solutions of the Dirac equation; the unitary 
propagator U (t, t ) describes this evolution. The two-point function at a later time 
t, cJt(x,y), can be found from the two-point function at the initial time t , u; to (x, y), 
with ^(x) replaced by U(t, t )V ; k( x )- We get 

tu t (f,g)=cu(U;f,U;g) = [uoat]{f,g), (11.5) 

where a t is the automorphism of the algebra described earlier. Both descriptions are 
equivalent; they correspond to the Heisenberg and Schrodinger pictures of quantum 
mechanics, respectively. 

In the literature there was a strong tendency to relate all the quantities to the 
no-external-field (call it Minkowski) situation, partially because one has regarded the 
unique Minkowski vacuum as an anchor to experimentally observable phenomena like 
the existence of particles 10 . The following question has been addressed: under what 
conditions on the time dependence of the external potential is the time evolution uni- 
tarily implementable in the GNS Fock space of the initial state? In other words, can the 
state u t for all times t be expressed as a vector state in the initial Hilbert space, or even 
in the Minkowski Fock space? The restriction of the allowed types of time-dependent 
backgrounds obtained in this way was quite dramatic |Sch96L ITha91L IRu77j . The 
time evolution can be implemented in the Fock space based upon the " Minkowski vac- 
uum" only for some 11 external fields of electric type, A = A (t,x), = 0. One way 
out of this trouble is to use the asymptotic notions. The state in the far future (the ini- 
tial state that evolved in the presence of external background) can be expressed in the 
initial Hilbert space under much weaker assumptions on the time dependent potential. 
In order to enlighten the way in which such questions are investigated we present the 
the Shale-Stinespring criterion (theorem 111.5(1 which settles the issue of quantum field 
theoretical implementability of classical unitary transformations. 



Implementation of unitary transformations. Consider a representation of 
the CAR defined in the Hilbert space of some pure "base state", which is uniquely 



In view of local quantum physics, particularly with the examples from quantum field theory on 
curved spacetimes (where there are no reasons to relate anything to the Minkowski vacuum), such a 
view must be regarded as obsolete. 

11 In |Ru77| . chapter 3D, which is our main reference for this result, the potential Aq is assumed to 
be a smooth function with rapid decay. 



described by its projection P + . Let U : H, — ► Ti denote a unitary transformation of 
the classical Dirac field 12 . It can be decomposed as follows: 

U = U+ + U+- + [/_+ + U- 

with 

C/ + = P+C/P+, \J_ = P_UP_, 

U + _ = P + UP-, U-+ = P_[/P + . 

On the algebraic level U can be promoted to an automorphism of the algebra, because 
the operators 

W) = VW) 

also fulfill the CAR. Suppose ip(f) denotes already the representative of the field ope- 
rator on a Hilbert space. The question arises, whether the field operator ip(f) can be 
represented as a unitary (Bogoliubov) transformation of ip(f), 

4>(f) = mf) = W(/)ir, 

for a given U. In this case, the original representation and the tilde representation are 
globally equivalent. The answer when this is the case is provided by the criterion of 
D. Shale and W.F.Stinespring [SS65J: 

Theorem II. 5 (Shale-Stinespring criterion) 

If the operators £/+_, £/_+ are Hilbert-Schmidt, then U is unitarily implementable. 
The state wo«[/ can be expressed in the Fock space of the state w as a vector state: 

n = un, 

where U is the implementation of U. 

Remark. The content of the Shale-Stinespring criterion is the following. We can 
always define the tilde representation of the CAR, 

where if)(f) is the base representation with the vacuum Q. We could also define the 
creation/annihilation operators a, b which would fulfill the same commutation relations 
as a and b: 

a(f)=i>(P + f), b(f)=r(P„f). 



U may denote the unitary propagator U(t,to) or the scattering operator S. 



If there existed a vector in the Fock space based upon f2, annihilated by a, b, 
then the existence of the unitary implementation would be automatic, namely the 
implementation of U could be defined on a dense subspace of T via 

U [V(/l) • • • HfnW(gi) ■ ■ ■ H9mM = mfl) ■ ■ ■ mfnW(U 9l ) . . . ^{U 9m )h. 

The Shale-Stinespring criterion tells us when the vector Q can be found. 

We note that the existence of a unitary operator U which intertwines between 
two representations of the CAR algebra is usually referred to as global equivalence of 
these representations. The Shale-Stinespring criterion tells us when two representations 
constructed upon pure "base states" whose projections are connected by a classical 
unitary operation are globally equivalent. 

Local equivalence. 

In local quantum physics global equivalence of states is a very strong criterion. 
There is absolutely no physical reason why all the allowed states should be globally 
equivalent to the, for instance, Minkowski vacuum. Indeed, phenomena like the thermal 
radiation present in the future of a gravitational collapse indicate that the condition 
of global equivalence is too strong. For instance, the Minkowski vacuum and the KMS 
(thermal equilibrium) state are not globally equivalent; if, due to some interaction 
with external sources (for instance cosmic background radiation), the Dirac field were 
to thermalize, then it would leave the folium of the Minkowski vacuum. 

In QFT one should, however, always restrict the folium 13 of allowed states. Instead 
of insisting on global equivalence one can investigate local equivalence. It says that in 
the bounded open regions the states should be normal to one another. In other words, 
there would always exist a density operator which expresses one of the states in the Fock 
(GNS) space of the other. Why to postulate local equivalence? Throughout this paper 
in order to define the non-linear observables it is necessary to restrict the investigations 
to the class of Hadamard states. Each two Hadamard states are locally normal 14 . In 
other words: only for locally normal states we are able do define (relative) non-linear 
observables like the current operator or Wick products or time-ordered products. 

The issues of global equivalence of states on the CAR algebra have been studied by 
R.Powers and E.Stormer in |PS70j and also by H.Araki [A ra71j . They have considered 



Folium means some reference state Q and all the states obtained as density operators in the repre- 
sentation connected with O. 

14 This has been proven in the context of curved spacetimes by R.Verch |Ver94] for scalar fields and 
for Dirac fields by C.D'Antoni and S.Hollands in [?]. 



the CAR on a Cauchy surface and have used their characterization with the help of a 
selfadjoint operator B, < B < 1 the way we have outlined it above. Among other 
things they have obtained the following result: 

Theorem II. 6 

Two states on the CAR algebra which are parameterized by the operators B\ and B 2 
(not necessarily projections) are globally equivalent if 

\\VB~i - Vb 2 \\h.s. < oo, || Vl - Bi ~ V 1 -B 2 \\h.s. < oo. 

Remark. If the operators B are projections then global equivalence reproduces the 
Shale-Stinespring criterion. To see this, let P + denote one of the projections and 
let B 2 = SP+S*, i.e. the second projection is connected to the first by means of a 
unitary transformation (Bogoliubov transformation). The Shale-Stinespring criterion 
says, that there exists a unitary transformation § on the whole Fock space related to 
P + which implements S, if and only if 

\\S + -\\ 2 BmSm =tr[(P + SP.y(P + SP.)] < oo. 

Using the properties of the trace one verifies that the above criterion is equivalent to 

tr(P+ - S*P + SP + ) < oo. 

On the other hand the criterion of global equivalence, in case the P's are projections 
as above, is 

tr[(P + - SP + S*){P + - SP+S*)} = tr(P+ - P + SP + S* - P+SP+S* + P+) < oo, 

and so it is the same as the Shale-Stinespring criterion. 

Let us now look at the local equivalence of two states tv 1 and uo 2 . We restrict those 
states to an open ball C on a certain Cauchy surface E t . The restriction of a state to 
C means that we multiply its projection P 12 by a characteristic function of the region 
C from both sides. We also restrict the underlying Hilbert space to H \c— L 2 (C, d 3 x) 4 . 
The states are locally equivalent, if the representations induced by the restrictions of 
the states to C are unitarily equivalent 15 . 

Remark. Even if the state u is pure, so that it is characterized by a projection B the 
restriction to C, u \c, will not be a pure state, i.e. will not be given by a projection. On 
the other hand it is possible do define a pure state on a given region C and extend it to 



'For all possible C on all Cauchy surfaces. 



the full algebra. Such a state, however, will not be locally equivalent to any Hadamard 
state 16 . 

The following theorem of |PS70j gives a condition for local equivalence of states 
Theorem II. 7 

Two states on the CAR algebra which are parameterized by the operators B\ and B 2 
are locally normal, if 

\\VB~i ~ ^/~B~ 2 \\h.s. < oo. 
The Hilbert-Schmidt norm is evaluated on the restricted Hilbert space Ti, \c] the above 
condition should hold for any open bounded region C. 

On the mathematical level it is clear that local equivalence is much weaker than 
the global one. As an example of locally equivalent states which fail to be globally 
equivalent let us consider the following 

Example. Let us take for simplicity two pure states described by projection operators. 
Moreover, let us assume both of them can be expressed as functions of some selfadjoint 
operator which possesses a continuous spectrum. The eigenf unctions of this operator 
will be denoted by t/> p (x). The second two-point function u 2 be built out of the first, u)\ 
by subtracting a projection operator A which is smaller than More precisely, the 
projection will be characterized by a characteristic function 17 x(p) °f compact support. 
We set 18 



w 2 (x,y) = wi(x,y) -J dp x(p)^ P (x)^*(y) = w x (x,y) - A(x,y). 
Global equivalence would imply 

II -Bi - B 2 \\ H .s. < oo, 

where B's denote the projections which characterize both states. Clearly the integral 
kernel of the difference B\ — B 2 is just A(x, y). Its Hilbert-Schmidt norm is 

\\B 1 -B 2 \\ H . 3 . = J d 3 xd 3 y |A(x,y)| 2 = J d 3 p d 3 k x(k)**(p) {ip k ,ip p ) {ip p ,ip k ) , 



Statements of this sort, although without a proof here, are a common lesson of local quantum 
physics Ha96, Fre94 , and are related to the fact, that local algebras Ct(O) are factors (possess a 
trivial center) of the type II I\. This ensures that there are no pure states on OL(0) locally normal 
to any quasi-free Hadamard state. 

17 The possible values of x(p) are an d b because we want it to be a projection, moreover, we assume 

X is such that the operator wi is positive (i.e. A < u>i). 

18 This must be done in such a way, that uji as an operator remains positive. 



where (ip p , ipk) is the scalar product, from 7i, of two scattering states, which is equal 
to S(p — k). Thus the above Hilbert-Schmidt norm is infinite. 

On the other hand, the question of local equivalence, may be answered with the 
help of the inequality proved by Powers and Stormer: 

II \[~B\ — v^H/f-.s. < 1 1 -Si — B 2 \\tr- 
Let g(x) denote the characteristic function of the region C. We estimate the RHS: 

||#i - B 2 \\tr,c = \\goA og\\ tr = || A O gr||^ s , 
because A is a projection. The Hilbert-Schmidt norm of A o g can be evaluated: 
II ■ • • Wh.s. = J d 3 P d 3 k x(p)x(k) |(^ P ,^ fc )| 2 < 

< J d 3 p d 3 k x(p)x(k) \\gip p \\ 2 ||# fe || 2 = J d 3 p x (p) ||#p|| 

by the Schwarz inequality. As the norm of the scattering states ip p restricted to the 
region C is a bounded function of p (at least if the external fields are not singular) and 
x(p) is a function of compact support, we infer that the states u)\ and uj 2 are locally 
equivalent in every bounded open region C. 

The preceding discussion has made it clear that in the external-field electrodynam- 
ics there is an enormous freedom of choosing "a basis state" for the construction of the 
representation. In what follows we shall investigate an example, where we compare 
two natural candidates for such "basis states", namely the state which evolved out of 
a preferred state in the past (when for instance the external fields were absent) and 
the state termed "instantaneous ground state" which is prescribed by a projection on 
the positive-frequency subspace of the instantaneous Hamiltonian. 

II. 2.4. An example (equivalence of states, instantaneous vacua). 

Example. We are going to consider a model of the Dirac field on such external field 
background that there would be positive- and negative-frequency bound states 19 . In 
order to simplify the computational side we assume the energies of the n-th electronic 
bound state to be equal in absolute to the energies of positrons in the bound state 
of the same quantum number. Additionally, there would be a time dependence of 
the Hamiltonian. Again, for simplicity, we assume that the interaction couples only 



We think of this model as of fermions bound by gravity. Indeed the model will resemble the 
construction of representations of the CAR given for Robertson- Walker spacetimes by S.Hollands 
|Hol99| . 



the bound states of the same n (but positive and negative subspaces are allowed to 
interact). The Hilbert space of the system, H, is 

Then a convenient orthonormal basis of this space is formed by the vectors 
where n G N. The time evolution of the model is governed by the Hamiltonian 



where u's are Pauli matrices and /(£) denotes an arbitrary function of time. 

Let the first state, o>i, be described by a projection on the positive frequencies at 
t — 0. Suppose that the function / vanished at that surface. Thus at that surface we 
have 

n 

where 

* - ( 1 °^ 

\0 0, 

and 20 ^ is the orthonormal basis introduced earlier which is also the basis of the 
eigenvectors of the Hamiltonian with / = 0. 

The goal is to compare the state uji(t), evolved from to t according to the time 
evolution governed by the Hamiltonian, with some instantaneous ground state defined 
via projection on what at t appears to be the positive-frequency subspace. Thus, the 
second state is defined as 

n 

where 

B, = ( 1 °^ 

and ip[ n are instantaneous eigenvectors of the Hamiltonian at time t: 
W sin(o;/2) / V — % cos(o;/2) / 



u Here we introduce a new index r for brevity. Its possible values are + and — which correspond to 
the upper and lower vector components respectively. 



where a = a(n,t) and 

sin(o;/2) 
cos(a/2) 



1 / t r 



\ 2 \ vw+m. 



n 



\ 2 \ vwrm. 



What is needed in order to check the local equivalence 21 of both states is to express 
the (evolved) ip r n (t) in the new basis ipl . As a matter of fact, it is generally impossible 
to find the unitary time evolution operator explicitly - we have to solve the Schrodinger 
equation: 

[idt-H(t)]if> n (t) = 0, 

which in the interaction picture 
leads to the equation 

n 

that is 

f(t) 

iip n = — ■ [cos(2E n t)a 2 + sin(2£' n i)cri] (p n . 

As the operator in square brackets depends on time 22 , it is not possible to write down 
an explicit solution of this equation generally, i.e. without investigating the concrete 
structure of f(t) and E n . 

Due to the above observation, we are forced to consider special cases in hope of 
gaining some insight into the construction of various states on time-dependent external 
fields. 

Suppose the state is defined by a projection operator B(0) and that the interaction 
term vanishes, / = 0. Then at a later time the state is defined via a time-evolved 
projection which may be written as 



B (t) = z~ lEntas r n B rs r: e 



Locality means that we assume finiteness of the norm of each ip^- 
22 Thus, the generator of the evolution cannot be written as a certain function of time times some 
time-independent operator. 



B+ B,e 2iEnt 



AB 



In other words, B(t) restricted to the subspace n is 

B ^ = [^B^+e- 2 ^ B 

Is B(t) equivalent to 5(0)? As both operators, B(t) and B(0), are projections the 
condition on the equivalence of the respective representations, 

\\y/W) ~ VWWh.S. <OC, 

may be replaced by 

\\B(0) -B(t) \\ 2 H . S . <oo. 
The difference inside the norm is easily seen to be equal to 

B + _(e 2iEnt - 1)' 

B_ + (e~ 2iE ^ -I) 

with B + - = B*_ + (they are numbers). We find 

(AB)*(AB) = 4sin 2 (£ n t) 

on the subspace n. The Hilbert-Schmidt norm may now be evaluated, say, w.r.t. the 
basis i\) r n of the Hilbert space H. We find 

II • Wh.s. = tr ( 4 sin2 i E nt)\B + -\ 2 \n) (n\ 

\n=0 

The trace is then calculated, and the condition for the equivalence of B(t) and B(0) 
becomes 

y^sin 2 (E ra t) sin 2 [2a(n)] < oo, 

n 

where we have employed the explicit parametrization of B at the subspace n, namely: 

cos 2 [a(n)] \ sin[2a(n)] 



B(0) = Bit) 




|sin[2a(n)] sin 2 [a(n)] 



where a is an angle which may vary with n. 

We observe the following: not every projection operator is compatible even with 
the free time evolution. Suppose the spectrum of the Hamiltonian is arbitrary, then 
the angles a(n) have to be chosen in such a way that the operator is Hilbert 
Schmidt which means that they have to decrease for large n quicker than 1/n. 

On the other hand, if the energies E n accumulate around for large n quickly 
(similarly, E n has to decrease quicker than 1/n) than the equivalence of B(t) and B(0) 
will be guaranteed irrespective of B(0). 



We see, therefore, that both concrete ingredients, the projection operator B(0) and 
the spectrum of the free Hamiltonian, may play a role in the equivalence of states at 
different times. 

II. 2. 5. Physical meaning of local and global equivalence. Physically the 
notions discussed in this section mean the following: global equivalence means that 
the vacuum corresponding to the state 0U2, when expressed in the Fock space of ui, 
contains finitely many "particles"; local equivalence means that the relative energy 
density (defined by a point-splitting procedure) is finite. 



II. 3. Quantization on the static external field backgrounds 

This section contains an important special case of the general construction of states, 
namely, the construction of a ground state of the Dirac field on static external back- 
grounds. An even greater specialization, namely the case of absence of any external 
backgrounds has been included in appendix for comparison. We recall the con- 
struction of this representation, because the static external fields provide the simplest 
non-trivial example of the external-field problem. 



II. 3.1. Negative and positive frequency subspaces of H. If the external 
fields are time-independent, it is possible to define a splitting of Tt into positive- and 
negative-frequency parts which will be invariant under the time evolution. More pre- 
cisely, let us define the operator 

r = sgn# = — , (II.6) 

if the Hamiltonian H has a vanishing kernel. With the help of r we define the projec- 
tions 

P± = \(1±t), (U.7) 
and the splitting of the Hilbert space Ti into the positive/negative frequency subspaces 

P±H± = H± 

H = H + ®H- 

If the zero belongs to the spectrum of H, it is necessary to settle the splitting of ker H 
by a separate prescription. The splitting presented above is invariant under the time 
evolution, because [P + ,H] = and the Hilbert spaces P + H defined at different times 
are the same. 

Remark. The operator r does not preserve the localization; namely if one takes if) ETt 
to have a compact support within a certain region N C M 3 , then in general if) + = 
P+ip G Ti+ does not have a compact support. Consequently, the term "relativistic 
quantum mechanics" has little sense, as the electron which should be described as a 
positive-frequency solution of the Dirac equation propagates acausally. 

On the Hilbert space TC we may also investigate the charge conjugation operator 
which is defined (in the spinor representation) by 

Cif> = i^ 2 if)\ 



for ip G Ti (ij)' denotes the Hermitian conjugation of ip, and 7 2 the second Dirac 
matrix). The operator C possesses the following properties: 

• it is anti-unitary: (Cip,Cx) = (x, V0i 

• double application of C is the identity 23 : C{C^) = ip 

• C "conjugates the charge", that is 

CH^C' 1 = -H(-e). 



II. 3. 2. Representation of the CAR algebra. 

In the following we shall recall the construction of the Fock representation which is 
based upon the splitting prescribed by P + . The Hilbert space on which the operators 
will be represented, the Fock space, is constructed with the help of creation operators. 
Let fi denote the orthonormal basis vectors of the space 7i + and <?j the orthonormal 
basis vectors of TC-. One defines a Hilbert space jF( n,m ) as a space spanned by the 
vectors of the form 

a*(f 1 )...a*(f n )b*(g 1 )...b*(g m )n. 

The space J r< y n < rn ) is equipped with the scalar product, whose form follows from the 
classical scalar product (., .) and the anti-commutation relations: 

{«(/,), a*(/,))} = (fijj), {&(<*),&*(<&))} = (9j,9il 

{a(/,), &*(<&))} = 0, Hfi),b(gj))} = 0, 

{a(/ i ),a(/ i ))} = 0, {%,),%,))} = 0. 



The Fock space is a direct sum of all the spaces 

oo 

T = jr(n,m) 
n,m=0 

where the case m, n = denotes the vacuum. 

The representation of the field operators can now be given: 

Definition II. 8 

The field operator is an operator valued distribution which takes as a test function 
all the 4-component Qj°(lR 3 ) functions on the surface of constant time (such functions 
are dense in H). If h e Qf(lR 3 ) then we express the field operator in terms of the 



'In the spinor representation this is due to {ij 2 )* = i"f 2 . 



creation/annihilation operators: 

ip(h) = a(P + h) + b*(P-h). (II.8a) 
The (Fock space) adjoint of the field operator is given by 

if>*(h) = a*{P + h) + b(P-h). (II.8b) 

The field operators are bounded, antilinear and fulfill the canonical anti-commutation 
relations (CAR), 

Mf),r(g)} = (f,g)i (n.9) 

where (., .) denotes the usual scalar product in TC, i.e. 



(f,g) = / d 3 x f(x) 1 °g{x) = / d 3 x /(x) 1 " g(x). 



The definition of the field operators completes the construction of the vacuum 
representation of the CAR algebra (|II.9J) 24 . 



II. 3. 3. Implementability of the unitary evolution. 

The next question we are going to discuss is whether the free time evolution which 
is given on the classical level by a unitary group U t can be implemented as a unitary 
group on the Fock space. This question is important because we would like to relate 
the field operator at the time t, ip t to the filed operator at the time t = 0, if) by an 
action of a unitary operator: 

Mf) = vMf)w t , 

which is just the Bogoliubov transformation. We may define 

where Ut denotes the classical unitary evolution operator, and the star denotes the 
H-adjoint operation. In order to answer the question under consideration we just have 
to apply the Shale-Stinespring criterion (theorem 111.5(1 . In the case of time-indepen- 
dent potentials, if the splitting of the Hilbert space TC is given by some function of the 
Hamiltonian, for instance 

P± = ^(l±sgntf), 



The abstract symbols ip(f) have acquired a concrete realization as operators on a Hilbert (Fock) 
space. 



it is particularly easy to implement the time evolution. The classical evolution is given 
by the one-parameter unitary group 

U(t) = e- itH , 

and, consequently (because P + commutes with H) 

u + _ = = u_+. 

The Shale-Stinespring tells us that the classical evolution in the case under considera- 
tion is always unitarily implementable. 

II. 3.4. Ground state. Let us consider the situation of a vacuum state in static 
potentials vanishing at infinity. For such potentials there are finitely many points of 
the spectrum of the Dirac-Hamiltonian H in the vicinity of E = 0. 

If P + denotes the projection operator on the subspace of positive eigenvalues of 
H, then, according to the preceding construction, we may obtain a vacuum state. 
However, there is a freedom in the choice of the splitting point (here E — 0). In what 
follows we will argue that different splittings lead to a globally equivalent representation 
(Bogoliubov transformation) and that only E = leads to the state of minimal energy 
which will be called the ground state. 

Lemma II. 9 

Let Ti = liffij be the Hilbert space spanned by the orthonormal basis /j, i G Z. Let 
the splitting Ti = Ti + © Ti- be prescribed by the projection operator on the vectors 
with positive indices: 

P+fi = fi Vz>0, 

both Ti + and 7i_ are assumed to be infinite-dimensional. A splitting of this sort may 
correspond, for instance, to a splitting into positive- and negative-frequency subspaces 
of a Hamiltonian H . Suppose, furthermore, that the system evolved in time and the 
classical evolution is given by the unitary operator 

Uf i = f i+1 V^GZ. 

Then the unitary transformation U can be implemented in the original Fock space. 
In other words the state u corresponding to P + and the state 25 Cj = lo o a t are globally 
equivalent. 



See equation (|II.5fl • 



Proof. The vacuum fl of the representation associated with uj is defined by the 
relations 

a{p + m = 0, b(P-f)Sl = 

for all /6H. The new representation of the CAR is defined via 

Hf) = mf) 

for all / e Ti. Specifically we obtain the following Bogoliubov mapping of the cre- 
ation/annihilaiton operators: 

a(P + f) = a(U + f) + V(U- + f), 
b(P_f) = a*(U + _f) + b(U_f), 

which are the Bogoliubov relations. From the definition of U we infer 

U+ _ = P + UP- = P , 
U-+ = P-UP+ = 0, 

where P is the projection onto f . We may summarize this as follows: 

a(P + f) = a(P + f), 
b(P.f) = a*(P f) + b{P.f). 
Thus the only nontrivial Bogoliubov transformation is 

6(/o) = <*•(/„). 

In this way we obtain the mapping of the two representations. What remains to be 
shown is that the state u can be expressed as a vector in the Fock space constructed 
upon the vacuum Q, which is annihilated by a and b: 

~a(P + f)Q = 0, ~b(P-m = 0. 

That this is the case may easily be verified, namely 

Q = a*(f )n, 

the operator b(f ) annihilates fi, because it is equal to a*(/ ) and we have Pauli's 
exclusion principle. □ 

If the states with different splitting of TL are expressed in the same Fock space (as 
in the preceding lemma) we may ask: which of them has the minimal energy? The 
energy operator (the self-adjoint generator of the time evolution unitary group) up to 



an arbitrary constant equals 

H = J2 l E +( n ) «*(/»)"(/») " 6*(/-n)6(/-n)] + ^+(0) a*(/o)a(/o), 

n 

where E±{n) denote the eigenvalues of the time-independent Hamiltonian H corre- 
sponding to the eigenvectors f n , n e N. We assume that E + (n) are positive numbers 
in contrast to E-(—n). If we compare the expectation value of EI in different states 
(regarded, due to the global equivalence, as different vectors in the same Fock space), 
we arrive at 

Lemma 11.10 

Of all possible states characterized by projection operators on the subspaces of fre- 
quencies greater then E the state with E = has the lowest energy. 

Proof. In the case, where the splitting is done with respect to E — 0, the ex- 
pectation value of EI evidently vanishes. Suppose that we investigate Q where one 
positive-frequency vector f is included in P_. We calculate 

(P,mn) = E(f ) (o*(/ ) n,a*(/ )fl(/„) a*(/ ) fi) = E(f ) 
which is greater than zero. □ 

Remark. The above lemma tells us that it is not at all arbitrary where to put the cut 
into positive/negative frequencies. Indeed, any choice other than E = leads to states 
of the Dirac field which are not the states of lowest energy in the Fock space. All the 
states of lower energy, however, differ in charge from Cl so that it would be difficult to 
imagine an electrodynamic process which could extract this energy. 

II. 3. 5. Time-dependent external fields. Suppose that we have external fields 
which vanish in the far future and past. In such a case, instead of trying to imple- 
ment the (classical) unitary propagator U (t, s) in the free (no external field) Fock 
space, which is rarely possible, one might try to implement the scattering matrix S 
only. The Shale-Stinespring criterion (theorem 111.5(1 . when applied to the scattering 
matrix, provides a (rather weak) tool to decide for which external fields the unitary 
implementation § of S exists: 

Theorem 11.11 (Theorem 8.25 of |Tha91j and 5.1 of |Sch96j l 

Let W n (t,x) denote the strong derivatives of the potential 

V(t^) = e 1 °[ 1 a A a (t^)} 



with n = 0,1,2, that is W n (t,x) = dV n /dt n . Let the family W n (t,x) be strongly 
continuous in t. Furthermore, let the Fourier transform of W n (t,x.) satisfy 

/oo 
\\W (t,. )\\ k dt< oc 
-oo 

for each n — 0, 1, 2 and all k = 1, 2. Then S + - and are Hilbert-Schmidt operators, 
and thus the scattering operator S is unitarily implementable. 

Remark. In the course of the proof of the above criterion it is easily seen that, given 
the assumptions of the theorem, the first-order operators 

= -ieP + J ds e^ 0S 7 °7 • A(s, x) e~ iHoS P-, 

S w + = -ieP_ J ds e^ 0S 7 °7 • A(s, x) e~ iH ° s P+, 
are also of the Hilbert-Schmidt type. 



CHAPTER III 



Quantization of the electromagnetic field 

The theory of the free, quantum radiation field (i.e. electromagnetic field without 
sources) will be presented in this chapter. We shall follow the standard method of 
Gupta and Bleuler |Gup50, Ble50; explained in a brief manner in |IZ78j . 



III.l. Quantization of the vector potential 

The vector potential A^(x) is the quantity that will be quantized. We introduce a 
sesquilinear form (., .) on the space of functions of fast decrease (Schwartz space) 

= "(2^)5 /^»(pV' 

This form is not positive-definite. The space of four-component Schwartz functions 
together with the sesquilinear form (., .) will be called the Krein space JC. 

The standard decomposition of the four- vector electromagnetic potential (classical), 
which is a basis of the subsequent quantization, reads 

1 C rl 3 k 

A,(x) = —3 / -= e«(k) R(k) e «* + a a (k) e~**} , (III.l) 

where ko = |k|, and e a (k) denote four real polarization vectors which are supposed 
to be orthogonal to one another and normalized so that three of them are spatial, 
namely e l ,e 2 ,e 3 , and e° is timelike. Moreover, as all spatial polarization vectors will 
be k-dependent, we choose e , e 2 to be of the form e ,2 = (0, e ,2 ) and to be orthogonal 
to (0,k): 

e Q (k) e /3 (k) = 5 at3 , 
e a {k) k = 0, 

for a,j3 = 1,2. Altogether (e a ) M form a four-tetrad: 

Vap (e°)^(e Q ) 1 , = r}^. 
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We can now quantize the theory turning a a (k) into operators 1 which are supposed 
to fulfill the commutation relations 

[a a (p),a;(k)] = -77^<J(p-k) (III.2) 

and to act on the Fock space constructed out of the Krein 2 space /C. Such a Fock- 
Krein space Tk, is constructed in a standard manner. Its n-particle subspace is the 
symmetrized tensor product 



= S n K ® . . . <8> K. 



r(n) 
K 

The operators a a , a* act on as standard creation operators: 



a* if) Snh ® . . . ® h n = Vn+ 1 S'n+x/ ® h x ® . . . ® /i n , (III.3) 

1 n 

a(/) Sn/ii ® . . . ® h n = —= )(f, hi) Sn-ih ® . . . ® ® . . . ® /i„, (III.4) 

V n ^ 

where the tilde denotes the omission of h{ in the symmetrized product. The vacuum 
is defined via 

a a {k)Q = 0. 

This is the standard Gupta-Bleuler representation of the field operators 3 . Note 
that the adjoint * means the Krein-adjoint. With such notations the operator Ao(f), 



Mf) = J d 4 x f(x)A (x), 

is Krein symmetric: 

(Mf)F,G) = (F,A {f)G) VF, G e T K . 



Those symbols will be operators when smeared with appropriate test functions. 

2 Here we can understand why the indefinite product space K. has been introduced. In quantum theory 
we always search for a representation of the observable quantities as selfadjoint operators on a Hilbert 
space. Here, however, such a representation cannot be constructed covariantly (i.e. to respect the 
covariant CCR l|III.5|l l. The reason for this is that the relation [oo,Oq] = — 1 cannot be fulfilled by 
selfadjoint operators on a positive product (Hilbert) space. 
3 The Fock-Krein space is sometimes called Gupta-Bleuler space. 



III. 2. The Lorentz condition and the physical Hilbert space 

The Krein product on /C is not positive-definite 4 . However, up to now the vector 
potential does not fulfill the Maxwell equations. It does fulfill 



which however follows from 



UA^x) = 



8 V F» V = 



with 

F — 8 A — 8 A 
only if the Lorentz condition is imposed 

dpA" = 0. 

It is very fortunate that we can remedy two shortcomings in one stroke, namely, a 
way to impose the Lorentz condition as a strong operator equation at the same time 
deletes the negative "norm" states from the Krein space. More specifically, the physical 
Hilbert space H p is defined with the help of 

1 f rl 3 k 

X (x) = 8»AM +) = -=s / ~m <2(k)*" *a(k)e-*° 

V27T J V^fc 



via 5 



H p = {Fe!C: x (x)F = 0}. 

The condition 6 

d^A^x)^ n = 
with our choice of the polarization vectors means 

(oq - a 3 ) F = VF G H p . 



Which means there are vectors F for which 

(F,F) <0. 

5 This condition is (as x( x )) Lorentz-invariant. 

6 As we have already said this is a strong operator equation. It is weaker than an operator equation 
(which would assert d' 1 A ll (xY + ^ = on an algebraic level). It is, however, stronger than the weak 
operator equation (F, d^A^x)^ G) for all F,G £ T K . 



A way to see the structure of T~C P is to introduce a new set of creation/annihilation 
operators: 

b = -^( a 3 ~ do), 
c = ~^( fl 3 + fl o), 

which fulfill 7 

[Ml = 0, [c,c*]=0, 
[b,c*] = l. 

The Lorentz gauge condition is a strong operator equation bF = for all F G 7Y P . 
Clearly the states from Ti, p may still contain longitudinal and scalar excitations in 
coherent pairs as 

b {b*) n n = 0. 

On the other hand no excitations created by c* can be present, due to 

b (c*) m n ^ 0. 

Therefore the physical Hilbert space H p is a Fock space generated by a*, a^, 6*. There 
are still zero norm vectors in Tip. In fact all vectors which contain scalar-longitudinal 
excitations have zero norm, which follows from 

(b* n n,b* n n) = (n,b n b* n Q) = o. 

Those states can be factored out, which means that we may assume them to be equal 
to the zero vector. In that way we obtain the transversal Hilbert space H tr in which 
only the physical, transversal photon polarizations are present. The Krein product, 
when restricted to H tr , is positive-definite. Henceforth, we may work in a standard 
Fock-Hilbert space of transversal photons. However, one should not forget the H tr 
is a factor space: a Lorentz transformation may lead to the appearance of zero-norm 
components which, due to factorization, are equivalent to the zero vector. 

The other important reason not to forget the non-physical polarizations is the fol- 
lowing: the commutation relations (|III.2j) lead to the standard commutation relations 
for the potentials: 

[A^x),A v {y)] = >A(i, y), (111.5) 



Here we omit the arguments of the creation/ annihilation operators as we only intend to illustrate 
the structure which appears. 



only if we do not omit eo and in the expressions for the quantized potential; here 
A(x,y) is the massless Pauli- Jordan distribution. The same happens with the two- 
point function, 

(ft, A^x)A v {y) ft) = -rj^A+ix, y), (III.6) 
which contains not only the transversal part: 

A£?(z,y) = 0, 
but the full scalar two-point function: 

= p^/f «*~>. (III.7) 

On the other hand, if the expectation values of gauge- invariant quantities are investi- 
gated, then the non-physical polarizations make no contribution. For instance, 

(E\x)E>(y)) Q = d° x d° y (A\x)A\y))n + fi£flj (^ (z)^ (y))n = fl^A^a;, y), 
even if we make use of (IIII.6J) . 



CHAPTER IV 



Parametrices of the Dirac equation on external field 

backgrounds 

This chapter contains an exposition of a general method of finding parametrices of 
partial differential operators 1 . In particular, we shall consider the scalar Klein-Gordon 
operator and the Dirac operator. 

Parametrices are distributions which, when acted upon by the wave operator, vanish 
or produce a smooth function. They are useful tools in the construction of fundamental 
solutions 2 of partial differential operators. It is therefore natural that much of this 
chapter is based upon the standard textbooks on partial differential equations [ HC661 
Fri74 . It is the aim of this chapter to investigate the structure of singular solutions of 
the Dirac equation. The two-point functions of the quantum Dirac field are examples 
of such solutions. Being distributions in two variables x and z, they are always singular 
at the coincidence point x — z. We shall present a way to deduce the singular part 
of such solutions which will be local, that is, it will depend only on the external field 
backgrounds in a neighborhood containing x and z. 

Before we begin, let us introduce the geometrical setting (figure HV.IJ) . The point z 
is fixed, it is the origin of the construction of the parametrix. The point x varies, and 
the parametrix fulfills the wave equation with respect to this point by construction. 
There is a unique geodesic curve which joins x and z, which is just the straight line. 
The tangent vector of this line is denoted by £. This vector is normalized to ±1, if x 
and z are not null-related. Moreover, the (positive-valued) geodesic distance of x and 
z is denoted by s. The square of the Lorentzian geodesic distance is denoted by 

r = (x - zf = (x - z ) 2 - (x - z) 2 . 

We have 

x a = z a + sC- 



Many formulas will be present in this chapter, some of which are easy to derive. Because of that we 
have decided to emphasize some of the more important equations with a box in order to distinguish 
them from what may serve as a reminder of their derivation. 
2 Such as retarded/advanced solutions of the d'Alembert equation. 
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(x-zf>0 



(x-zf=0 
light cone 




(x-z) 2 <0 



Figure IV. 1. Geometrical setting. The point z is fixed, and there the 
initial condition is specified. The parametrix fulfills the wave equation 
in the variable x. 

IV. 1. Scalar field case 

Before we concentrate on the Dirac field, let us look for the fundamental solutions 
of the scalar wave operator 



If the functions B a (x), c(x) exhibit little or no symmetry, then it is difficult to find the 
fundamental solutions explicitly. We shall present here a method of Hadamard which 
allows us to find approximate fundamental solutions of the wave equations. They are 
approximate in the sense that 

• they differ from the exact fundamental solutions by a residual term which is 
different iable up to an arbitrarily high order, 

• the differential operator T> x acting on the residual term of order N gives a 
distribution which vanishes as T N with r — > 0. 

IV. 1.1. Progressing wave expansion. In order to find approximate fundamen- 
tal solutions we use the method called 'the progressing wave expansion' 3 . The solutions 
of order N will be parameterized as follows: 

N 



n=0 

where U n (x) are C°° amplitudes, S n are distributions of one variable called the wave 
forms, and x is a smooth function, the phase. Any discontinuity of (p comes from the 



[□* + B a (x)d^ + c(x)} ip{x, z) = V x ip{x, z) = 0. 



(IV.l) 




3 See |Fri74j section 3.6.; [RC66 . 



distributions S n , and thus surfaces of constant phase are possible surfaces of disconti- 
nuity. 

It is clear from the outset that the progressing wave expansion introduces a con- 
siderable redundancy. Not only it is possible to describe different types of singular 
solutions of the wave equations (with different surfaces of singularity), but also the 
same solutions may be expressed in different ways . 

We use this redundancy in order to specify the concrete situation (we look for a 
singularity on the light cone T = 0) and to simplify the search for the amplitudes. In 
the construction we fix from the beginning the wave forms S n and the phase function 
and ask, whether it is possible to choose the amplitudes so that the distribution ip 
fulfills the wave equation. 

The wave operator T> has the same principal symbol as \3 X ; because of that it will 
turn out advantageous to require the wave forms to fulfill 

dS n ^ 

Moreover, if we choose the phase function to be equal to the square geodesic distance 5 
between x and some fixed basis point z, an additional equation between the wave forms 
must be fulfilled: 

rs_ 2 (r) + 2£_x(r) = o, (iv.2) 

for otherwise no choice of amplitudes can lead to a solution of the wave equation. 

If we choose a set of wave forms S n which fulfill both equations and if we find for 
this set the amplitudes U n {x) (which we will do in the subsequent sections), then the 
problem of finding ip is solved by means of a series (of progressing waves) . 

The equation (jIV.2|) prescribes the degree of homogeneity of 5*_i. Formally, the 
following distributions fulfill this requirement: 

So = - or S = 8{T). 

Irrespective of this choice and also of the choice of the regularization of T, the method of 
progressing waves will yield the same amplitude functions. For instance, the retarded 
and advanced solutions as well as the Feynman parametrix have the same progressing 
wave amplitudes and differ only by the wave forms. 



We may, for instance, change the distributions S n away from their singular support and absorb this 
change into a redefinition of the U n 's. 

5 This is only possible after a regularization of T, see the next section. 



From now on we shall concentrate on 



s (r) = l/r, 
5i(r) = inr. 

In this case the progressing wave expansion coincides with the ansatz of |DB60j which 
will be later recognized as the universal form of the singularity of the two-point function 
of the quantum Dirac field (the Hadamard form): 

(p(x, z) = + V ( X) z ) i n (r) + w (x, z), (IV.3) 

where v and w are sums of higher order wave forms and their coefficients: 

v = v (x, z) + v x (x, z)T + v 2 (x, z)T 2 + . . . , 
w = Wq(x, z) + Wi(x, z)Y + W2(x, z)T 2 + . . . . 

The amplitude functions u, v n , w n will be constructed locally (in the smallest causal 
normal neighborhood containing x and z) and will fulfill a set of (first-order differential) 
transport equations. 

IV. 1.2. Regularization of the phase function. In the method of progressing 
waves the wave forms were distributions of one variable taking as an argument the 
value of the phase function x( x )- Such a composition should be a distribution on the 
space of spacetime- valued test functions. Its action on g G C^°(M 4 ) can be obtained 
from 

(S[x],g) = J dr s ( T ) s b~ ~ x{x)} g(x) d 4 x. 

The operation is justified, only if the integration over d 4 x yields a smooth function of 
compact support in r. This point is in general non-trivial. A detailed investigation 6 
reported in |Fri74j section 2.9 assures that, as long as xi x ) is smooth and the gradient 
of x is non-vanishing, the function 

f(r) = J S[ X (x) - r]0(x) d A x 

is indeed a C£°(K.) function for all test functions G C~(1R 4 ) so that the composite 
distributions S = S[x(x)] are well-defined. The phase function x( x ) = T, however, 
possesses a vanishing gradient at the origin 7 x = z (for T = 0), and so it is not directly 



6 Notably the theorem 2.9.1. 

A similar difficulty arises in the attempts to define distributions of the radial coordinate r = |x| 
in M 3 . If x = belongs to the singular support of (the one dimensional distribution) S(x), then S 
cannot be extended to a composite distribution of the radial coordinate S l [r(x)]. Indeed, expressions 
like 5(r) have an obscure meaning. 



suitable for the method of progressing waves, although we know from simple examples 
that the fundamental solutions of the wave equations are precisely singular solutions 
whose singularities lie on the light cone T = 0. One finds them with the help of the 
Fourier transform. There the transformed fundamental solution (eg. of the d'Alembert 



The inverse of T>(p) has a nontrivial manifold of zeros which makes H(p) not locally 
integrable. Such expressions have to be regularized. There are different possibilities 
of defining the inverse, and they lead to different distributions; as we know, there are 
different fundamental solutions of the d'Alembert equation. 

In the case at hand we regularize the phase function x = T by means of a weak 
limit ie — > 0. Specifically, we define 



Both the above regularizations of x{ x -> z) = T and the distributions of them coincide, 
if the test functions with which they are integrated are not supported at the origin. 

As is the case with momentum space regularizations of 1/p 2 which have very differ- 
ent properties in the coordinate space (eg. retarded vs. advanced solution), different 
regularizations of T have different properties in momentum space. Regarded as distri- 
butions of both arguments, x and z, they have different wave front sets. 

The regularization of the distribution at the origin does not alter its homogeneity 
and therefore all regularizations of 8(T) and 1/T fulfill (jIV.2|) . All their progressing 
wave amplitudes are therefore the same. 

IV. 1.3. Construction of the parametrix; the transport equations. 

Consider the wave equation T> x (p(x, z) = together with the particular progressing 
wave expansion (taken with some regularization of T): 



operator) H(p) fulfills 



V(p)H(p) = e 



T e = (a? - z + ie) 2 - (x - z) 2 , 
T F = (x - z ) 2 - (x - z) 2 + ie. 



ip(x, z) = — ^ h v(x, z) In T + w(x, z), 

where the smooth functions v, w are given by the power series: 




(IV.4) 



oo 




n=0 



oo 




n=l 



The method of progressing waves is a way to find fundamental solutions which are 
singular. There always exists a possibility of adding a smooth solution of the wave 
equation to tp. This freedom is expressed here as a freedom to choose an arbitrary, 
smooth amplitude wq(x,z) which (as we shall see in the moment) influences all the 



w n (x, z), n > 0. 
With 



d x a T = 2(x-z) a = sU 

<9 a r<9 a r = 4r, 

QT = 8, 



we find 



oo / oo oo \ 

d aV = ~ d a T + f d a u + d a T Vn?"- 1 + In r ld a T n v n T n - 1 + ^ d a v n T n 

V 1 / 

oo 

+ {9aW n T n + d a F U WnF"- 1 ) , 

1 

where all the differentiations are executed with respect to x, the variable z being merely 
a parameter. We also find 

1 1 1 °° ,00 

U V = — (-2 d a Yd a u) + Uu- + 8v - + 8 ^r"" 1 + d a Td a v - + d a T ^ d a v n Y n - l + 

1 1 

1 00 1 00 00 



+ lnT 



8 J2 n ^r"" 1 + «9 a r ^ n 0°u n r"" 1 + 4 ^ n(n - 1) v n V n ~ l + ^ Dv n T n + 



+ d a Tj2n d a v n Y n ~ l 



i Dw n r ™ + 2n d a Td a w n r™" 1 + 4n(n - 1) u^r™" 1 + 8n i^r" -1 ] . 



Ordering the above expression so that appropriate progressing waves appear together, 
we find: 

Op = h-2d a Td a u) + hou + 4v + 2d a Td a v ) + 
1 z 1 

oo 

+ ^lnT r n [2(n + 1) d a Td a v n+1 + 4(n + l)(n + 2) v n+1 + + 
o 

oo 

+J2 T n [2 d a Yd a v n+1 + 4(2n + 3) v n+1 + 2{n + 1) d a Yd a w n+l + 4(n + l)(n + 2) w n+1 + Dw n ] , 



(we have set Wq = 0). Finally, by setting all the coefficients in front of the factors 
r n and T n In T to zero, we obtain the following system of differential equations of first 
order for u, v and v n ,w n : 



2 d a Yd a u + uB a d a Y = 0, (IV.5a) 

2 d a Td a v + v (4 + B a d a T) = -V x u, (IV.5b) 

2 d a Yd a v n + v n [4(n + 1) + B a d a T] = --V x v n . u (IV.5c) 

n 

and 



2 d a Yd a w n + w n [A(n + 1) + B a d a T] = 

- - K[4(2n + 1) + B a d a T] + 2^0^ + V x w n _ x } . (IV.5d) 



This is a system of partial differential equations. If we recognize that the contrac- 
tions of the partial derivatives on the left-hand side can be written as derivatives along 
the geodesic, Jj, 

d a Td a f = 2sCd a f = 2s f, 

as 

we realize that we are dealing with a recursive system of ordinary differential equations. 
They will be solved in the following order: first u, then v n with growing n, and finally 
w n . At each step we have an ordinary differential equation to solve, yet the whole 
system is a system of partial differential equations, because the knowledge of all the 
previous amplitudes in the neighborhood of the geodesic (there are normal derivatives 
in T>) is necessary in order to establish the equation for the next amplitude. 

Apart from the ambiguity of choosing wq, the system has a unique smooth solution 
with initial condition u(z,z) = 1. Setting the function u equal to 1 on the diagonal 
x = z and demanding the smoothness of vo,v n , w n determines these functions uniquely. 
The proof of the above statement will be apparent from the construction which follows. 
The role of smoothness as an initial condition is the following: the derivatives with 
respect to s are all multiplied by s and therefore have to vanish at s = 0, if we want 
the solutions to be smooth. The transport equations themselves taken at s = give 
us the initial conditions for v ,v n and w n in a recursive manner. 



IV. 1.4. Solution of the transport equations. 



The system of differential equations from the previous section will now be solved 
constructively. 

Equation for u(z,x): 
This equation allows us to find a smooth coefficient u(z,x), i.e. the coefficient which 
appears in front of the most singular term in the parametrix (jIV.4|) with appropriate 
initial condition u(z, z) = 1. Due to d a T = 2s£ a , we find 



L * Jo J 

where only B a is r-dependent. It is important to realize that the above solution is 
smooth in the limit s — > 0. 

Equation for Vq(z,x): 
The equation (jlV.Sbj) for t> may be integrated as soon as the function u(z, x) is already 
known. We find 



which is a smooth function if u(x, z) is smooth. 

Equations for v n (z,x): 

In order to find smooth solutions of (|IV.5c|) one proceeds recursively. Here we shall 
present the typical method of finding such solutions in a compact form. Noting that 



u 



exp — — 





B a d a T = — 2 d a Td a u 




1 du 



u 



it is possible to transform (|IV.5c|) into 




and then to obtain straightforwardly, 




The function v n (x, z) is smooth, if v n -\ is also smooth. 



Equations for w n (z,x): 
The equations (|IV.5d|) can be solved similarly. According to (|IV.5cj) the RHS of (JIV.Bdj) 



can be transformed into 



RHS rrvjzn = — ^rr 1 - 4fn - -V x w n -x. 

vr n 

Analogously to the f n -case one finds: 

w " = -^r r~sr- + —r- + 4m » ) dT - 

Remark. It is perhaps valuable to note that all the coefficients, u(x, z), v n (x, z), w n (x, z) 
are functionals of the external field. Moreover, what counts is the external field in an 
infinitesimal neighborhood of the geodesic (straight line) which connects 8 x and z. Due 
to the powers of s in the denominator, however, those functionals do not vanish in the 
limit s — > 0; they rather have a form of certain averages of the external field on the 
line x — z. 



IV.2. Dirac field 

Now we shall consider the case of a free Dirac field propagating in a given, fixed 
electromagnetic environment. Our goal is to find the parametrix exactly the way 
we did in the last section in the scalar case. By doing so, we shall also investigate 
the transport equations for the expansion coefficients in much greater detail, in part 
because the Dirac parametrix will be relied on this thesis. 

We consider the Dirac equation: 

i^ida - ieA a {x))H(x, z) - mH(x, z) = D xA c H CB (x, z) = 0, (IV.6) 

where we have emphasized the bi-spinorial character of H(x,z). We define also an 
auxiliary differential operator, 

D x = %i a d a + erA a {x) + m, (IV.7) 

and make use of a standard ansatz, namely, the parametrix of the Dirac equation 
(|IV.6J) is given by 

H AB (x,z) = D xA c </> CB (x,z), (IV.8) 

where <ft(x, z) is the (bi-spinorial) parametrix fulfilling the following second order partial 
differential equation: 

V(f) = DD(f) = [i(p- ieA) - m][i(fi- ieA) + m](j) = 



Not only the external field on the geodesic, because there are partial derivatives of A a (y) in all 
directions involved. 



modulo smooth terms ; here 

V=[U- lie A a d a + m 2 + e 2 A 2 } ■ 5 A B - |( 7 V - tV) A B d a A b (IV.9) 

is the spinorial differential operator which acts on the x-variable (and its spinorial 
index "A"). From now on the external electromagnetic field is considered to be in the 
Lorentz gauge OA = 0. Making use of the definition |[7 a , j b ] = & ab we obtain 

V=[U- 2ie A a d a + m 2 + e 2 A 2 } -6 A B -^ o abA B F ab . 

Now we introduce certain abbreviations which will make the following considerations 
much more apparent 10 : 

V = [U + B\x)d a + c(x)) ■ 5 A B - 6 -a abA B F ab (x), (IV AO) 

where 

B a {x) = -2ieA a {x), 
c(x) = m 2 + e 2 A a (x)A a (x), 

and so the Lorentz gauge condition reads: d a B a = 0. 

Clearly the operator (jlV.lOj) differs from the scalar operator considered in the 
previous section. Although the first part of it is multiplied with a spinorial delta S A B , 
i.e. effectively acts as a scalar operator , the second part (which does not contain any 
differentiation) contains a term a abA B which mixes spinorial indices. In other words, 
it is only because of the second term that we have to consider the spinorial character 
of the expansion coefficients. 

The bi-spinorial 11 parametrix will be found with the help of the progressing waves. 
The issue to be investigated throughout this section is the gauge covariance of the 
system of the transport equations. They are: 





s i a {d a - ieA a ) x u 




*) 


= o, 




-ieA a ) x v (x,z) +v 


>> 


*) 


= -^V x u(x,z), 


si a {d a -ieA a ] 


x v n (x,z) + (n + l)v n 




A 


V x v n _ x {x,z) 
An 



There will always remain a freedom of choosing an arbitrary smooth function wo^i z). 
10 The introduction of B a serves only as a link to the scalar field case of the previous sections. We 
shall not use it frequently, as it tends to obscure the issue of gauge invariance. 
11 Each spinor index is attached to a different spacetime point. 



and 



s C(d a - ieA a ) x w n (x, z) + (n + l)w n (x, z) = 

X> x 1 

s C{d a - ieA a ) x v n (x, z) + (2n + l)v n (x, z) + w n -i(x, z) . 

By writing the transport equations in this way, we have illuminated their covariance 
w.r.t. the gauge transformations, namely, if the external field is changed by 

A' a = A a + d a A, 

then the primed coefficients of the Hadamard expansion are related to the initial ones 
by 

u '(x) = e iek{x) u{x) 

and fulfill the (primed) transport equations. 

A couple of general remarks about the system of differential equations for the 
parametrix coefficients are in order. Firstly, they are supplied with only one boundary 
condition u(z,z) = 1 (times the identity matrix). All the boundary conditions for v n 
are to be derived recursively from the condition that they remain finite as s — > 0. This 
is only possible because s = is the singular point of each of those equations. Secondly, 
there is an impression that we deal with ordinary differential equations. Indeed, as the 
geodesic tangent vector £ is fixed, at least the homogeneous (i.e. LHS) part of those 
equations contains only derivatives w.r.t. the geodesic distance s. However, the RHS, 
in particular the operator T> x contains also derivatives in directions orthogonal to £. 
As a consequence the expansion coefficients v n (x,z) depend not only on the external 
field on the geodesic which joins x and z, but also on its values in an arbitrarily small 
neighborhood of this geodesic. The third important issue is the dependence of the 
coefficients on z. Up to now z was regarded as "origin" of the expansion and played 
no active role. However, we may pose the boundary value problem at a different point, 
say z, and ask to compare the expansion coefficients v(x,z) and v(x,z). We shall 
investigate this issue in a separate section (jIV.4j) . 

The equation for u A B (x, z) reads 

j-u A B = ieA a t a u A B 

with the initial condition u A B (z,z) = S A B . As the operator multiplying u A B on the 
RHS of the above equation is a scalar, the solution of the above transport equation 




must have form of a spinorial delta 5 A B multiplying a scalar, smooth function u, that 
is 

u = 5 A B u(z, x). 

We find, as in the scalar case 









u(z, x) = exp 


ie(x - z) a [ A a (y) dr 






Jo 





(IV.11) 



where y = z + r(x — z). Now the gauge covariance of this coefficient may be fully 
understood: although we see that the potential along the whole geodesic seems to 
count for u, the gauge term d a A influences the value of u only at the boundary points, 
that is only the A(x) and A(z) are important. To see this we recall that 

Cd a F(x) = ±F(x), 
for any function F(x), where r denotes the distance along the geodesic. We have 

r r d 

u'(x, z) = u(x, z) exp 



and so 

u'(x,z) = e teA(x) u(x,z)e- ieA(z \ 

which is the covariant transformation of a bi-scalar. 

The properties of the scalar distribution u are all that is needed to establish an 
equation for v Q A B : 

2 d a Td a v + v (A + B a d a T) = -V x u. 
Using the properties of u we find 



2d a Y 



9 a v A B - v Q A B \d a u 

XL 



+ Av A B = -V x u A B 



and thus 



v A b(b,z) = 



4s 



V y u(y,z)^ 
u(y,z) 



dr, 



(IV.12) 



where y(r) = z + £t. We may also write 



J B 



[x,z) = 



u(x,z) f 1 V y u(y,z)j 



u(y,z) 



dr, 



where y = z + (x — z)r. All the remaining v's can be determined recursively: 



V n b = 



U(X, Z) f 1 ^ VyV n - X A B {y,Z) 

nu(y,z) 



dr. 



(IV.13) 



Now that the formulae for the coefficients of the progressing wave expansion are 
explicitly written, we may argue, that the following theorem holds true: 



Theorem IV. 1 

If the external field is smooth then all the amplitudes of the progressing wave expansion 
are smooth functions for all their arguments. 



Before proving the theorem let us first denote 



so that 

V u 

and prove: 



;□ + B a d a ) X + {d X f + c - e -a ab F ab 



Lemma IV. 2 

If the external field A^(x) is real and smooth then the function z) is also smooth. 

PROOF. The function \ can be rewritten in the form 

1 f 1 
X(x, z) = --{x a - z a ) / B a [z + r{x - z)\ dr, 
1 Jo 

from which the conclusion is evident, as B a is proportional to the smooth A a . □ 



IV. 3. Explicit form of the singularity of the Dirac parametrix 

The parametrix of the Dirac equation (jIV.8|) is given by (p(x, z) which is acted upon 
by D x : 

Hab{x, z) = (eA+ m)(j)(x, z) + ift(fi(x, z). 

Suppose that some quantity, non-linear in the Dirac field, must be regularized in order 
to make it an operator-valued distribution. If the quantity of interest does not involve 
differentiations of the field operators, then in the point-splitting limit only finitely many 
terms of the parametrix will be important (because e.g. rinT — > 0). In what follows 
we shall find all those terms explicitly for a good purpose, namely, in the rich literature 
on QED in external fields there have been many attempts to regularize such nonlinear 
quantities; all of them use some subtractions. The singular part of the parametrix 
gives the universal singularity structure of all such subtractions. In other words: what 
does not match the singularity behavior which will be presented below certainly does 



not lead to operator valued-distributions (at least if we consider the representation of 
the CAR in some Hadamard state). 

If we let D act on the scalar field parametrix 



u 

r 



+ v in r + ln r + w ° r "' 

i o 
with the abbreviation a = ipY = ry a 2(x a — z a ), we will find 



H AB (x,z) = -— + 



au Du + avo 



oo 

+ J2 [Dv n + a(n + l)v n+1 T n \nT+ 
o 

oo 

+ [ Dw n + a(n + l)w n+1 + av n+l T n . (IV.14) 



We will characterize explicitly the part of the Dirac parametrix which is singular in 
the coincidence limit. Evidently, 



h ab 9 ( x i z ) 



au Du + avn 
\ + 

r 2 r 



Dv o + av i 



lnT. 



(IV.15) 



In what follows we will deal with the coefficient of the lnT-term, which is the least 
singular one. Its computation, however, requires the knowledge of v 1 which is the third 
coefficient in the recursive hierarchy of the progressing wave expansion and thus, by 
calculating it, we will simultaneously illustrate the practical side of the progressing 
wave expansion. The appropriate expansions of the other coefficients can also be 
obtained, but since their derivation is less illuminating we shall not present them here. 
We recall that 



u(x, z) = exp[ie x( x , z )}, 



X 



(x-z) a ['AM dr, 
Jo 



where as usual y = z + t(x — z) is the point which lies on the geodesic from z to x at 
the distance r away from z. The action of D, 



D x =il a d a + e 1 a A a {x)+m 



, leads to 



f a = d aX = [ A a {y) dr+(x-z) b [ d y a A b {y) rdr 
Jo Jo 



d a u = ie f a , 

Du = [m + e(A-f)\, 



where from now on we will keep all the partial differentiations under the integral to be 
with respect to y. The formula for Vq(x, z), (|IV.12j) . 

U f 1 VyU , 

4 Jo u 

requires the knowledge of T>u = [□ — 2ieA ■ d + m 2 + e 2 A 2 — |ct ■ F\u. With the help 
of the Lorentz gauge condition OA = we find 

-i 



Uu 



V y u(y, z) 



te 



te 



( ./• I DA a r 2 dr 

>o 

(y - z) a I UA a r 2 dr 
Jo 

e 



-e 2 f a f a ju, 
-e 2 f a f a }u(y,z) + 



+ (m 2 + e 2 A 2 — -a ■ F) u{y, z) + 2e 2 A a f a u{y, z). 

The above expression divided by u(x, z), that is without u at the end of it, will be 
called i[)(y, z) for brevity 12 : 



z) = e 2 {2A a f a + A a A a - f a f a } + xe 


(y - zf j 


f DA a r 2 dr 




+ {m 2 - 6 -a-F). 



With that abbreviation we have 



v {y,z) 



u(y,z) 
4 



ip{w, z) ds 



(IV.16) 



with w = z + s(y — z). 

The difficult task which lies ahead is to determine the least singular term of the 
parametrix, the one which comes with the lnT. What needs to be computed is 
Dvq(x, z) and even vi(x, z). The task will be significantly simplified, as we are only 
interested in the point-splitting limit x — > z. If we expand the coefficient of In V into 
a Taylor series in (x — z), then only the zeroth term, i.e. v±(z, z) + Dvq(x, z) , 

L J x=z 

is important in the point-splitting limit 13 . In other words, if two different smooth 
functions Vi(x, z) and vi(x, z) are investigated, then the point-splitting limit vanishes, 

lim In T [vi(x, z) — Vi(x, z)] = 0, 

x—*z 

as long as Vi(z, z) = Vi(z, z) irrespective of the direction of the limit (space- ,time- or 
light-like) . 



12 Again we stress that we do not divide by a bi-spinor u A B — S A B u but rather by a scalar u. 

13 Appropriately the first five terms of the Taylor expansion of the coefficient of 1/r 2 , ljlV.15|) . are 

important in the point-splitting limit and also the first three terms of the coefficient of 1/r. 



Before we proceed further, it is appropriate to note that 

d a f c = ! (d a A c + d c A a ) tcLt + (x - z) b ! d c d a A b r 2 dr, 
Jo Jo 

Df c = [ DA C r 2 dr + (x- z) b [ d c UA b r 3 dr. 
Jo Jo 

The calculation of t>i(z, z) + Dv (x, z)\ x=z will proceed much easier, if we first give 
the coincidence limit of various expressions: 



lim/a = A a , 
\imDf a = lDA a , 



lim d a f b = \(d a A b + d b A a ) 



lim-0 = e A a A a + m 
limdV = & 2 A a d b A a + f UA b - | a cd d b F cd , 
limQ0 = e 2 (3 «9 a A 6 <9M 6 + d a A b d b A a + AA b DA b ) - ferQF, 

where all quantities on the right-hand side of the equations are taken at the point z. 
We can now evaluate the coincidence limit of Dvq(x, z) (cf. (|IV.16|) ): 



lim —4 Dvq(x, z) — (m + eA)ip + i^ a [{ie)A a %l) + \d a ip\ = mip + 



and so 



Dv (x,z) 



mtp(z) - -7 a d a i/;(z) 



As far as V\ is concerned, from (jIV.13|) we have 



4 Vi(x, z) = — u(x, z) I rdr " ' 



u(y,z) 



thus in the coincidence limit 

32 Vi(z, z) = V y 



u(y,z) I dr ip(w,z) 
o 



y=z 



A calculation shows that in the coincidence limit 



lim(D - 2ieA a d a ) ( u(x, z) ■ J ip{y, z)dr ) = e 2 A 2 {z)i){z) + \U^{z), 



and so finally 



Vi(z,z) 



32 



1 | + (2e 2 A 2 + m 2 - % ■ f) } 



IV. 4. Left/right parametrices of the Dirac operator 

In the following we make a comment on the issue, whether what we have found, 
namely, the left parametrix of the Dirac operator is also a right parametrix of some 
other differential operator. As later on the parametrix we construct will be used to 
regularize the Wick square of the Dirac field operators : ip(x)ifj(z) :, this problem is of 
great importance. Let us recall that in the case of the scalar Klein-Gordon operator 
the left parametrix is also the right parametrix of the adjoint Klein-Gordon operator. 
This means that from 

[LT + B a (x)d x a + c(x)} <p(x, z) = V x <p(x, z) = modulo C°°, 

it follows 

[LT + B a (z)d z a + c(z)] ip{x, z) = ip{x, z)V z = modulo C°°; 
here the Klein-Gordon operator is self-adjoint in the generalized sense 

(V x( p,f) = ((p,V x f). 

However, the Dirac-operator parametrix is formed from the Klein-Gordon parametrix 
with an action of the D x on the left variable of <p(x, z). The adjoint Dirac operator is 
given by 

Dj = —ip z + e^L — m, 
so that for a classical solution 4>{x) of the Dirac equation there also holds 

4>{z) Sj = o, 

where the bar denotes the Dirac conjugation. The question therefore is: 

Does the Dirac-operator parametrix H(x,z) (cf. (|iV.14J) ) fulfill the adjoint Dirac 
equation 

H(x, z)E? = modulo C°° ? 



We do not possess, at present, any answer to this important question. As a conse- 
quence we do not know whether H(x, z) should be used in section IVl. 6. 31 as a regular- 
izing distribution of : ip(x)ip(z) : or : ip(x)ip*(z) : with both quantities differing by a 
right multiplication with a 7°-matrix 14 . 

14 For undifferentiated Wick products this issue can be attacked with laborious " brute force" calcula- 
tion of the right action of on the H (x, z) defined with the help of only few terms of the Hadamard 
expansion. 



CHAPTER V 
Hadamard form 

The singularity structure of the two-point function of the Dirac quantum field on 
an external background has a central role in the theory under development. On the 
one hand it provides a selection criterion for the allowed class of states (the Hadamard 
states), on the other hand the particular form of this singularity allows for a pointwise 
multiplication of the distributions which posses it and in turn allows for the develop- 
ment of the causal perturbation theory. 

There are two possible formulations of the Hadamard property, one which uses 
the Hadamard series elaborated upon in chapter ITVl and the other which employs the 
notions of microlocal analysis 1 and characterizes the Hadamard distributions in terms 
of their wave front sets. 

The chapter which follows reflects our struggle to translate the results obtained in 
the context of the quantum field theory on curved spacetimes to the external-field case. 
In particular, we shall attempt to translate the following theorems/propositions: 

(i) The definition of the Hadamard form in terms of its wave front set and in 
terms of the Hadamard series are equivalent, theorem IV. 41 

(ii) Every two Hadamard states are locally equivalent. 

(iii) The ground state on a static background is a Hadamard state. 

(iv) If u is a Hadamard state, then also w ( =wo a t , i.e. the state composed with 
a unitary time evolution, is also a Hadamard state. 

All the above theorems provide the crucial intuition as to which states of the Dirac 
field are allowed in our investigations. Indeed this class is much broader than anything 
considered before in the context of the external-field QED. 

Before we proceed further, let us note the connection between the two-point func- 
tion of a state of the CAR on a certain Cauchy surface and the two-point function of 
a state of the CAR on the spacetime as a whole. The latter algebra is the algebra of 
symbols ip(f), if)*(f) and their polynomials, where / £ Cq°(R 4 ) 4 are complex, smooth 



1 See appendix IbI 
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and rapidly decaying spinor functions. The algebraic relations are 



M/),#/)} = o, 



where S(x, y) is the unique anti-commutator distribution. This distribution solves 
the Dirac equation in the first variable and the adjoint Dirac equation in the second 
variable. If x° = t = y°, then 

S AB (t, x, t, y) = 5(x - y) 5 AB . 

We have an important 
Lemma V.l 

The two-point function cj + (x, y) of a state of the Dirac field on a Cauchy surface 
(of constant time, to) extends naturally to a two-point function at different times, 
u + (t,x, s,y), via 



where the composition means an integration with respect to z and w supplemented 
by an appropriate contraction of the spinor indices. The restriction of cu+(t, x, s, y) to 
t = t = s gives cu + (x,y). 

Proof. We shall establish the above equalities instead of saying that they follows 
trivially from the uniqueness property of the (weak) solutions of the Cauchy problem for 
the Dirac equation. Let ^(x) denote the (positive- and negative-frequency) generalized 
eigenf unctions of the Hamilton operator at time t = to- They form a complete set: 



where d\i v denotes the spectral measure of the Hamiltonian and we have omitted the 
spinor indices for brevity. If we replace in the integral the eigenf unctions ip p {?t) by 
their time-evolved versions, 



where U denotes the unitary propagator, then what results is a bi-distribution which 
solves the Dirac equation in the first variable and the adjoint equation in the second 



cu + (t,x, s,y) = S(t,x,t ,z) ocj + (z,w) o S(t ,w,s,y), 




rl> ^ U(t,to)rl>, 
^*-[l/(MoM 



variable: 

J dn P ip p (t,x)ip*(s,y) = S(t,x, s,y). 

This is the solution of the Cauchy problem, and it is unique, because the unitary 
propagator U(t,t ) is unique. We can also write 

^(*>x) = J d 3 y S(t,x,s,y)ip(s,y) 

for all solutions ip of the Dirac equation. 

The two-point function on the Cauchy surface t = to, uj + can be expressed in terms 
of the generalized eigenfunctions of the Hamiltonian: 

^ + (x,y) = J dfx p d[i k Vv(x) 5(p,k) ip* k (y). 

A composition with S(t,to) on the left side and S(t ,s) on the right side with the 
completeness relation for the ip p 's yields the desired result. □ 

V.l. Two definitions of Hadamard states 

Let a state be defined via its two-point function 

u\ B {x,z) = u{il) A {x)ijj* B {z)) = iG\ c {x,z) 7% 



in the notation of chapter IVI1 understood as a bi-distribution on the double copy of 
the spacetime as a whole, not just the Cauchy surface. 

Definition V.2 (Hadamard series definition) 

A quasifree state of the Dirac field defined with the help of the two-point function 
uj ~ab( x -> z ) i s a Hadamard state, if for each N 

is a continuous function and all its derivatives up to the order iV are also continuous. 
Here H N denotes the Hadamard parametrix (jiV.14j) with the series cut at iV and with 
the u>o-term absent. 

The other definition of Hadamard states is: 
Definition V.3 (Microlocal definition) 

A quasifree state of the Dirac field defined with the help of the two-point function 
u ~ab( x i z ) i s a Hadamard state, if the primed wave front set of uj + is 

WF\to + ) = {(x,C,y,0 -x ~ y,£ > 0}, 



where x ~ y means that the points x and y can be joined by the lightlike geodesic 
(straight line) with tangent vector £. 

Remark. Note that in the second definition the normalization of uj + is important as 
auj + for an arbitrary constant a has the same wave front set. Another important issue 
is that io\ B is a bi-distribution, i.e. a 16-component matrix. The wave front set of 
such an object is defined as the union: 

WF(uj + ) = |J WF(u+ B ). 

AB 

In the light of that, it is truly remarkable that the micro local definition is equivalent 
to the definition given in terms of the Hadamard series. The situation is similar in 
quantum field theory on a curved spacetime 2 . With the identity of Lichnerowicz 3 , 

(if + m)(-if + m) = □ - ^R + m 2 , (V.l) 

we may also write 

u~X B {x, z) = + m)y? + (x, z), 

where y? + (x, z) is a spinorial bi-distribution, the Hadamard property of which needs to 
be investigated. 

The equivalence of both definitions, essential for the further development of the 
theory, has been proved for the scalar field on a curved spacetime by M.Radzikowski 
|K.a96j . For spinor fields the appropriate modification has been given by S. Hollands 
[Hol99 and independently by K.Kratzert KaOO . The generalization to other vector- 
valued fields has subsequently been given by H.Sahlmann and R.Verch in |SV2j . where 
the authors also remove a gap present in the previously cited works. We follow the 
proof of S.Hollands in the sequel. 

Theorem V.4 (Equivalence of both definitions) 

The definition of Hadamard states in terms of their wave front sets and the definition 
with the help of the Hadamard series are equivalent. 

Proof. (Microlocal definition =>- Hadamard series) 
In the proof we shall make use of two plausible properties: 

• The squared (spinorial) Dirac operator V x , eq. (jIV.9|) . possesses four dis- 
tinguished parametrices Ar, A a, Ap, A-p, which are called the retarded, 



2 See |FV01j section 2.3. 

3 Hcre y denote a covariant derivative which contains an appropriate spin connection. 



advanced, Feynman and anti-Feynman distributions. They are distinguished 
by their wave front sets: 

WF'(A A ) ={(x,Z,y,£): x~p6 J + (y)}, 

WF'{A R ) = {(x,C,y,C) ■ x~y-xe J-(y)}, 

WF'(A F ) = {(x,£,y,£): x ~ y; [£ Q > if y G J+(x)}, [£ Q < if y e J~(x)}}, 

WF'(A ¥ ) = {(x,£,y,£): x~y, iZ >Oii y eJ-(x)},[Zo<OiiyeJ + (x)}}. 

Those parametrices are unique in the sense that, if a bi-distribution is a 
solution of the inhomogeneous equation, 

T> x A(x, y) = 5(x — y) + smooth function, 

and possesses one of the wave front sets named above, then it must be equal 
to the respective parametrix (modulo smooth function) . 
• Up to smooth functions there holds 

A F + Ap = A A + A R . (V.2) 

Both of the above properties are plausible, because they hold for the (squared) Dirac 
operator (jV.l|) on a curved spacetime and the proof of this does not depend on the 
particular structure of what other terms apart from the principal symbol appear in 
the differential operator. The principal symbol is, however, the same on background 
external fields. Note that the first property already exhibits almost the same strength 
as the theorem does. It says that in each case the singularity structure of one of the 
components already fixes (up to smooth terms) the whole 16-component matrix. 

Let us denote the (also distinguished) parametrices of the Dirac operator by Sa, 
Sr, S Fj Sp-: 

S t (x, y) = D x A$(x, y) = [n a d a + e 7 a A a (x) + m] A$(x, y) 

for all parametrices (the subscript jj denotes either one of A,R, F, F). 

With those properties we may now follow the proof, the idea being to show that 
u + can be expressed (up to smooth terms) by Sa and Sp. We define 

ujp = iu + + Sa (V.3) 



It is important that we speak of inhomogeneous solutions, otherwise, if T)A = for some singular 
bi-spinorial distribution A, then also, for instance X>(A • 7 ) = 0. 



and show that 

u F - S F = C °° function. (V.4) 

Clearly, then 

iu + = S F - S A , (V.5) 

and what appears on the right-hand side is distinguished up to smooth functions and 
therefore possesses the universal Hadamard expansion into powers of V. 

In order to show ()V.4j) we investigate the wave front set of u F . We notice that, due 

to 

u + + u~ = S = S A - S R , 

we also have 

UJp = —UJ~ + Sr. (V.6) 

For x €/J~~(y) the advanced solution S A vanishes and, therefore, for such points 
WF{uj F ) coincides with that of WF(u + ), which is known by definition (cf. IV. 3|) . 
On the other hand for x £f_J + (y) the retarded solution vanishes and from (|V.6J) and 
definition IV. 31 we also know the wave front set. Hence we conclude 

WF'(u F ) = {(x,Z,y,0: x ~ y; > if y G J + (x)], [£„ < if y e J-{x)}}, 
which is the same as the wave front set of S F 5 . An analogous property holds for 

UJp = -UJ + + Sr. 

We therefore conclude that oj ff - have the same wave front sets as S F -p. From equation 
(|V.2|) and u F + iv-p = Sr + S-j we infer 

uj f — S F = uj-p — Sp- 

modulo smooth ingredients. Thus 

WF(S F ) D WF(uj f - S F ) = WF(up-Sp) C WF(Sp). 

But the sets WF(S F ) and WF(S-p) are disjoint. Therefore u F — S F as well as u-p — Sp- 
are smooth bi-spinors which is what we intended to show. The Hadamard series is 
the necessary expansion of any two-point function uo + which fulfills the microlocal 
definition because up to smooth terms it can be expressed by the unique parametrices 
S F and Sr, the expansion of which does have precisely the Hadamard form. □ 

With the equivalence at hand we obtain an immediate corollary: 



The application of the differential operator D to the parametrix does not enlarge the wave front 
set, due to the pseudolocal property l|B.5(l . 



Theorem V.5 

Every two Hadamard states u>i, u>2 are locally equivalent. 



Proof. With theorem III . 71 we only need to show that 

|| \ABi - \/B^Wh.s. < oo, 

where Bi, B 2 are the positive operators corresponding to uix, u 2 restricted to a bounded 
open region C. With the inequality of Powers and Stormer [PS 70 , 



it suffices to show that the trace norm \\B\ — B 2 \\t r taken on C is finite. This, however, 
follows trivially, because, by definition IV. 21 the difference B\ — B 2 is an operator with 



Suppose that the external field was static for some time interval and that the state 
of the Dirac field in this region coincided with the ground state. It is very plausible that 
such a state is a Hadamard state (see the next section). An important question arises 
whether the state remains Hadamard, if we switch on in a smooth way some time- 
dependent external field. The answer to this question is affirmative. More precisely, 
the following theorem has been proven |SV2j in a (much more complex) context of 
curved spacetimes: 

Theorem V.6 

Let u>(x, y) be a bi-solution (modulo smooth function) for the wave operator T> . More- 
over, assume that there is a causal normal neighborhood iV of a Cauchy surface S f 
so that <jj is of Hadamard form for the wave operator on N. Then, if N' is a causal 
normal neighborhood of any other Cauchy surface, E s , u is also of Hadamard form for 
the wave operator T> on N' . 

V.3. Static ground states and the Hadamard form 

In what follows we shall attempt to answer the question as to whether ground states 
in time-independent external fields are Hadamard. 

Conjecture V.7 

Suppose that the external field is time-independent and smooth. Then the two-point 
functions of the ground state of the Dirac field are of Hadamard form. 




smooth integral kernel. Such operators are trace-class. 



□ 



V.2. Time evolution preserves the Hadamard form 



Remark. The only ground state is the state for which the splitting of the one-particle 
Hilbert space TC selects functions of positive energy, that is the corresponding projec- 
tions (jll.7|) are 



The above operator, as a function of the self-adjoint H 0} may be defined in terms of the 
functional calculus. Indeed, in many references there appear expressions of the form 



where ^ p (x) denote the (possibly generalized) eigendistributions of H parameterized 
by p and d/j,(p) denotes a measure 6 in p, roughly 



where dfi s is the singular part corresponding to the bound states. If the integral kernel 
of the two-point function is written in that way, it may be useful in some cases, but it 
obscures very well the singularity structure of G + , in particular the leading divergent 
terms in the limit x — ► y are difficult to obtain. 

justification OF the conjecture IV. 71 We will present an attempt to prove 
this conjecture. Its statement is plausible, because the ground states (even all passive 
states) on the static external gravitational background are known to be Hadamard 
states [SVT| . 

The approach we shall present makes use of the techniques of microlocal analysis 
(see appendix |Bj) . We include it here, because we believe this may lead to a complete 
proof in the future. In overview the proof may proceed as follows: 

(i) Definition of the two-point functions of the ground states in terms of the 
causal anti-commutator bi-distribution S(x,y) and P + , lemma IVTT1 

(ii) Calculation of the wave front set of with the help of microlocal techniques. 

(iii) Equivalence of the description of Hadamard states by the Hadamard series 
or the wave front set, theorem IV. 41 

Let us begin with the first point. The two-point function is a bi-distribution which 
takes as arguments the four-component test functions, i.e. 



P ± = (l±sgn# )/2. 





dfi(p) = d 3 



'p + dfi s (p) 



iG + (x,z) = (n,il)(x)ifj(z)Q), 

iG~(x, z) = (n,^(x)^(z)n) 



'We suppress the spinorial indices for brevity. 



act on C^° 4 in both of their variables. They are also solutions of the Dirac and the 
adjoint-Dirac equations in the appropriate variables. 
Due to lemma we have 

G + (x, z) = S(x, y) o P + (y) o S(y, z). 

In the following we shall investigate the wave front set of this distribution. We shall 
begin with the investigation of P + o S. As the bi-distribution S is not properly sup- 
ported (and neither is the Schwarz kernel of P + as we will see), we should explain, why 
the composition of P + and S is well-defined. 

It is not difficult to see that the above problem causes no difficulties. For a test 
function / G C^° 4 the general composition of distributions would require S(x, f) to be 
a test function, which is not the case, because it does not decay rapidly in the time 
direction. However, P+{y) is supposed to act on the L 2 wave functions on the Cauchy 
surface t — y°. Certainly, S(x, f) is such a function (as / is smooth and S : — > C°° 
and even compactly supported at x° = y°, due to the causal support of S). Therefore, 
the composition P + o S is well-defined. 

A standard fact, which we shall not elaborate on, is that the unique anti-commutator 
distribution S has the wave front set 

WF , (S) = C = {(x,£,y,0: x~y}, 

where x ~ y means that x and y are light-like related with £ being the tangent vector 
of the geodesic (null ray) which joins them. The wave front set of P + o S will be found 
from the propagation-of-singularities theorem which states that 

WF[P + (y) o S(y, z)] = C n [^(O) \ {0} x 1], 

a being the symbol of P + . 

There are two difficulties, however: the projection operator P+{y) which may be 
written as P+(y) = 9(po) is not a \l/DO (pseudo-differential operator), because its 
"symbol" is not continuous at po = 0. The second pitfall is that the theorem we would 
like to invoke is present in the literature only for properly supported ^DO, and P + 
certainly is not properly supported 7 . Both of the above difficulties will be cured in 
what follows. 

The problem of the lack of differentiability at Pq = may be cured by means of 
a redefinition of P + is such a way that it is smoothed out in a compact vicinity of 
Pq = 0, say on the interval [— e,e]. If zero does not belong to the spectrum of the 

Its Schwarz kernel is proportional to x0 _y 0+ie ■ 



Hamilton operator (which we assume; if not true, the splitting with respect to P + 
would result in peculiarities, anyway), then the composition P + o S does not change 
as a bi-distribution, as is easily seen from its Fourier transform w.r.t time. This fact 
already settles the issue. In general, the definition of a discontinuous symbol may be 
given in terms of a convergence of smooth symbols, as indicated by exercise 6, p. 4 of 

As to the second issue, whether the conclusion of the propagation-of-singularities 
theorem remains true even for the P + , not properly supported, we may proceed in the 
following manner: we can decompose this \I/DO into a singular part, which is properly 
supported, and a regular part, which is not. The regular part possesses a C°° x C°° 
integral kernel. We have 

p+(v) = P+(.v) + p+(v)- 

At this point we stumble upon two difficulties. The first is that one would need to 
argue that P£° o S as a bi-distribution is actually a smooth function in both variables. 

If this is true the whole wave front set of P + o S comes from P£ o S. It is not 
difficult to show that the principal symbol of Pf is equal to that of P + . It follows that 
the inverse of this symbol, cr -1 (x,p), is given by 

a-\0) = {(x,p): Po >0}. 

However, we would like to invoke the propagation-of-singularities theorem IB.6| 
which would give the desired result that the wave front set of G + is: 

WF'(P + o S) = {(x, £, y, :x ~ y, £ 2 = 0, £o > 0} 

However, the second difficulty occurs, that for pseudo-differential operators which are 
non-polynomial functions of only few momenta this theorem is not valid 8 . The difficulty 
lies in the crucial decay property of its symbol, which is not fulfilled by operators like 
0(po) - see pu£[2]. 

□ 



'We are grateful to W. Junker for pointing out this important obstruction to us. 



CHAPTER VI 



Construction of local non-linear observables 

VI. 1. Causal perturbation theory - an overview 

The causal perturbation theory (CPT) is a mathematically precise way to construct 
interacting quantum field theories. In this section we shall give a brief overview of this 
method in order to put what comes in later sections in a firm context. Our exposition 
is based upon [BFOO . 

Instead of starting from the (non-linear) Maxwell-Dirac field equations, which are 
difficult to interpret as operator equations (at least at the beginning), the CPT first 
attempts to construct the evolution operator § as an operator-valued distribution with 
values in an algebra W to be defined. This operator is a functional of g and of the 
interaction Lagrangian C, symbolically given by the Dyson series 

§[£,<?] = J2^\ / d ' x ^-- d i x n g(x 1 )...g(x n )T[£{x 1 )...C(x n )}, (VI. 1) 

where 

C(x) = e : ^(x)7>(x) : A a (x) 

is the interaction Lagrangian of the quantum electrodynamics. The functions g(x) are 
the test functions from a convenient test space; they are intended to be equal to 1 in 
the region where we wish the interaction to take place 1 and to vanish outside of it. 
The series (| VI. 1)1 should be understood as a formal power series in e. The terms in the 
Dyson series 2 

T[C{x 1 )...C{x n )\ 



1 The results of the finite-order time dependent perturbation theory of ordinary quantum mechanics 
are in general trustworthy only for a limited time scope; the strong limit \t\ — > oo of the finite- 
order interacting evolution operators Ui(t,to), even if it exists, only in special cases (like stationary 
perturbations) leads to results of clear physical meaning. Thus, we shall think of g(t,x.) as having 
its support on a finite (possibly small) time interval t G 7; the compactness of the support in spatial 
dimensions will be less important, as we shall see in concrete applications. 
2 The time-ordered products can be obtained from S[£, p] by means of a functional derivative 

S n 

T[C{x x ) . . . £{x n )} = — — - §[£, g] 

i n og(xi) . . .og[x n ) 
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are called time-ordered products of n-th order. Obviously the whole information about 
the QED is contained in the time-ordered products, and their construction also contains 
the renormalization difficulties. 

The main idea of CPT is to obtain the time-ordered products by means of an 
inductive procedure, namely, one attempts to construct T n , the n-th order time-ordered 
product with the knowledge of all the time-ordered products of lower orders, T <n . This 
is possible with the assumption of causality. 

Definition VI. 1 (Causality of time-ordered products) 

The hierarchy of time-ordered products satisfies the causality requirement, if for any 
set of points, xi, . . . , x n , yi ■ ■ ■ , y m , such that the points Xj do not lie in the past of the 
points yi (cf. figure IVXTj) the factorization property holds: 

T[C{x x ) . . . C(x n )C( yi ) . . . C(y m )] = T[C( Xl ) . . . C(x n )} T[C{y x ) . . . C(y m )}. 

Here the LHS is a (m + n)-th order TOP, and the RHS is a product of n-th and m-th 
order TOP. The factorization condition is equivalent to the (more intuitive) condition 
on the evolution operators, namely 

%i + <? 2 ) = % 2 ) %i), 

whenever 

supp g 2 £J~ (supp #1 ). 




Figure VI. 1. Relative position of the point y for the factorization of 
the TOP and of the supp g\ for the factorization of 8. 

Given a set of points x\, . . . , x n we ask whether there is a partition of it such that 
at least one of the points lies in the causal past of all the others. Evidently, this is 
possible, only if the points do not all coincide 3 (x x — x-i — . . . — x n ). Therefore, if 
at least one pair of points does not coincide, it is possible to construct the n-th order 

3 The coincidence of all points will be called the diagonal of (R 4 ) xn . 



TOP. The renormalization problem, in this scheme, manifests itself as the problem of 
extending the ra-th order TOP to the diagonal. The TOP are, however, operator-valued 
distributions with values in an algebra W (to be defined later). Instead of searching for 
an extension of OVDs, one reduces the problem to the extension of c-number-valued 
distributions with the help of the Wick expansion. 

Theorem VI. 2 (Wick expansion) 

The product of Wick polynomials can be expressed as a sum of Wick polynomials. 
Without loss of generality we have 

n+1 

: 4)*{xi)^){xi) : ... : ip* (x n )?p(x n ) : • : ip*(y)^(y) : = ^ d m (xi, ...,x n ,y) : ■0'V (m) :, 

m=0 

where : ip* : is the Wick product of m-th order which contains m field operators 
and m adjoint field operators. The d m denote number-valued distributions 4 . 

Focusing attention on the TOP of the order n, T n , we may factorize it to become 
a product of lower order TOP's and, after Wick expansion we obtain the n-th order 
TOP in the form 

n+1 

T n (xi, ...,x n ) = T[C(xi) . . . C(x n )} = 2J tm(xi, ■ ■ ■ , x n ) : 4>*^ (m) :, 

where the distributions vL( known for all points with the exception of 

the diagonal x\ = ... = x m . We need to extend it to the diagonal, because the 
evolution operator S(g) smears the time-ordered product T n (x\, . . . , x n ) with the same 
test function g(x) in every variable, and hence the diagonal belongs to the support of 
the product of the test functions. 

If relative coordinates are employed, for instance xf = Xq + then the diagonal 
coincides with the origin of the (M 4 ) xn space, where = 0. There is an elegant way to 
extend distributions singular at the origin to that point, if only the distributions are 
homogeneous. If we require that the scaling of the distribution must be preserved by 
the extension, then the following possibilities arise: 

• the extension to the diagonal is unique (this is usually the case in the so-called 
tree Feynman diagrams) and may be obtained, for instance, by multiplying 
the Wick expansions for different configurations of points with an appropriate 
step (Heaviside) distributions; 



See <|VI.9() for a concrete example, where the product of two first-order Wick products has been 
expanded. 



• the extension to the diagonal is not unique (which is the case in the loop- 
containing Feynman diagrams) and is characterized by a finite number of free 
parameters, which correspond to the renormalization constants. 

The truncated evolution operator §, constructed with the help of the time-ordered 
products of order not greater than n and smeared with test functions g(x) which 
specify the interaction region, describes the quantum electrodynamics in finite-order 
perturbation theory. The results obtained in this way correspond [Sch96 to those 
obtained in the other formulations of the perturbation theory (which in our opinion 
have a less transparent structure). 



VI. 2. Algebra of Wick polynomials 

Let from now on OL{0) denote the CAR algebra of the free Dirac field. The 
region O is an open bounded region in spacetime which additionally is geodesically 
convex (a double-cone). In order to incorporate the interacting field we shall investigate 
algebras larger than (X(0). In particular, we will define the algebra W(0) which also 
contains the Wick polynomials 5 of the free fields. Before we define W(0) consider the 
Wick square of the free Dirac field. 

Example. The Wick squares are defined by the point-splitting limit: 

: ip*ij) : (x) = lim ip*(x)ip(y) — d(x,y) , (VI. 2a) 

y^x L J 

: ipip* : (x) = lim - d{x, y)} , (VI.2b) 

y^x 

with appropriate distributions d(x, y), d(x, y). The smeared Wick square is a functional 
of one test function: 

:V>V : (/) = y d 4 x d 4 y[ip*(x)ip(y) - d(x,y)} 5(x-y)f{x). 

The above formula looks as if we have smeared the ordinary operator-valued distribu- 
tion 

ij)*(x)il>(y) - d(x,y), 

(which when smeared with two test functions becomes an operator belonging to (X 
and is well-defined for all d(x,y)) with the distribution 6 5(x — y)f(x). The result is 



5 For the purposes of quantum electrodynamics it is sufficient to restrict the investigations to the Wick 
polynomials of undifferentiated field operators with an equal number of field operators ip(x) and their 
adjoints ip*{x). This would, however, be insufficient if we wanted to define the energy-momentum 
operator. 

6 This is the nomenclature frequently employed in H W2| . 



an OVD in the algebra of Wick polynomials W irrespective of d(x,y). However, in 
order to find a representation of this algebra it will be necessary to restrict the class 
of allowed d(x, y). 

Definition VI. 3 (Wick products of n-th order) 

The Wick product of n-th order is a regularization of the product of n field operators 
and n adjoint filed operators (grouped in pairs). It is defined inductively in terms 
of Wick products of order n — 1 and lower. Let d(x, y) and d(x, y) denote the two 
distributions which have the same wave front set: 

WF'(d) = {(x,Z,y,Z):x~y,Z o >0}. 

The n-th order product is given by 

: i/j*(x 1 )ip(y 1 ) • ■ -^*(x n )^(y n ) : = i/j*(x 1 )ip(y 1 ) • • • ^*{x n )^{y n )- 

-^2 l A{x k ,y i ) : i/j*(x 1 )ip(y 1 ) • ..ip*(xi)ip(yi) . . .^*(x k )^(y k ) . . . ij>* (x n )if>(y n ) : - 

k,l 

(contraction of two pairs ipip*) : (product of remaining operators) : — 

n 

— (contractions of k pairs) : product of (n — k) pairs :, 

k=3 

where A(x k , yi) denotes either d(x k , yi) in case I < k or d(x k , yi) in the other case. The 
tilde denotes the operators which do not appear (have been contracted). The Wick 
product of an equal number of field operators and their adjoints in arbitrary order 
is defined similarly, the only difference is that in front of each contraction there is a 
factor (— iy , where j is the number of permutations necessary in order to bring the 
contracted operators together. 

Remark. The definition of Wick products makes use of only two distributions d, d. It 
is assumed that they regularize the two-point functions of Hadamard states, namely, 
the distributions 

u(i) A (x)i)* B (z)) -d AB (x,z), 
u(ij* B (z)ij A (x)) -(1ba(z,x) 

are assumed to be smooth for every Hadamard state uj. The pointwise product of d 
and d is well-defined as a distribution, due to the properties of their wave front sets. 
The multiplication enlarges the wave front set; for instance, d(x, y)d(x, y) has the wave 



front set 

WF'{dd) = {(x,Z 1 ,y,Z 2 ):Z i eJ + }. 

This set is stable under further multiplications with distributions having the Hadamard 
form of sine; ularity [BFOO . 

Lemma VI. 4 

Let T be the Hilbert space constructed upon a Hadamard state u. The Wick products 
defined above, even upon multiplication with 5(xi — yj), are operator- valued distribu- 
tions on T . 

Remark. We do not prove this lemma here, although for second-order Wick products 
it is not difficult. The generalization to higher order products for the scalar field has 
been given m [BFOO . In order to use these results here we would need a characteri- 
zation of the Wick products in terms of derivatives of certain functionals which we do 
not have at the moment. 

In what follows we shall prove a lemma which greatly facilitates the Wick ordering 
of the time-ordered products. 

Lemma VI. 5 

If the regularizing distributions fulfill the condition 

d AB (x,y) + d BA (y,x) = S AB (x,y) = {i/j A (x),i/j* B (y)}, 
then the Wick products change sign upon commutation of adjacent field operators: 

: V* • • • ip*(a)ip(b) . . . ip*(x n )ip(y n ) : = 

= - : ^*(a?i)^(j/i) . . . ip(b)ip*(a) . . . i/;*(x n )ip(y n ) : . 

Proof. The proof is standard and can be completed with combinatoric methods. 
We shall consider the case when the commutation occurs right at the end of the Wick 
product 7 . The LHS in that case is 

LHS = ip* (xi)ip(yi) . . . ip*(x n )ip(y n ) ip*(a)ip{b) + contractions. 



This makes the proof more evident. The general case is more laborious. 



The RHS is similarly 



RHS = i/;*(xi)i>(yi) ■ -fW^W r(a)ij(b) + 



S(b, a) ip* (xi)ip(yi) ■ ■ ■ V ; *( a; n)V ; (z/n) + contractions 



where "contractions" means all the due terms from the expansion of 



: ^(x^iyi) . . .ip*(x n )ip(y n ) ip*{a)ip(b) : . 



Similarly, "contractions" means all the due terms from the expansion of 



(VI.3) 



Evidently, the products of field operators are equal on both sides; what remains to be 
done is to show that the terms with contractions are also equal. Let us take the RHS 
and consider the two possible cases: 

• There is a contraction of ip*{a) or ip(b) with the other points. Then in the 
Wick product there is a minus sign due to the even number of field operators 
between the contracted quantities. This minus compensates the minus of the 
whole expansion (jVI.3|) . and what remains after the contraction corresponds 
precisely to the respective term of the expansion on the LHS. 

• The operators ip(b)ip*(a) are contracted together. This leads to the term 

d(b,a) [all contractions of : ip* (xi)ip{iii) . . .ip*{x n )ip{ii n ) :], 
which is equal to 



d(b, a) [V>*(xi)V>(yi) • • • ^*{x n )il){y n )\. 



If we combine this with the other summand on the RHS (which involves 
S(b, a)) it yields 



d{a, b) [^(x^ipiyi) . . . i/;*(x n )i/;(y n )}, 



which is precisely the sum of contractions on the LHS in the case of ip{a)ip*(b) 
contracted together. 



□ 



Remark. As a simplest corollary we note that 



: ip(x)ip*(y) : 



: ip(y)ip*(x) : . 



Also 

: ip*(x)ip(x)il)*(y)ip(y) :=: ip*(y)tl;(y)ip*(x)ip(x) : . 

Both the above expressions facilitate greatly the causal Wick expansion of the second 
order TOP (see section IVI.6.2j) . 

With this preparatory definitions we may now define the larger algebra W(0) 
which will be large enough to contain the non-linear field quantities, in particular the 
interacting fields. 

Definition VI. 6 

On the distributional level the algebra W(0) is the algebra generated by the Wick 
polynomials of an even number of elements 

W 2 (x,y) = : *p*(x)iP{y) :, 
W 4 (x 1 ,x 2 ,y 1 ,y 2 ) = : ^(xiMyiWfatyfa) 

possibly multiplied with delta-distributions in order to define Wick powers. The ele- 
ments of the algebra W(0) are formed, when the resulting quantity is integrated with 
the test functions of compact support in O: 

W{0) 3 W(gx, . . . , g n ) = J d 4 X! . . . d 4 x n gi(xt) ■ ... ■ g n {x n ) W(x u ...,x n ), 

where for brevity we have omitted the remaining arguments y%, . . . , y n which should 
also be smeared with the test functions. 

Remark. In applications the expectation values of the elements of W will be studied. 
The expectation values of the Wick products W n (xi, . . . , x n ) will become harmless 
smooth functions, but their distributional coefficients c n will require investigation (see 
the expression (|VI.9|) for a concrete example). 



VI. 3. Locality in causal perturbation theory 

In this section we will introduce a new condition of locality for quantum field 
theory on external backgrounds. This new locality requirement will be stronger than 
the notion of locality usually employed in quantum field theory without external field 
backgrounds. In what follows we will define this new notion and explain why it is 
necessary to incorporate it in the construction of the quantum electrodynamics. We 



shall put less emphasis on the mathematical rigor of the definition of locality 8 in favor 
of a clarification of its physical meaning. 

The standard notion of locality in quantum field theory means that the observables 
in spatially separated regions should be commensurate. More precisely, if we focus on 
the algebras of observables Ot(0) parameterized by spacetime regions, we have the 
following 

Definition VI. 7 (Quantum field theoretical locality) 

The net of algebras &(0) is local, if each A 6 OL{0) commutes with all the elements 
of the algebras of the regions spatially related to O, that is 

[a, b} = o v a e a(o), b e a(o'), o \6. 

The notation O \0 means, that the regions O and O cannot be connected by a causal 
curve. 

As far as the standard notion of locality is concerned, there is the following, simple 
Observation VI. 8 

The Wick squares defined above are local in the field-theoretical sense (def. IVI.7)) . and 
this property is independent of the choice of the regularizing distribution d(x, y). 

Proof. We investigate the commutator 

[: Va^b :(x),: tpc^D ■ (y)], 

where the points x and y are spacelike related. If we undo the point-splitting limit and 
order the expression appropriately we obtain 

[■,■]= S B c(x\ y)ip* A (x)i/)(y') - S DA (y', x)i)* c (y)i) B {x'), 

which is independent of d(x, y) and vanishes, due to the properties of the anti-commuta- 
tor distribution. □ 

VI. 3.1. Local quantum field theory on external field backgrounds. 

At this point we may ask, whether it follows form observation IVI.8I that any choice 
of the regularizing distribution d(x, y) leads to a local theory. If the definit ion I VI . 71 is all 
we ask for, then indeed the net of Wick polynomial algebras W(0) are local. However, 



8 A precise definition, which can be found in jBFVOlj . uses the theory of categories. 



there is a stronger locality requirement that is in the spirit of general relativity. This 
requirement, called the Local Position Invariance (LPI), lies at the foundation of the 
theory of relativity; it is instructive to recall its definition 9 and the way it has been 
tested experimentally. Local position invariance states that: 

The outcome of any non- gravitational experiment is independent of where and when 
in the universe it is performed. The fundamental constants of non- gravitational physics 
should be constants in space and time. 

In experiments the LPI is measured with the help of the following procedure (see 
figure IVL2|) : two precise frequency standards (the atomic clocks) are employed. They 
are synchronized with a light signal. They follow their wordlines 10 and continuously 
send the light signals which carry the information about their states. Those signals 
are compared at a single event. The result is scrutinized against general relativistic 
predictions, namely, one calculates the geometric lengths of both wordlines. If the time 
lapse the atomic clocks have measured is proportional to the length of their wordlines, 
then indeed the non-gravitational experiments (here the quantum optical experiments) 
are independent of the position in the universe. 

Let us view the experiment described above from a different perspective and ask 
what does the result of the measurement (of the comparison) depend on? Clearly, if 
the LPI is fulfilled, then the result depends only on the gravitational field in the region 
of spacetime which contains the wordlines of the two detectors and all the geodesies 
which join them (the clocks must be synchronized and their state must be compared). 
Therefore, the LPI which at the moment is supported by strong experimental evidence 
[WiOl implies that the results of quantum optical experiments can only depend on 
the external field in the causal neighborhood containing the entire measurement setup. 

We now turn back to the definition of local quantum field theory. Recall that the 
quantum field theory consists of three main ingredients: 

• an algebra of observables W(0), on which 

• the time evolution acts as a group of automorphisms; 

• the states are functional on W(O); they describe the expectation values of 
the observables. 

In this paper we have adapted the following way to incorporate the LPI into quan- 
tum field theory. 



The discussion of the equivalence principle (which has various ingredients) as well as its current 
experimental status can be found in an excellent review article by C.Will [WiOlj . 
10 We can assume them to be geodesic in order to rule out the influence of acceleration on the clocks. 




light cone 



Figure VI. 2. GPS-like test of local position invariance. The experi- 
mental result is allowed to depend only on the gravitational field in the 
causal normal neighborhood containing the whole measurement appara- 
tus (together with the final signal read-out). 

Definition VI. 9 (Strong locality) 

The quantum field theory is local in a strong sense (LPI) if: 

• The bi-distribution employed in the definition of the Wick square (d(x, y)) and 
of all Wick polynomials (see Wick expansion, theorem IVI.2|) . is a functional 
of the external field (and its derivatives). In particular 

d(x,y) = d[A,dA,d 2 A,...}. 

• The distributional coefficients of the Wick powers, c n (xi, . . . ,x n ), are also 
functionals of the external field. 

• All the above distributions depend functionally only on the external field 
in the causal normal neighborhood which contains all of their arguments. 
Specifically, if g(x) denotes an arbitrary test function of compact support, 
then the functional derivatives of (d,g®g) and (c n , g® n ), with respect to the 
external field A(z), must necessarily vanish: 

8(d,g®g) *(cn,g® n ) 
5A(z) ' 5A(z) 

if the point z does not belong to the smallest causal normal neighborhood 
containing the support of g (see figure IVI~3|) . 




Remark. As a consequence of the strong locality requirement we may compare two 
quantum field theories on different external fields, say A and A. If those external 
fields match in a geodesically convex neighborhood 0, 11 the Wick algebras Wa{0) 
and W^(0) may be identified (are naturally isomorphic). This leads to the picture of 
a local quantum field theory as a functor between the categories of all possible external 
field configurations and the category of unital C*-algebras 12 . In this framework it is 
possible to investigate the change of algebraic elements under the compactly supported 
variation of the external field (relative Cauchy evolution) BFV01J. It is important 
to stress, however, that the notion of locality takes into account only the two first 
ingredients of a quantum field theory: the observables and their time evolution. There 
is an intrinsic non-locality in the construction of states of the quantum field (this is 
discussed in the next section). We should also note the new insight into the LPI 
of general relativity which is gained from the quantum field theory on external field 
backgrounds, namely, it is not true that there exist "non-gravitational experiments" 
at all. Even the local influence of the external field may lead to the alternation of 
non-gravitational phenomena such as the atomic dynamics. This together with the 
non-locality of the two point function of the radiation field, which in the view of the 
appendix [D] influences the atomic emission process, is a possible way to try to explain 
the decrease of the fine structure constant at distant epochs recently observed with the 
help of very modern methods by J.K.Webb and collaborators [Web02j. 



The causal normal neighborhood of such a region coincides with itself. Typically, regions of this 
sort are diamonds of the form O = J + (y) H J~{x) for some x € J + {y). 
The morphisms of the respective categories are explained in jBFVOlj . 



VI. 4. Non-locality of the two-point functions 

In this section we show that the two-point functions G±(x,y) contain information 
about the external field not only in the causal neighborhood of x and y, and thus they 
do not fulfill the strong locality requirement IVI.9I 

We construct our argument for the massless, hermitian, scalar field, as in that case 
the argument is transparent and simple. Let the external potential be static and only 
dependent on one spatial variable, say x 3 , called "x" in the sequel. Investigate the 
classical scattering, stationary wave functions in order to define the field operator with 
their help. 




Figure VI. 4. Two independent solutions e\(x) and C2{x). The other 
dimensions (t, xj_) are suppressed. 



The x\ and x 2 (called "xi" for brevity) dependence factorizes: 

w(t,x) = e iE(p ^ p)t ~ ip ^ -e(p,x). 

As there is no external field for small x, both independent solutions ei^{x,p) will be 
linear combinations of the running waves: outgoing e~ ipx and incoming e ipx . We may 
take as the first solution the one for which e\ = e~ ipx (which is a purely outgoing wave) 
for small x. The second solution, orthogonal to ei, will have the form 

e 2 = ae~ ipx + be ipx , 

where a / 0, because there is a potential scattering, 6^0 because it needs to be 
orthogonal to W\. Now the field can be quantized in a standard manner. We introduce 
the object (later, the field-operator-valued-distribution) 

^(t,x,x ± )= [ d 2 P± J* P ■{e lEt - ip ^[e l (x,p)a 1 (p) + e 2 (x,p)a 2 (p)]+H.c.} , 



with the CCR for the a's 



Mp), al(k)] = 5(p - k), Mp), a*(k)] = 0, 

[a 2 (p),a*(k)] = «J(p-k), [a.(p),a.(k)] = 0. 

Note that we have two sets of creation/annihilation operators, because we have sepa- 
rated the positive and negative values of p = p$. 

The vacuum representation is defined by requiring 

a.(p)fl = o. 

Now the two point function (O, ^(x)^(y) Q) (for equal times) is given by 

/ d 2 p± f°% •{e-^ ( ^[ei(x3)e 1 ( 2 /3) + e 2 (x 3 )e 2 ( 2/ 3)]}. 
Jm. 2 Jo ltj 

The square bracket is evaluated to be 

[...] = (1 + \a\ 2 )e- ip{x - y) + \b\ 2 e ip{x - y) + ab e ip(x+y) + ab e - ip{x+v) . 
On the other hand in the absence of the external field we have 

[ 1 —ip(x-y) I ip(x—y) 

Both distributions are not equal. □ 

We infer from the above example that the two-point function contains non-local 
information. We cannot implement the strong locality requirement IVI.9l bv employing 
objects constructed with the help of the two-point functions. 

VI. 5. Non-localities of the extensions of local distributions 

The construction of interacting field theory, as outlined in section IVl.H requires an 
extension of the time-ordered products to the coincidence point, where the appropriate 
c-number-valued distributions are singular. In this section, by means of an example, 
we discuss the difficulties arising in the extension process. 

Suppose that the appropriate distribution t°(x) away from the diagonal already 
fulfills the strong locality requirement 13 . Moreover, suppose that the singularity of the 
distribution t° is such that its action is well-defined for all test functions that vanish 
at the origin 14 x = 0. The following simple method of extension of t° to the origin is 
generically employed in the usual causal perturbation theory: define the subtraction 

13 For the sake of simplicity, we investigate here distributions of one argument. The variable x in 
concrete applications has the meaning of the relative coordinate of two points of the time-ordered 
product. 

14 That means that the derivatives of the test functions do not necessarily have to vanish there. 



to on the space of test functions 



(tof)(x) = f(x)-w(x)f(0) V/eCo°°(R) 



where w(x) is an arbitrary C^°-function equal to 1 in the vicinity of x — 0. The 
extended distribution (denoted by t) is defined by 



Evidently, the extension t coincides with the distribution t°, if the test function vanishes 
at the origin. Moreover, the difference of two different extensions, characterized by w\ 
and w 2 , is given by 



which is a distribution supported only at the origin, and therefore (which is a general 
property - see HC66 ) can be expressed as a sum of the delta-distribution and its 
derivatives. If we additionally require the extension to have the same scaling behavior 15 
as to, we can say which derivatives of S appear. In the case of t°(x) = 1/x the result is 

t 2 (x) = t%(x) + c S(x) 

with an arbitrary constant c. For instance, the Cauchy principal value: 



and the ie prescription 

t° = — *— (VI.5) 

x + ie 

are precisely two different extensions of the distribution t° = 1/x which differ by a 
c • 6(x) distribution (with c = —in). 

Let us consider the difficulties with the strong locality requirement (definition |VL9j) 
which arise in the presence of external fields. Suppose the distribution t° was local in 
the strong sense and had the same singularity structure as 1/x. An example of such a 
distribution is 



which shares some similarity with the first term of the parametrix of the Klein-Gordon 
operator on an external background A(x). The extension of this distribution by means 



(t,f) = (t°,wf). 



(VIA) 



(ti-h,f) = (t ,w 2 -w 1 )f(0), 




15 A discussion of the extension of distributions (in the context of interacting field theories on back- 
ground manifolds) to a point can be found in |BF00| section 5.2, where the construction of the 
extended distribution with appropriate scaling is carried out explicitly. 



of the ro subtraction is 

(t, f) = (t°, tof) = Jdxds ^A(sx) [f(x) - w(x)f(0)}. (VI.6) 

The conflict with the strong locality principle is now evident. The general extension 
depends not only on the external field in the support of f(x), which would be the 
manifestation of strong locality, but also on the external field in the region of support 
of w(x), which spoils the locality (if suppw c/supp/). 

A further investigation is therefore needed in order to make the renormalization 
method of Epstein and Glaser DM7 5, BFOO strongly local. 



VI.6. Local causal perturbation theory in the lowest orders 

In this section we shall construct local Wick and time-ordered products of first 
and second order. They are the building blocks of the evolution operator § (when 
smeared out with test functions), and it is essential to investigate them carefully, as 
any alternation of their form can potentially lead to the modification of measurable 
quantities. The second order time-ordered product already describes a number of 
important phenomena. 

We begin this section with some considerations of the lowest two orders of Wick 
and time-ordered products of the standard approach to the external field quantum 
electrodynamics. They are not local in the strong sense, as has been explained before. 
We present them nonetheless, because this discussion will be of help in the construction 
of the local Wick and time-ordered products, and because it may help the reader to 
recognize the more familiar objects and realize where the necessary modification will 
be done, in order to make these objects local. 

In quantum electrodynamics the first order time-ordered product is given by the 
interaction term Cf. 

T x {x) = Cj(x) = le : ^ A (x)^ AB ^ B (x) : A»(x). 

The second order time-ordered product satisfies the causality relation 

T 2 (x,y) = T 1 (x)T 1 (y), 

if only y ^J + (x). Thus, for x 7^ y the fermionic part of the second order TOP can be 
brought into the form: 

2 

T 2 (x,y) = J2 tn(x,y) : « (n) : 

n=0 



with the help of the Wick expansion I VI. 21 The Wick products on the RHS are inde- 
pendent of the relative position of the points, but the distributions t n (x, y) are at first 
defined only for non- coinciding points. 



VI. 6.1. Usual Wick product. Now we will recall some facts about the Wick 
product of two Dirac field operators. Let there be a fixed representation of the free 
CAR algebra with the base state described by a vector Q in the Hilbert space. We 
introduce the commonly employed abbreviations 

iG + AB (x,y) = (n,ip A (x)ip B (y)n), 
iG~ AB (x,y) = (Q,ip A (x)ip B (y)Q). 

The (non-local) Wick products are defined by 

: xp A (x)ip B (y) : = ip A (x)ip B (y) - iG\ B (x, y), (VI.7a) 
: ^a( x )^b(v) ■ = ip A (x)ip B (y) -iG AB (x,y). (VI. 7b) 

Note that the Wick products are not linear and the definition (jVI.7b|) cannot be 
derived from (jVI.7aj) with the help of the anti-commutation relations 16 . In the case of 
absent external field the distributions G ± {x,y) are well-known: 

>C-^,y) = ^/;gtf- m )e--*-»> 

We may write 

iG\ B (x,y) = (ip x + m) AB ■ A + (x,y), 
iG AB (x, y) = {ip x - m) B A ■ A + (x, y), 



"Indeed, consider 

: i/j A (x)tp B (y) : + : ^ B {y)ip A (x) := {ip A (x) ^ B {y)} - iG\ B {x,y) - iG BA (y,x) = 0, 

which follows from the definition of the Wick products. On the other hand linearity would give 

. . . =: {ipA[x),ip B (y)} ■-■ iG A b{x, y) :^ 0. 

Therefore, these definitions are logically independent. If we take (as was done so far) G~ and G + as 
the expectation values of the appropriate field operators, then they are connected, as both of them 
are in one-to-one correspondence with their projections, which however sum up to an identity on Ti., 
that is P + + P- =1. It follows from l|II.8(l that we have (in case f,g € H are test functions) 

G+(f,g) - (n,a(P + f)a*(P + g)Sl) - (f,P+g), 

G~(f,g) = (n,b(P_f)b*(P„g)n) = (f,P-g). 



where A + (x,y) denotes the massive (m), scalar, vacuum two-point function. 17 It is 
important to stress that both G + and G~ are derivatives of the same distribution. 
Consequently, the wave front set of A + contains the wave front sets of G ± . 

The anti-commutator distribution G AB (x,y) can be expressed in terms of G : 

iG AB (x, y) = {ip A (x),ip B (y)} = {Mx),ip B (y)} fi ) = iG^ B (x, y) + iG BA (y, x), 

(VI.8) 

which assures the validity of the lemma on the (anti) commutation of Wick products 
(lemma IVI.5|) . We also note that 

iG AB (x,y) = {i/> A (x),il>c(y)} 7° C b = S AC {x,yh 0C B , 
and we find, in accordance with lC.4| the equal-time anti-commutation relation 

iG AB (x, y) = i AB 5(x - y). 
The Wick product of the free electromagnetic fields is defined by 
: Af,(x)A v (y) : = Ap(x)A v (y) - r^Ajfoy), 
where Aq (x, y) denotes the massless, scalar field two-point function (|III.7|) . 



VI. 6. 2. Second-order time-ordered product. Turning now towards the sec- 
ond order TOP, T 2 (x,y), in the case y jtJ + (x), we will use the causal Wick expansion 
to order of the expression 

T 2 (x,y) = (ie) 2 : ^ A (x)^ AB M^) tfM^Mv) ■ M*)Mv)- 

As T\ already involved a limiting procedure (in the case of the Dirac field operators), 
in order to Wick-expand T 2 we will step back this limit and consider 

T 2 (x,y) = (ie) 2 [: ^ A (x)M*') Mv)My') :] Y AB Y° D M*)Mv)> 



As the free (from the external field) case is extensively studied in |Sch96) , it is perhaps valuable to 
establish a correspondence with the notation used there. For the fermionic operators we have 

S i a t\x~y) = G+ B (x,y), 

whereas in the electromagnetic case 

iD ( +\x-y) = A+(x,y) 
(the expressions on the LHS are in the notation of G.Scharf |Sch96j ). 



where only the square-bracket term requires further attention. It is seen to be equal 
to 

[...]= ip A (x)<ip B (x')ip c (y)ip D (y r ) - iG AB (x,x') : i/> c (v)il>(l/) '■ ~ 

- iG CD(y,y') '■ M x )Ma/) : +G AB (x,x')G CD (y,y'). 

The first term should now be Wick-expanded. This leads to the expression in which 
the limit x' — > x and y' —>■ y can be taken. One obtains: 



T 2 (x,y) = -e 2 [: i) A (x)i) B (x)i) c (y)i) D (y) : +iG AD (x,y) : i/; B (x)i/; c (y) : + 
+iG+ c (x,y) : $a(x)Mv) ■ -G AD (x,y)G + BC (x,y)} ^ AB ^ CD 

[. A^x)AM : +r^A+(;r, y)] . (VI.9) 



The pointwise product of G AD (x,y)G BC (x,y) is well-defined, and its wave front set is 
not greater than that of A + (x, y) 2 . 

In the other case, when x ^J + (y), there is 

T 2 (x,y) = ri(y)7i(x), 

and the similar Wick expansion can be performed. We only have to exchange x with 
y (and also A with C, B with D and \i with v). The causal Wick expansion requires 
,however, to order T\(y)T\(x) in such a way that the Wick products in T 2 (x, y) are the 
same as in (|VI.9|) . In order to do that we employ the (anti) commutation lemma IVl.51 
We arrive at 

Tx{y)T x {x) = -e 2 [: xp A (x)ip B (x)ip c (y)ip D (y) : -iG^ A (y,x) : ip B (x)ip c (y) : + 
-iG CB (y,x) : $ A (x)Mv) ■ -G + DA {y,x)G~ CB {y,x)] Y AB Y CD 

[: Ap{x)Mv) ■ +V^A+{y,x)] . (VI.10) 

A comparison of the above expression with ()VI.9j) is presented in the table 1. 

What remains to be done in order to define the T 2 (x,y) for all its arguments is 
to extend the distributional factors of Wick products to the diagonal x = y. Thus, 
in accordance with the general procedure (see section IVl.l|) . the definition of T 2 for 
all x, y reduces to the problem of the extension of number- valued distributions, which 
are known for x 7^ y, to the diagonal x = y. Some of the distributions appearing in 



y iJ + {x) 


x £J + (y) 


G AD {x,y) 


-Gl A (y,x) 


Gh c [x,y) 


-G CB (y> x ) 


G AD {x,y)G^ c {x,y) 


GbA\yi x )G C B\yi x ) 


Ao + (^Z/) 





Table 1. Correspondence of factors of Wick products with different 
relative position of the points a; and 

T 2 (x,y) can be extended to the diagonal uniquely; the other, however, have a family 
of extensions parameterized by what is known as the renormalization constants. 

VI. 6. 3. Local definition of the Wick product. The Wick products defined 
with the help of the G ± are not local in the strong sense. We shall repair this deficiency 
in what follows. 

In order to make the Wick square an operator-valued distribution it is necessary 
only to remedy the short-distance singularities of the product of the field operators. 
We proceed therefore and give a local expression for the Wick product, the unique- 
ness of which will be analyzed in a separate section (see section |VI.7.2|) . The leading 
short-distance singularities of the expectation value of the product of the Dirac-field op- 
erators in some Hadamard state are all contained in the (local) Hadamard parametrix 
H A s(x,y) constructed in chapter ITVl We therefore define 

: ^A(x)^ B (y) := Mx)^ B (y) ~ H A B(x,y) 

and, similarly, 

: tPa( x )$b(v) ■= M x )^b(v) - H AB (x,y), 
where the Hadamard parametrix H must be separately derived 18 with the methods of 
chapter llVl for the case of the adjoint Dirac operator acting on the variable x: 

H AB (x,y)Dl A c = 0, 

with 

Dl A c = {-Cp -m + efi) A c . 

18 The parametrix H is not just the adjoint of H; although j°H^° would fulfill the appropriate 
equations, it would not have the appropriate wave front set (this can be seen even on the no-external 
field distributions G ± ). 



We therefore establish the following substitution rule for the Wick products: 

H AB (x,y) <h> iG\ B (x,y), 
H AB (x,y) «-> iGj B (x,y). 

The definition of the Wick product with the help of the Hadamard parametrix 
encounters, however, two obstacles: on the one hand the full series defining Hab{x,v) 
is convergent, only if the external fields are analytic functions 19 (which is a narrow 
case), on the other hand the full H A b{x, y) is of no practical use, as it is unknown how 
to sum up the Hadamard series. 

If the coincidence limit y — > x is of interest (that is if only powers of the field are 
of interest), then an observation of V.Moretti, |Mor01j . is helpful, namely: 

Observation VI. 10 

The Wick product : ip A (x)il) 'b{ x ) '■ can be defined with finitely many terms of the 
Hadamard expansion (in this case all what is needed are u, Vq and V\). The so-defined 
OVD fulfills the wave equation (the Dirac equation) in its variable modulo a constant: 

Y\mD x : i[) A (x)il) B (y) := 3V x v (x)\ x=y 
(with abbreviations of chapter HVj note that D x and T> x are different). 

Proof. As the field operator ip A {x) fulfills the field equation by definition, what 
remains to be done is the inspection of 

lim D x H N (x,y), 

x^y 

where H N (x,y) is the Hadamard parametrix without the smooth term w(x,y) and cut 
at Vjsr. Not surprisingly the calculation which preceded the derivation of the transport 
equations (JlV.5dj) will be helpful, for 

lim D x H N (x,y) = lim V x (j) N (x,y). 

x^ry x—ty 

By inspection, we conclude that T> x (j) (x, y) consists of terms proportional to 1/r 2 , 1 /T, 
In T as well as T n and T n In T. All three first, singular terms will vanish as long as iV > 1 
due to the transport equations fulfilled by u,v and V\. The coefficients of the powers 
of T will have no effect in the limit x — > y. What will remain, however, is the coefficient 

19 The coefficients vn(x,z) are functionals of the external field dependent on the partial derivatives 
of it up to order 2 (TV + 1) (and therefore require an appropriately high differentiability of A(z)). 
What F.G.Friedlander claims in Fri74 chapter 4.3 is that also smoothness of the external field is not 
sufficient for the convergence of the series v n (x, z)r™8(r) (where 0(f) is the analogue of lnT). 



of r° which defined the transport equation for Wi(x,y) set to zero in our case. One 
finds 



lim V x (f> N (x, y) = lim V x (f)\x, y) = -Vx(12 + B a d a T) - 2d a Td aVl = -4v x + V x v Q , 



x^y 



x^y 



as a consequence of (jlV.Bcj) . Now, from the recursion relation, we have 



u(w,y) 



w=y+r(x—y) 



in the limit x — > y 



vi{y,y) = -^xM x ^y)\x= y - 



□ 



VI. 7. Local definition of the current operator; back-reaction effects 

In this section we formulate the problem of the back reaction of the quantum fields 
on the background arena. Let us recall the field equations of the theory developed in 
this paper (cf. chapter HJ): 

[i^dl + eA c J ass (x) - m] ijj{x) = -eA a (x)^(x), (Vl.lla) 
d b x F£(x) = 4vr #r) 7( ^(x), (Vl.llb) 
d b x F* ass (x)=4nJ a (x), (VI.llc) 

where A class = A is the classical electromagnetic field which is the background field, 
A is the quantum electromagnetic field, ip is the quantum Dirac field and J is the 
external classical current which is the source of the A-field. The interaction of the 
quantum fields (the RHS of the first two equations) is treated in a perturbative way. 

The main assumption of the theory is that the external current J is not influenced 
by the state of the quantized radiation field A. We propose the following way to assess 
this influence: if the current produced by the quantized Dirac field is comparable 
(locally) to J, then the external field approximation is not reliable anymore. 

Suppose that we focus on the current density J°(x) and the current-density operator 
p(x) (to be defined later). There are three issues arising, if we want to make the above 
criterion more precise. Firstly, for a given state \F) of the Dirac field it is insufficient 
to compare the expectation value (F,p(x)F) to J°(x). The current-density operator 



is not positive, and therefore one should also investigate higher momenta of it, for 
instance 20 (F, p(x)p(x)F). 

The second issue is that we have only the statistical interpretation of quantum 
theory, and thus we do not really know what happens in a single experiment. This 
is related to the problem of measurement. If the results of quantum theory such as 
(p 2 (x)) do correspond only to the average over many experimental realizations, why 
then should in a single experiment the back-reaction effect be governed by such an 
average? Let us put this question in a physical context: Imagine a weakly localized 
electron which, upon an interaction with a measurement device, for instance a CCD 
camera, produces a sharply localized, strong effect - a macroscopically visible trace. 
The CCD camera is composed of various charges and could have been regarded, for 
instance, as the source of the classical current density J°. Therefore, even arbitrarily 
low current densities can produce intensive back-reaction effects. This shows that the 
phenomenological criterion we have put forward in order to judge the possible back- 
reaction effects can fail for some sources of the external current J. 

In quantum field theory there is an additional difficulty The current- density ope- 
rator has to be defined as a coincidence limit of a difference between the two-point 
function of the state of the quantum field and a certain regularizing bi-distribution. 
This bi-distribution in the no-external-field case is taken simply to be the vacuum ex- 
pectation value of the appropriate (bi-local) current-density operator. However, as we 
have discussed in the chapter |H] there is no privileged state of the Dirac field on exter- 
nal field backgrounds. How should the current-operator be defined? 21 It is precisely 
the strong locality criterion (definition IVI.9|) that allows for a definition of the abso- 
lute current density. In other words, the strong locality criterion allows us to separate 
the (state-independent) singularity of the two-point function from the state-dependent 
information which may cause the back reaction. 

20 All operators appearing here are understood as follows: we choose a weighing function /, 



in such a way that it is supported in a small neighborhood of the point x. Then the current density 
and its fluctuation are defined as 



The need to investigate higher momenta of the charge distribution is clearly visible in the following 
example: Suppose that we want to examine the electric field of a photon. The expectation value 
of the electric-field operator (which is also not a positive operator) in the one-photon state is zero; 
on the other hand, the expectation value of the square of the electric-field operator is non-zero, and 
precisely this quantity produces an effect on a photodiode. 

21 We are not interested here in a relative current density between two states, but rather in the absolute 
current density of a given state. 



fix) > 0, 




(F,p(x)F) = (F,p(f)F) 



(F,p(x)p(x)F) = (F,p(f)p(f)F). 



What follows in this section deals with the definition of the current operator j^{x). 
We shall define it for the free Dirac field by the point-splitting procedure. The quantum 
version of j M (x) should correspond to the classical current 

J class = 6 ^(aOWfa)- 

It should also posses the following properties: 

• In the time- independent case, the integration of j°(x) over a Cauchy surface, 

Q = [ rf 3 xj°(x), 



should correspond to the total charge operator Q = e(N + — N~), where 
are electron-/positron-number operators. 
• The current operator should be conserved 22 

d„f{x) = 0. 

We recall that the classical current is conserved as a consequence of the Dirac 
equation and its adjoint, namely the (co)-spinors ip(x), ip(x) fulfill 

[ip+ e$x) — m]ip(x) = 0, 
ip(x) [—i ft + eft — m] = 0, 

from which conservation of the current follows immediately. 

As a first candidate for the current- density operator we investigate 

f{x) =: ^(x)-fTP(x) : u = Km^M-f^x) - iG~ AB {y,x)r A % (VI.12) 

where iG AB (y, x) = (fl, ip{y)%l){x) Q). The above operator is not local in that it depends 
on the external field in the regions other than the infinitesimal neighborhood of x\ 
nonetheless we shall consider its features. Let us look at the total charge. We find 



d 3 x j°(x) = J d 3 x J dfi p dfik [u\x, p)u(x, k)a*(p)a(k) — v\x, p)v(x, k)6*(k)6(p)] + 
+ / d 3 x / dfi p dfi k [u^x, p)v(x, k)a*(p)6*(k) +v\x, p)u(x, k)6(k)a(p)] . 



With regard to the current conservation wc note that the energy-momentum operator is covariantly 
conserved in the QFT on a curved spacetime as a consequence of the diffcomorphism invariance (see 
BFV01 theorems 4.2 and 4.3); the conservation of the quantized current should follow from the 
gauge invariance of the theory. 



Now, due to 

(u(k)Xp)) = <Hk- P ), 

(u(k),t;(p))=0, 



(«(k),t;(p)) = tf(k-p), 



we obtain 

| ci 3 a; = | d// p [o*(p)o(p) - 6*(p)6(p)] , 
which is just the total current Q = e(N + — N~). 

VI. 7.1. Charge conservation, local definition of the current density. The 

non-locality of (|VI.12|) was not the only shortcoming of this expression. It is easy to 
see that the so-defined current density is not conserved. Indeed, as the differentia- 
tion is performed before the coincidence limit y — > x is taken, the -0-terms are not 
differentiated and, consequently, the current is not conserved. 

With regard to the point-splitting procedure of defining j M there are therefore two 
issues: 

• The conservation of the current density. 

• The appropriate behavior of the coincidence limit. 

Various procedures have been devised in the literature to remedy both of the above 
problems 23 . The main theme of those procedures is: 

• to remedy the non-conservation by means of the symmetrization of j^{x,y) 
w.r.t. its variables; 

• to assure the finiteness of the x — > y limit by a subtraction of a certain local 
bi-distribution. 

A current density which fulfills all our requirements and which is constructed as indi- 
cated above is: 

f(x) = hm^h^Or) - n AB (y,x)Y AB + Hx)Y^(y) - n A B(x,y)Y AB ]. (vi.13) 

VI. 7. 2. Uniqueness of the current operator. The uniqueness of the current 
operator is connected with the uniqueness of the Wick square of the field operators. 
Suppose that we alter the definition of the Wick product by terms of lower order in 



23 See for instance |DM75| . where the ideas originally due to J.Schwinger together with their possible 
generalizations are discussed. 



the field operators: 

: ipip := ijjijj — d(x, y)7° + c(x)tp + d(x)?p* + f(x, y)l. 

If we require the current to be selfadjoint and to fulfill the standard locality requirement 
IVI.71 then we must have c = = d. In order to further restrict the remaining ambiguity 
we use the locality requirement, namely, the function f(x,y) can depend functionally 
only on the external field in the smallest causal neighborhood of x and y. In the point 
splitting limit this should become a functional of the external field at the point x. Such 
functionals are known to be functions of the external field A(x) and its derivatives 
d x A(x). If we further require these functions to be analytic in the external field and 
the mass m of the Dirac field, we can expand f(x,x) as a power series in A, d n A for 
all n, and in the mass m. Furthermore, the Wick product must have a dimension of 
1/cm 3 and both the external field A and the mass m have dimension 1 /cm. Due to 
this and the Lorentz condition OA = 0, we obtain the result 

Lemma VI. 11 

The local current- density operator is unique up four arbitrary numbers a% . . . a±. Any 
two local definitions differ at most by 

A : f : (x) = ai A u d v A* i + a 2 m 2 A^ + a 3 A V A V A» + a 4 UA" 



VI. 8. Scaling transformations for local observables 

In this section the (re)scaled non-linear observables will be investigated. The pur- 
pose of the scaling is to approach the short scale-limit of the theory. 

VI. 8.1. Scaling transformation in the scalar case. In the following we shall 
investigate the behavior of local observables costructed out of the scalar massive fields 
under scaling. 

In analogy to the extraction of the long-range part of the electromagnetic field 
introduced by Gervais and Zwanziger |GZ80j . 



A l ™ 9 (x) = lim XA a (Xx), 

A— >oo 

we shall associate with a (classical) scalar massive field tp(x) a field <fi(x) rescaled by 
A: 

<p(x) = Xip(Xx). 
We define the scaling transformation for classical quantities: 



Definition VI. 12 (classical rescaling) 

A scaling transformation a is a prescription to obtain various quantities at a different 
scale A, where A — > corresponds to the short-distance limit. In order to investigate 
the rescaled theory we 

• rescale all quantities, that is 

a : f(x,y,...) -> f(Xx, Xy,...); 

• interpret the result in terms of the rescaled field 

<p(x) = Xip(Xx). 

As an example, let us find the field equation fulfilled by (p. We rescale 
= a [(n x + m 2 )^{x)] = (X~ 2 n x + m 2 )^(Xx) = X~ 3 {D X + X 2 m 2 )<p{x) = 0, 
and therefore (p{x) is a massive scalar field with the smaller mass 24 Am. 

Now we shall define the scaling of the (massive, hermitean, scalar) quantum field 

Definition VI. 13 (quantum rescaling) 

The scaling transformation is a prescription which transforms operator- valued distribu- 
tions constructed with the help of the (linear) quantum field cp(x) on the GNS Hilbert 
space of the state uj into an operator-valued distributions of the scaled quantum field 
(p(x) on the GNS Hilbert space of the scaled state u. The two-point functions of the 
respective states are related by: 

u(x,y) = X 2 u(Xx, Xy). 

Remark. The scaling transformation relates the quantities of one quantum dynamical 
system to the other. The dynamical laws in both systems are different, as the field 
operators tp(x) and <p(x) fulfill different field equations. Both systems (on a single 
Cauchy surface) are representations of the same CCR algebra: 

[v( x )> V?(y)] = <J(x - y), [<p(x), <f(y)} = 5(x - y). 

Even on this single surface, however, the states uj and uj are not even locally equivalent; 
as we shall see, their short distance singularity structure is different. 



In the interesting case A < 1. 



As an example let us investigate the scaling of the Wick product which is given by 

: <p(x)<p(y) := <p(x)<p(y) - H(x,y), 

where H (x, y) is a certain bi-distribution, for instance 

H (x, y) = — + m 2 ln(rm 2 ). 

Let us now perform the scaling a on the Wick product. We have 

<7 A [: (p(x)(p(y) :] = ^[A 2 '(p(Xx)(p(Xy) - X 2 H(Xx, \y)] = ^[<p(x)<p(y) - H x (x, y)}, 

where we have denoted the rescaled regularizing bi-distribution by H\. 

H x (x, y) = X 2 H(Xx, Xy) = i + (Am) 2 In [T(Xm) 2 ] . 

The singularity structure of H\ clearly differs form H (as long as the field is massive). 
It is also evident that in the short distance limit A — > the less singular term in H\ 
will disappear, 

\im H x (x,y) = - = A+(x,y), 

A — >{j 1 

where Aq(x, y) is the vacuum, massless two-point function. 

VI. 9. Scaling transformation for the Dirac field in external potentials 

The action for the classical Dirac field on a background external potential A a (x) is 
given by 

S = J d 4 x | [ih ii^dai) - mc ^] - [^ 7 > e -A a \ } 

Here the action is dimensionless and the Dirac field tp(x) has a dimension of cm" 3 / 2 . 
Analogously to the scalar case, the scaling is performed by a substitution x — > Xx. 
Here, however, the rescaled field (classical and quantum) is defined via 

^( x ) = X 3/2 ip(Xx). 

Again the scaling transformation transforms operator- valued distributions on the GNS 
Hilbert space T into the operator-valued distributions on the GNS Hilbert space T 
constructed with the help of 

u(x, y) = X 3 lu(Xx, Xy). 



Remark. Similarly to the scalar-field case, we note that on a fixed Cauchy surface the 
OVD ip(x) and ip(x) fulfill the same ant i- commutation relations. However, 

• they obey a different dynamical law, and 



• the state us is not locally equivalent to uj. 

It is interesting to note 25 that if we introduce the rescaled external field A and the 
rescaled mass m, 

fa = Am, A a (x) = XA a (Xx), 

then for all operator-valued distributions F[x, y, ■ ■ ■ if>(x), A, m] there is 
a{F[x,y,...,i/j(x),A,m]} = F[x, y, . . . , $(x), A, fh\. 
An investigation of the Wick product 

: ^A(x)ip B (y) := ip A (x)ip B (y) - H AB (x,y), 

with Hab chosen to be a parametrix (with certain Wq(x, yj) dependent on m and the 
external field A a (x) as explained in chapter IVI.6| reveals 

— 1 — 

ipA(x)ip B (y) :} = ^ A {x)^ B {y) ■ -H x AB(x,y)} 

with 

H x ab(x,v) = \ 3 H AB (\x, \y), 

being the scaled Hadamard parametrix. In what follows Hx ab will be investigated in 
greater detail. 

VI. 9.1. Scaling of the Hadamard parametrix for the Dirac field. In chap- 
ter \ve have found the Hadamard parametrix Hab{x,z) of the Dirac field. It was 
given by 

Hab(x, z) = (i>y a d a + e>y a A a (x) + m) A C (j) CB (x, z), 
where <Pcb{x, z) is the Hadamard parametrix of the operator 

V=\U- 2ieA a (x)d a + m 2 + e 2 A 2 (x)] - ^a abA B F ab (x), 

which was determined with the help of the progressing wave expansion: 

/ x oo oo 

4>cb(x, z) = — ^ v cb(x, z) lnT + ^VnCBjx, z)T n \nT + y^ j w nC B(x, z)T n 

1 1 

where the spinor indices have been restored in order to emphasize the bi-spinorial 

character of the smooth coefficients u, v n ,w n . 

Remark. The term lnT cannot stand as it is, because T = {x — z) 2 has a dimension. 
This dimension can be cancelled with the only length scale available in the theory - 



'This (in part) will be the conclusion of the next subsection. 



the Compton wavelength of the electron. Therefore, everywhere In T should be under- 
stood as ln(m 2 T). It is another matter that a different scale (if present) would give a 
parametrix differing from the above by a smooth function. 

Because the parametrix (or certain parts of it ) have been used to define nonlinear 
quantities, it is important to investigate the scaling properties of the diagonal (x = z) 
values of the smooth coefficients under the transformation a. We find the following 

Lemma VI. 14 

The diagonal values v n (z, z),w n (z, z) of the coefficients of the Hadamard parametrix 
scale under the transformation a\ as: 

^M*)] = v n (z), 

where n — 0, 1, ... for the i>'s and n — 1, 2, ... for the w's. The v and w denote v and 
w with A, m replaced by A,rh. The coefficient of the highest singularity u(z,z) = 1 
by assumption for all A. 

Proof. Before considering the case of arbitrary n, let us first illustrate the reason- 
ing in case of n = and n — \. Everywhere we consider z as a fixed point of reference 
so that all the differential operators are being taken with respect to the variables other 
than z. We use the abbreviation 



u(x, z) = exp [iex(x, z)\ = exp ie(x - z) a \ A a {y)\ y=z+T{x _ z) dr 

Jo 

One easily verifies that in the limit x — > z there is 

(d aX )=A a (z), 

□x = o, 

with the last equation being an indirect consequence of the Lorentz gauge of the ex- 
ternal field. 

The coefficient Vq is found from 26 

" x V„u, 



u f V y u 

V (X,Z) = ~- / \y =z+T (x-z) dT . 

4 Jo u 



'In the following considerations, for brevity, we drop the tilde above the scalar part of u. 



Clearly, in the limit of interest the distributions under the integral become r -independent 
and can be pulled out of the integral, which then assumes the value 1. Therefore, 

v (z,z) = -~V y u(y,z)\y =z . 

Now the differentiation of u(y, z) results in terms of the sort D x x(x, z), A a (x)d a x( x i z )> 
[m 2 + e 2 A 2 (x)} - x( x i z ), as e -g- explained in and above the formula (jIV.16|) . The scaling 
means a replacement of (x, z) by (Ax, \z). We shall also make use of the rescaled mass 
and the external field 

m = Am A a (x) = \A a (\x). 

By the explicit formulas for z) and all its derivatives, it is evident that every term 
of T>u transforms under <j\ homogeneously with degree —2. Therefore, we have 

a\[v {z, z)\ = -r^v {z, z), 

where vq contains the rescaled A and m. 

Let us consider then the next coefficient, v 1(2, z), which is the prototype of a general 
situation. According to the formulas of the chapter HVl we have: 



1 



yi 



u(yi,z) u (2/0, z)V yo u(y , z) dr 



v t (x, z) = — u(x,z) \ r x dr x u {y x , z)V 
16 Jo 

where y\ = z + T\(x — z) and y = z + T (yi — z). Now the differential operator V yi 
either acts on u{y±, z) in front of the integral or can be pulled under the integral, where 
(as all the expressions there depend only on y ) a chain rule has to be invoked which 
e.g. transforms d yi into d yo ■ tq. In the first case there will be derivatives of x{yii z ) 
appearing in front of the ro-integral, and in the second the derivatives of xivo-i z ) an d of 
A a (yo) or c(yo). In the limit x — > z all those factors will become To- and Ti-independent 
(which also assures that they will be functions of the external field at z - and thus - 
local expressions) and can be pulled out of all the T-integrals. All the derivatives of x 
are easily seen to have the homogeneous scaling, namely, 

v\[d' ai ■■■d' an x(y, z ')\y=z< = ~^x(z,z). 

As each differential operator D contains two such differentiations (or a factor with 
identical scaling), it follows in particular that 

ax[v[(z',z')] = —vi(z,z), 



and in general 

a x [v' n (z',z')} = X2(n+1) v n (z,z), 

a x [w' n (z' , z')} = x ^ n+1) w n (z,z). 
The general expressions follow from the recursion relations of chapter IIVI □ 



CHAPTER VII 



Physical applications 

The considerable modification and generalization of existing quantum electrody- 
namics in strong external fields described in the previous chapters raises two important 
issues. The first is whether the framework presented in this thesis is concrete enough 
to describe physical situations. One may ask how to apply the general construction of 
states of chapter [TT1 together with the local causal perturbation theory of chapter IVII 
The other issue is whether our improvements lead to predictions which differ from the 
conventional ones. In this chapter we will address the first issue while the second must 
be deferred to future investigations. 

The local causal quantum electrodynamics constructs the local 8-matrix which 
describes the evolution of observables. This S-matrix in each order of perturbation 
is composed of local Wick and time-ordered products, defined in previous chapters. 
When smeared out with a test function g(x), all these elements become operators 
belonging to the algebra W. In order to be able to investigate the physical properties 
of matter and radiation it is necessary to specify the representation of this algebra and 
identify in some way the states on this algebra with physical situations. In a concrete 
representation the non-linear quantities will become operators on a Hilbert space. 

The non-linear observables were constructed in such a way that the GNS represen- 
tation based upon a state u of the free field algebra (X can be extended to the enlarged 
algebra W. This extension is possible only, if the two-point function of the state u is 
of Hadamard form (chapter IVjl. 

The "base state" to carries with itself the concept of excitations. If to is prescribed by 
a projection operator P + , the states of the form ?/>(/_ ip*(f + )Q with f± = P±f 3 TC 
describe excitations which would be conventionally called particles. It is the aim of 
interacting quantum electrodynamics to investigate the evolution of such excitations. 

The identification of mathematical objects with physical reality always has its lim- 
its; they are particularly exposed in the interacting theory. Suppose, for instance, 
that the external field is static. Then the excitations can describe particles in bound 
states (if such states are allowed by the external field). But what are these states? 
The quantum electrodynamic interaction changes the dynamics of the Dirac field, and 
what was stationary in the free theory is not stationary any more in the interacting 
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theory. We would like to know the energy of the excited state 2P of an electron in a 
hydrogen atom. We describe this state by ip*(f 2 p)fl, where f^p is the wave function 
of the 2P state calculated with the help of the Dirac equation. However, because this 
state is not stationary, we will never be sure what physical situation it corresponds to 
or how to produce such a state in experiment. 

With this remark in mind we begin to address the first of the issues named at the 
beginning; in what follows we shall show, how the theory developed in this thesis is 
applied to describe the simplest non-trivial physical effects. 

VII. 1. Electrodynamics in the presence of a static background 

Static backgrounds are important, mainly because they provide a sufficiently well- 
characterized and isolated regime for a study of the effects of the QED developed 
in this thesis. Not only the fully ionized heavy elements (see |Sto98j) but also fine 
atomic/ionic traps (see [Deh90 and many others) allow for a study of the electro- 
dynamic effects with great precision. In many of these experiments a single electron 
probes the external field effects. Perhaps for those experimental reasons the con- 
ventional literature on the external field QED has focused on static external fields 1 
jMP898Ll5na02] . 

VII. 1.1. Vacuum representation, static background. The construction of 
the vacuum representation in the case of static external backgrounds has been treated 
in chapter III. 31 Here we review it from a more heuristic point of view and restore the 
appropriate spinor indices. 

The free- field operator in the static external fields "at a point x" is an operator- 
valued distribution which can be decomposed according to: 

^\ x ) = J d 3 p |^(x,p) e - lE+ip)t a s (p) +<Hx,p) e ~ lE ' {p)t 6* s (p)} + 

N 
n=0 

Here we have separated the bound states from the scattering states. The w(x, p), 
f(x, p) denote the positive/negative- frequency scattering states (the index p is con- 
tinuous). The w n (x) denote the spinorial wave functions of the bound states, which 
we have assumed to have positive energies (for simplicity); there may exist an infinite 
number of them (N — > oo). The symbols E ± (p) denote the energies (the eigenvalues of 



Another reason might be the sense of uniqueness that the ground state vacuum carries with itself. 



the time-independent Hamilton operator) of the positive/negative- frequency scatter- 
ing solutions, respectively. The creaton/annihilation operators obey the standard CAR 
relations. The vacuum Q is the state annihilated by all of the annihilation operators: 

a(p)fi = 0, a n Vl = 0, 

6(p)fi = 0. 

VII. 1.2. First-order processes, creation of the electron-positron pair. In 

quantum electrodynamics with the external field treated as a perturbation it is known 
that one photon (with sufficiently large energy) can create an electron-positron pair 
[BLP82 i . In what follows we will investigate this process within the framework of 
local QED. In the zero-th and first order the evolution operator § is given by 

d x g(x) : ip{x)^ tJ ''il){x) : A^x), 

where the test function g(x) equals one in the interaction region 2 . In order to investi- 
gate the process in which a one-photon state evolves into a state with fermionic content 
("creation of an electron-positron pair") we will calculate the matrix elements of the 
S-operator with respect to the following states of the free fields: 

the initial state : I = Q <S> \f), 
the final state : F = \u, v) Cg> Q, 

where / denotes the photon's wave packet and u, v denote the wave packets of the 
electron and positron respectively (they may, for instance, be wave functions of some 
bound states). We obtain 

(I,SF) = J d A xg{x) (n,:^(x)j^(x) : \u, v)) ■ (|/), A,(x)Q) 

with 

I f d 3 k „ s _ ikx 



1 f d 6 k — 
(\f),A,(x)n) = —3 / -= e«(k)/ a (k)e" 



What remains to be evaluated is the matrix element of the Wick product of Dirac field 
operators: 

(Dtrac) = (fi, : ^ A {x)^ B {x) : \u, v)) 7^, 

where we will have to employ the local version of the Wick product. The final state 
of the Dirac field is orthogonal to the initial state (the vacuum), and therefore the 
distribution which has been chosen to regularize il)(x)ip(y) will have no influence on 



Later we shall take g = 1 on the strip t € (0,T) (i.e. we shall abandon the spatial smearing) and 
assume g to vanish just outside of this strip. 



the result. Consequently, the amplitude will turn out to be equal to the standard one 
calculated with the help of usual formulations of external field QED AB65J. We get 

(Dirac) = v(t, x) 7^ u{t,x), 

where the classical Dirac Hamiltonian governs the time-evolution of the wave packets. 
Finally, 

(I,S F) = * f d*x g(x) [ Mt^) ^(k)«(t,x) 7 Q (k)e-^. 

Now we can explain, why it is possible to abandon the spatial smearing of g(x). The 
electron and positron wave packets are intended to be regular wave packets 3 . Wave 
functions of such packets vanish faster than any power of l/|x| at spatial infinity. There 
is therefore no need to further strengthen the decay behavior at spatial infinity. We 
take g(t,x) = 1 for each t G (0, T). 



VII. 1.3. Second-order processes, consequences of the redefinition of the 
Wick product. In free quantum electrodynamics on external backgrounds the elec- 
trons may reside in bound states. Such states correspond to some discrete eigenvalues 
of the free Dirac Hamiltonian which lie in the interval (—mc 2 ,mc 2 ). We shall investi- 
gate the process of annihilation of a bound electron with a scattering state of a positron 
in a process, where two photons are emitted. 

We shall use the local version of the second-order time-ordered product l|VI.9|) . As 
there will be two incoming fermions and two outgoing photons, only the summands 
with the Wick product of two Dirac and two electromagnetic field operators will have 
a non- vanishing expectation value. In the case y ^J + (x) they are 



-e 2 



H AD {x,y) : foWMv) : +H BC (x,y) : ^ A (x)^ D {y) : ^ AB ^ CD : A^x)A v (y) :, 



where we have already utilized the local version of the Wick product. The initial state 
is defined to be a two-particle state: the ground electron state and some scattering 
state of the positron: 



|I> = a* b*(m, 

where 

&*(/) = J d 3 p f* (p) 6* (p), a = 1, 2, fe P.H, 



We may also take as u or v the bound state wave functions (if such bound states exist in a given 
external field background). 



and the two polarization functions f a are normalized by 

J £ P F(P)/ Q (P) = 1 

The outgoing state contains two photons of with certain polarizations: 

\F) = I/1/2) = J d 3 p d 3 k (p)/?(kte(p)a£(k)n, 

where a* (k) denotes the creation operator for the a-th polarization. 

The amplitude of the process at hand requires us to calculate the following expec- 
tation value 4 : 

A 1 = (\I)®n,T 2 (x,y)Q®\F))= (VII.l) 

= (|J), : ^ A (x)(YH F Y)AB^ B (y) : O) ■ (O, : A^)A(y) : |F» (VII.2) 

(the other summand of T2 will be considered later; here Hp denotes the Feynman 
parametrix which is the Hadamard parametrix with IV). The electromagnetic expec- 
tation value clearly gives 

(n, : M*)Mv) ■ W) = ^rvz [ 4== e~^ + ^ + x ~ y ]- 

2 ( 27r ) J vMffl 

• /i, a (qi)/2,/3(q2)- 

The Dirac-field expectation value, 

(Dirac) = (\I),:lf A (x)i/; B (y): Q), 

turns out again to be indifferent to the modification of the Wick product (this is 
because the initial state is orthogonal to the final state). We obtain 

(\I),:^ A (x)i; B (y) : fl) = u A {x) J d 3 p f a (p)v B (y,p), 

where the time dependence of the wave functions again is governed by the Dirac equa- 
tion. 

The amplitude of the process differs from that of usual formulations of QED in that 
the Feynman parametrix Hp(x,y) stays in place of the Feynman propagator Sp{x,y). 
The former is strongly local, the latter is not. In order to calculate the amplitude one 
should use the expansion of Hp in powers of Y — (x — y) (x — y) + ie, similar to that 
of the parametrix described in chapter IIV1 



Here we denote the electromagnetic and the Dirac vacuum with the same letter, as this should cause 
no confusion. 



VII.2. Outlook 



The investigations of this thesis point into three directions which are in our opinion 
good candidates for future research. In what follows we shall describe them briefly. 

(i) The local definition of the current operator, eq. IVI.13[ allows for a calculation 
of the expectation value of the current in some given state of the free Dirac 
field. We expect a non-trivial result, if this quantity is computed for a ground 
state of some external field configuration. This bears some resemblance to the 
investigations of Casimir effect in the framework of strongly local quantum 
field theory |Ke03j . 

(ii) The topic of (perturbative) physical properties of interacting matter has only 
been touched by the investigations of this chapter. It is important to under- 
stand, how the usual predictions of external field quantum electrodynamics, 
derived for instance in chapter V of the textbook of Akhiezer and Berestet- 
ski [AB65 , are altered, if we use the local theory developed in this thesis. 
Furthermore, a detailed investigation of the dynamics of the excited states 
of an electron in a hydrogen atom appears possible. Issues like the shift of 
energies of the states or the time dependence of the excited state's amplitude 
and phase with today's methods are of direct experimental significance. 

(iii) One of the improvements our thesis has brought is the possibility to handle 
time-varying external fields. They can be investigated, because the states 
of the Dirac field on the surfaces of constant time are not assumed to be 
globally equivalent to one another anymore. Thus, if a deep potential well 
is turned on extremely quickly, the state of the free Dirac field, which was 
the vacuum state in the past, evolves into some state in the future 5 . It 
would be important to know in what way this picture is altered by quantum 
electrodynamic interaction. The interaction of the radiation field with the 
Dirac field presumably changes the final state of the Dirac field, but it is 
impossible to tell without calculation whether this change is significant. It 
is our hope that models of this type can shed some new light on the process 
of collapse of matter in the formation of black holes in the theory of general 
relativity. 



It would be premature to call such a state excited, because in time-dependent situations there is no 
preferred reference state. 
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APPENDIX A 



The electromagnetic units 

As is well known, it is possible to construct a dimensionless quantity out of the 
unit of electric charge, h and the velocity of light: 

e 2 1 



he 137.036' ^'"^ 
Consequently, it is possible to eliminate all electrical artificial units like Coulomb or 
Volt and express them in terms of mechanical "cgs" units only. In Gauss units the 
attraction force of two charges is 

f = ^. 

where the distance is measured in "cm" and charge in charge units i.e. one electron 
possesses the charge of 



e =W^- = 4.803- 10-° -/ gCm3 



137 V s 2 

The electromagnetic field potential, typically expressed in Volts, now has the unit of 



[v] = 



A.l. Action of the Maxwell-Dirac electrodynamics 

The unit of action is the same as the unit of Planck's constant: 

[S] = [h] = erg s = i^. 

s 

The Lagrange density C 

S = J dx° d 3 x C, 
where x° — c • t, has the dimension of 

g 



[4 



2 ' 

enr s 



which happens to be the dimension of the energy density. The Dirac field has the 
dimension of 



V cm 3 ' 
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so that the probability density is measured in 1/cm 3 : 

p = ifj*ijj. 

The action of the classical electrodynamics coupled to the classical Dirac field, 
which 

• contains the appropriate combination of physical constants, 

• has all the appropriate signs, 

• leads to the Maxwell-Diraca system, 

is 



S 



d 4 x 



r-r „ , e 



[ih %p^ a d a ip - mc ipip] + wfty -Ail - 

c 



167T C 



-F nh F ab 



It leads to the Maxwell-Dirac system 



d b F ba (x) = Anf(x) 



r 



ihd a A a (x) ip(x) = mc tp{x), 



where the electromagnetic current is 

j a (x) = eipj a ip. 

A. 2. H = 1 = c, particular combinations of electromagnetic quantities 

In order to use Plank's units, H = 1 = c we note that setting c = 1 is equivalent 
to s = cm. Subsequently, h = 1 means g = 1/cm. In such a way the electric charge 
unit becomes dimensionless. The interaction term in the action and, therefore, the 
first order time-ordered product has the dimension of 



ipj a ip -A a 
c 



[1 / cm 4 



due to 



[A] = 1/cm, [ip] = l/cm 3/2 . 
The two-point function has the same dimension as the current density 

[(fi,^0 fi)] = 1/cm 3 . 

Not surprisingly, all the dimensions above correspond to the rescaling powers in the 
transition to the rescaled quantities discussed in chapter IVI.8I 



APPENDIX B 



Microlocal analysis 



Microlocal analysis deals with pseudo-differential operators (\l/DO), distributions 
and their singularities. In this appendix we summarize the most important definitions 
and results. The purpose of our presentation is to exhibit the beauty and usefulness 
of microlocal analysis in a readable way. We shall thus not put much emphasis on the 
correctness of our notation; neither shall we present the most general versions of the 
theorems (which sometimes obscure their immediate value). 

Definition B.l 

A wave front set at x, WF x (u), of a distribution u is the closed, conic set of directions 
in the tangent space at x in which the Fourier transform of u, localized at x, does not 
decay rapidly. In other words 



Here <p(y) is a Qf-function with support in a small neighborhood around x. The wave 
front set WF(u) is simply 



the smooth function a(x, k) is called a symbol of A. The \l/DO is of order m, if its 
symbol fulfills 



for some constant C = C(f3) and all (3. The derivatives with respect to x do not affect 
the decay behavior in \k\. 

Remark. It follows from the above definition that the symbols have to be C°° in 
both of their variables. Quite often, however, one would like to use operators the 



k 0VF x (u) ((pu)(Xk) decays rapidly for A — > oo. 



WF(u) = {(x, k): ke WF x (u)}. 



Definition B.2 



A pseudo-differential operator A is a linear operator defined via 
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symbols of which possess a discontinuity at certain points (eg. symbols of the sort \k\ 
or A/|k| 2 + m 2 ). 

Definition B.3 

An asymptotic expansion of a symbol a(x, k) is a sum of symbols of the form 

p m (x, k) + p m -i(x, k) + p m - 2 (x, k) + ... 

such that 

— oo 

a(x,k) - ^2p n (x, k) 

n=m 

is a symbol of order — oo. Here p n (x, k) are symbols of order n homogeneous in k. The 
first of those symbols, p m (x, k), is called the principal symbol of A. 

Definition B.4 

The \l/DO A is properly supported, if its distribution kernel 1 A(x,y) is such that 
A(K,y) as well as A(x,K) for a compact K have compact support 2 . 

The following theorems are of great practical importance: 
Theorem B.5 (Pseudo local property) 

Let A be a ^DO of any order m and let u be a distribution. Then 

WF(Au) C WF(u). 

Theorem B.6 (Propagation of singularities) 

Let A be a ^DO of order m which is properly supported and has a real principal 
symbol p m (x, k). If u is a distribution which solves the (in)homogeneous equation 

Au = f 

for feC™, then 
(i) 

WF(u) C a -1 (0) \ {0} = {(x, k) : a(x, k) = 0}. 



From the Schwarz kernel theorem every 'J'DO may be represented as 

(Ai/))(x) = J dy A(x,y)ip(y), 
where A(x, y) is called the distribution kernel of A. 

2 The notation is obvious: A(K, y) has compact support, if A(x, y) has compact support for all x e K. 



(ii) If a certain point (xq, ko) is in the wave front set of u, then also the whole 
Hamiltonian trajectory of this point lies in WF(u). Such a trajectory is 
derived from the Hamilton equations of motion with p m (x,k) taken as the 
Hamiltonian. 

Quite often the wave front set of distributions has a form (x,p, y, —p); it is common 
to introduce a primed wave front set in order to get rid of the minus sign: 

Definition B.7 

The primed wave front set of a bi-distribution is defined via 

WF'(x,k,y,p) = WF(x,k,y,-p). 

With that definition we have: WF'[5(x,y)} = (x, k,x, k) C M x M. 
Definition B.8 

Let u(x,y) be a bi-distribution, x G X, y eY. Then we introduce the notation 

WF x (u) = {{x, k) : (x, k, y, 0) G WF(u) for some y G Y}, 
WF Y (u) = {(y, k) : (x, 0, y, k) G WF(u) for some x G X}. 

Theorem B.9 (Composition of bi-distributions) 

Let u, v be properly supported bi-distributions on X x X. The composition of them 

(u o v)(x, z) = J dy u(x, y)v(y, z), 

exists if 3 

WF^nWF^v) =0. 

In such a case we have 

WF\u ov)d WF\u) o WF'iy) U 

U [WFi(u) x (X x {0})] U [(X x {0}) x WF^v)} , 

where the composition of wave front sets above means: 

(x,k,y,p) G WF'(u)oWF'(v) ^3{z,q) : (x, k, z, q) G WF'(u) and (z,q,y,p) G WF'(v). 



3 Here WF^u) denotes WF' x {u) in the second variable. 



APPENDIX C 

Quantum Dirac field in the absence of any external potentials 

In the present appendix the free Dirac field in the absence of any external potentials 
will be quantized. The quantization procedure contains two steps: 

• Defining the CAR algebra. 

• Constructing a representation of the CAR. 

We will present both of them in a somewhat heuristic manner. 

C.l. CAR Algebra 

Consider the mode decomposition of the classical Dirac field: 



Here s — 1,2 are connected with the orientation of the spin of the electron/positron 
(they enumerate the basis spinors), and u A (p), v A (p) are the basis spinors 1 ; a s (p) and 
b* s (p) are ordinary functions. There holds 




>p K(P) « S (P) e~ ipx + vf(p) b* s (p) e**} . 



(C.l) 



$> s (p)>Mp)| = p+(p), 



s 



J2\vs(p))(v s (p)\=P-(-p), 



s 



which in spinorial indices means 




(C.2) 



s 




(C.3) 



s 



Capitular letters denote spinor indices. 
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The operators -P ± (p) project onto the positive/negative-frequency 2 subspaces of H. 
They can therefore also be written as 

where p = \H(p)\. The basis spinors supplied with the Fourier factor 

u(x, p) = e~ tpx u(p), v(x, p) = e ipx v(p), 

are solutions of the classical Dirac equation: 

{il a d a - m)i> = 0. 

From this equation we infer 

7°(2p )P-(p) «(P) = 

as well as 

7 (2po)P+(-p)^(p)=0. 

Thus 

P + (p)v(-p) = 0, 

p_(pHp) = o, 

which is in accordance with P ± being orthogonal projections. The basis spinors fulfill 
also: 

(u\u r ) =5 S \ 



(v 8 ,v r ) =5 sr . 
'., .) denotes here the classical scalar product 

(^,^) = J d 3 x ^7°^ = J d 3 x 



Thus far we have used the mathematical structure of the classical Dirac field. 
There is a Hilbert space 7i of four-component spinors, which are square-integrable with 
respect to the scalar product (., .) given above. The basis spinors are eigendistributions 
which correspond to the continuous spectrum of the free Dirac operator. 



The frequency is the eigenvalue of the Hamiltonian 

H = 7 °(^+m). 



The expression (IC.ljl is a basis for a more common quantization of the Dirac field. 
One introduces the ("sharp") creation/annihilation operators in place of the ordinary 
functions (Fourier components) a(p) and b(p). They are objects which fulfill the 
following anti-commutation relations, 

{a s (p), a*(q)} = 5 rs 6(p - q), 
(&*(p)A(q)} = S rs 5(p - q). 

When smeared out with square-integrable test functions 

a s (f) = J d 3 P P + (p)f(p)a s (p), 

bs(f) = J d 3 pP-(p)/(p)fc s (p), 

they fulfill 

{a s (f),a* r (g)} = 5 rs (f+,g + ), 

{K(f),b r (g)} = S r8 (g-,f-). 

The anti-commutators between a and b are assumed to vanish. The Fock space is 
constructed in that one assumes that there is a vacuum O: 

a(m = 0, b(f)Q = 0, 

and one constructs the many-particle subspaces by a successive application of the 
creation operators on Q. 

At this point the representation of the CAR algebra (X, that is the algebra of the 
polynomials of the smeared-out field operators (jC.l|) . has been constructed. We may 
verify the CAR property explicitly 

ftKx), V*(y)} = (2^8 / d 3 P J2{u s (pK(p)e^-^ + ^( P K( P y p (- x +^} = 

= J^fJ d 3 p{P + (p) + P-(p)}e^- y) = 5(x-y), 
where in the positron part the change of variables p — > — p was performed. Therefore, 

{^(x),^(y)} = 5(x-y), (C.4) 

which can be written as 

W),r(g)} = (f,9), 



where as before (., .) denotes the scalar product in the classical Hilbert space H 
is the analogue of (|II.9J) ). 



APPENDIX D 



Model of the spontaneous atomic emission of light 

This appendix contains a field-theoretical model of the emission of radiation from an 
excited state of an atom. Although this problem has been attacked by many authors in 
the past we believe that our model describes this fundamental process in a better way. 
A huge part of what has been done on this subject relies on the perturbation theory 
w.r.t. the interaction of the electromagnetic field with the electron under consideration. 
Such a treatment can only give reliable information about the beginning of the decay 
process, but that is not what we are aiming at. On the other hand, a different approach 
due to V.Weisskopf and E.Wigner establishes a closed integro-differential equation for 
the time-dependent amplitude of the excited-state part of the electron's wave function. 
Our model is a refinement of a version of this approach which has been described in 
the textbook of M.Scully and M.Zubairy 1 SZ97J. Our model goes further because: 

• it is constructed with the help of the field-theoretical methods, describes a 
unitary time evolution of the system; 

• in the derivation it is possible to keep a freedom of the initial state of the 
radiation field. It can be recognized at the end that the spontaneous decay 
is caused directly by the two-point (autocorrelation) function of this state. 

In particular, we regard the second point to be very important. In recent years many 
attempts to experimentally alter the spontaneous emission process have been devised 
(eg. the presence of a squeezed state PCK92J or a light reflecting cavity). Those ex- 
perimental setups influence directly the autocorrelation function of the electromagnetic 
field, and thus can be easily investigated in our framework. 

The spontaneous emission of radiation is important for the general context of this 
thesis for the following reasons. Firstly, it provides a testing device, a detector of the 
state of the quantum field. Indeed, with a help of a single isolated atom it is possible, 
for instance, to investigate the properties of the KMS (thermal) states of the radiation 
field. With an ensemble of such atoms it is even possible to measure its temperature 
[MaQl]. The model can be of help in the investigation the sub-vacuum fluctuations 

1 The original paper of V.Weisskopf and E.Wigner is not easily accessible (published in 1930). It 
makes use of the methods of the early days of quantum mechanics. The modern expositions of that 
approach often use dubious mathematical tricks as in S Z97j . 
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of certain exotic states of the radiation field. The other issue is that our model of the 
spontaneous emission gives evidence of what happens in an interacting field theory. 
Apart from the ground state (whose dynamics are also altered) all other "excited 
states" of the free Dirac field are unstable in the interacting theory. Accordingly, one 
should exercise caution with the claims about the properties of the excited states and 
their energy shifts M P^98| ITTFS99] . 

The problem of an interaction of non-relativistic bound-states with the radiation 
field has been investigated with great rigor by V.Bach and collaborators (see for in- 
stance BFS99 ) 2 . In comparison to those investigations our model is a modest attempt 
- a special case at most. It provides, however, a relatively simple evolution equation 
which gives an insight into the time dependence of the system. In our opinion such an 
insight complements the investigations of V.Bach. 

D.l. Hilbert space and the interaction 

Consider a system consisting of an electron with two bound states (non-relativistic, 
described by quantum mechanics) coupled to the quantum radiation field. The natural 
Hilbert space for such a system is 

where T denotes the physical (transversal) Fock space of the Maxwell field. Later we 
will restrict this Hilbert space to its subspace Hi spanned by the vectors 

V>i(x)<g>n, ^o(x)®|/), 

where V>o/i denote the wave functions of the two bound states, Q is the electromagnetic 
vacuum, and |/) denotes the one-photon state 

|/> = a *(/)n = J d 3 P r (pK(p) n. 

Here f a (p) G S is a test function which describes the a-th polarization component of 
the photon's wave packet. 

In order to describe the interaction we introduce the interaction Hamiltonian 3 , 

V(t,x) = -— Pi ®A l (t,x), (D.l) 



We are grateful to K.Fredenhagen for pointing out this reference to us. 

3 In what follows we shall use the units K = 1 = c (the speed of light in this chapter is denoted by 
a bold c). We shall also omit the factor — e/mc for brevity. We note that in this chapter p means 
Po = c|p|. The correct factors as well as the atomic units will be adapted at the end of the calculation. 



where pj = —ihdi and A l (t, x) denotes the electromagnetic field operator in the ra- 
diation gauge with the free time evolution already implemented 4 . This is a standard 
interaction Hamiltonian describing the interaction of the non-relativistic systems with 
the electromagnetic field under the assumption that A 2 can be neglected, what we 
hereby also assume. The electromagnetic field operator is expressed in terms of the 
creation/annihilation operators (cf. chapter ITTT|) : 

1 f d 3 k 

A(x) = —3 / -= ef (k) R(k) e ik * + a a (k) e~ ik *} . 

We adopt the interaction picture w.r.t. the free electronic Hamiltonian which is 
given by p 2 /2m, the unperturbed energies of the states ipo,i will be denoted by E ^. 
In order to develop the model we now restrict the Hilbert space and the Hamiltonian 
to the subspace TCi. By doing so, we shall obtain a closed quantum system with a 
bounded, selfadjoint Hamiltonian. It is another matter to what extent our restriction 
describes the physical situation well. We expect that, as long as ^ ,i are the lowest two 
bound states which differ by ±1 in angular momentum, our assumptions are reasonable. 

We start by writing the time-dependent state vector S(t), which is an element of 
TC%, with the help of c(t) and /*(p): 

S(t) =c(t) fa(t,x)®n + fa(t,x.)® \f t ). 

The interaction (jD.lj) leads to the following evolution equation 

ic(t) e' iElt fa ® n + ie~ lEot fa ® \f t ) = 
c(t) e~ im Pl fa g) |^(t,x)fi) + e- lEot Pl ij ® \A%x)f t ). 

Here, however, what stays on the RHS does not belong to 7i\. We restrict the RHS 
to Tii in order to obtain a closed system. We do it by contracting (taking the scalar 
product) the above equation with the vectors fa ® and ® | fl /3(k)fi) separately 5 . 
We obtain 

ic(t) = e lujt (</>!, a} (t,x) p^ ), 

where 6 

4(t,x) = (n,A%x)a*(j)n) = ^ +ipx / t a (p)4(p)> 



4 Which means that the time dependence of A 1 is generated by the purely electromagnetic Hamiltonian 
Jd 3 p p a*(p)a{p). 

5 The operator a(k) is taken sharp here; however, this is allowed, as the use of a(k) still leads to 
expressions which make sense as distributions. 

6 This comes from the contraction of the negative-frequency part of the electromagnetic field operator 
with the /((p)-smeared creation operator. 



and uj — Ex — E . In the calculation we have used the fact that (Q, A 1 Q) = 0, which 
is true for any quasi-free state Q, in particular also for the vacuum. 
The second equation reads: 

iff, t (k) = c(t)e-^ (^o, (n, ai3 (k)A%x) Q) p^o) 
= c(t)e-^ Vo, P^o , 



2fc 

here we have again utilized the assumption that Q is a quasi-free state: the expectation 
value of an odd number of creation/annihilation operators vanishes for such states. The 
particular form of the operator in the bracket above, 

gifct-ikx e W k ) 



'2k 

is true only, if (Q, ap(k)a* a (p) Jl) = 5(p — k)<5 Q/ 3, i.e. if f2 is the vacuum. For other 
states Q we would get some other expression, but the analysis which follows proceeds 
analogously. 

We introduce some convenient abbreviations: 



Xi(p) = J d 3 x e ipx Vi(x)p i Vo(x), 
y»i(p) = y d 3 a; e ipx ^ (x)p^i(x). 



If the external field which binds the electrons is smooth, then it follows from the 
elliptic regularity that Xi(p) i s a smooth function of rapid decay. The vector index i 
of XijPi is always contracted with polarization vectors e l a (p). If only the transversal 
photon polarizations are emitted (thus a = 1,2), then all the polarization vectors are 
orthogonal to the vector pj: 

Pi £(p) = 0. 

Then 



<MPJ = XAP)- 

The system of differential equations that needs to be solved is 
ic(t) = e™ J ^= e-* UP) efa) ff(p), 

pipt 

iU,t{p)-c{t)e-^-j= X i{p) 4(p). 



It is not difficult to see that it is the Schrodinger equation on Tli with the restricted 
interaction Hamiltonian 

H \ n = PVP, 



where P is the projection from T onto Moreover, the restricted Hamiltonian 
is bounded (because the function x(p) decays rapidly for large p) and symmetric, 
therefore selfadjoint. 

The initial value condition (no one-photon component at t = 0) is 

/ t=0 (p) = 0, c(0) = 1. 

The equation for /g,t(p) may be integrated: 

/t ips 
ds c(s)e-^ s ^= ^(p) e*,(p). 



Upon insertion into the first equation of the system, one obtains 



Suit 



c(t) = -e tuJt \ ds c(s)e~ tuJS S(t - s), (D.2) 

where 



o 



s« - ') = I d ^' lrl - ] ( - «X,(P) 



2p V P 2 

is the Xrsmeared transversal two-point function 7 . The equation (jD.2|) together with 
the initial condition c(0) = 1 describes the atomic emission of light. The equation 
(jD.2|) is well-defined. The function S(t — s) in the limit s — > t gives a finite number, 
as x is a smooth function 8 of rapid decay 9 . 

We also note, that S(t — s) decreases for large t — s; this is important in the 
investigation of the behavior of c(t) for large t. 

The equation (jD.2|) may also be brought to a more familiar form. We may integrate 
it, namely, from zero to T and (with an appropriate change of variables) obtain 

C (T) = 1 - I Z(T-s)c{s) ds, 
Jo 

where 

pT 

-iwt 



Z(t) = [ dt Stye- 
Jo 



The above function is a smooth function of its argument. Therefore we remain with a 
task of solving the Volterra type integral equation of second type with a smooth kernel. 



7 We have made use of x = The term (S^ — p'pVp 2 ) comes from the contraction of two transversal 
polarization vectors: ^2 a e* (p)e^(p). The assumption on the transversality of the emitted radiation 
leads to this factor. 

8 Thus the infrared point p = is an integrable singularity. 
9 The decay property makes the ultraviolet part p — > oo harmless. 
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Figure D.l. Spontaneous emission of light from an atom. The figure 
presents a numerical solution for some artificial values of various param- 
eters. The presence of oscillations at short times accompanied by decay 
for large times appears to be a general feature of the model (the number 
of visible oscillations and the decay rate are parameter-dependent). 



It is difficult to say anything general about the solutions of equation (jD.2|) . The 
numerical experience with it shows that the solutions seem to oscillate on a short time 
scale and then later decay to zero. Figure ID. II shows one of the solutions (for some 
artificial data, however). 

In order to investigate the concrete atomic emission processes we introduce atomic 
units. Formally this means that we set e, m, h equal to one and note that 

• r = 1 means r = 0.53 • 10~ 8 cm (N.Bohr radius), 

• E = 1 means E = 27.2 eV (typical atomic ionization energies), 

• t = 1 means t = 2.4 • 10 -17 s (typical atomic oscillation times), 

• c = 137 (result of the relation e 2 /he = 1/137). 

With those units it is possible to investigate realistic systems. One realizes that the 
only non-trivial constant is the speed of light, c which appears only in the expression 
of the smeared two point function S(t — s). We have 



D.2. Atomic units, comments, outlook 




It is not difficult to find the explicit expressions for Xj{p) i n case °f the harmonic 
oscillator or the hydrogen atom wave functions. The numerical analysis, however, 
stumbles upon the following difficulties: 

• The typical, experimental, decay times in case of the 2P — > IS transition are 
of the order of nanoseconds. Such a time corresponds to 10 8 in atomic units. 
This raises the question of the propagation of the numerical errors. 

• The appearance of the speed of light in the function S(t — s) makes it both 
small and supported only for very short \t — s\. This even further worsens 
the problem, as the time step must be accordingly small. 

Those difficulties, however, could have been anticipated from the very beginning: 
the model constructed here describes equally well the emission of X-rays from the 
atomic nuclei and the emission of rf radiation from the oscillatory levels of molecules. 
In terms of the decay times those processes are separated by a vast gap of more than 
ten orders of magnitude. 

In the case of the hydrogen atom bound states 15* and 2P the smeared electromag- 
netic two-point function is: 



where A = 3/2, all in atomic units. 

The main merit of the model constructed above (apart from its well-defined math- 
ematical status) lies in the appearance of the two-point function of the "base" state 
of the radiation field on the right-hand side of equation (|D.2|) . Thus, it is clear that 
the "vacuum fluctuations" cause the decay of the atomic energy levels. Indeed it is 
relatively easy in the presented framework to investigate the spontaneous emission in 
the presence of a thermal state or a squeezed state of the radiation field. 

For the latter type of states such an investigation may have important consequences, 
namely, it could clarify the issue why the squeezed states produce atomic spectra with 
lines narrower than the natural line width 10 [PCK92 . 



This arises presumably, because the squeezed states exhibit a sub- vacuum level of fluctuations of the 
electromagnetic field. However, the periods of reduced fluctuations seem to last no longer than f 0~ 14 s 
(they are bounded by quantum inequalities Ma02 ), which is far less than the usual spontaneous 
emission time of the energy levels. 




APPENDIX E 



GNS construction and thermo-field dynamics 

In the chapter devoted to the construction of representations of the CAR algebra Ol 
we have indicated that interesting things occur, if the integral kernel B of the two-point 
distribution uj is not a projection. 

We shall consider a simple example of this here. Let H = C 2 and 



B 



1 

If we look for the Gelfand ideal in the one-field operator subalgebra of Ol, we will find 
that 

V>[(a,0)] 

is the only (one-field) annihilation operator present there. 

In order to obtain a representation of Ot we may pretend that we have the cre- 
ation/annihilation operators, as usual, with the vacuum: 

a(VBf) n = 0, 



b(Vl-Bf) n = o, 

but those operators will in general no longer be elements of Ot. In the standard Fock 
space built with the help of a*, 6*'s we may construct the representation of Ol via 

^(f) = a(VBf) + b*(VT^Bf). 

The representation defined as above turns out, however, to be reducible; there are 
operators different from the identity which act in the representation space and commute 
with the representation of Ol given above. 

Much better insight into the nature of the problem is gained, if we extend H to 
H2 = H^H and choose 

t , b VbVT^" 

Bo 



By/T^B 1 - B 

which is a projection operator on 7i 2 irrespective of B. We shall now search for a 
representation of a doubled algebra generated by the elements ^(/), / G H2, and their 
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adjoints. In our example 

/l 

1/2 



\ 

1/2 







\0 1/2 



1/2/ 



If we introduce an orthonormal basis of Ti.2'- 



u 



V 



e 



e 



(1,0,0,0), 
(0,0,1,0), 
(0,1,0,1)/V2, 
(0,1,0,-1)/V2. 



then 



B 2 = P e + P L 



where P's denote the respective 1-dimensional projections. The Gelfand ideal reveals 
four annihilation operators: 



which obey the standard anti-commutation relations among themselves. An irreducible 
representation of the doubled algebra is thus given on the Fock space constructed with 
the help of the adjoints of the above operators. What is interesting is that on such a 
Hilbert space there acts also the representation n(Ot) of the original algebra (X with 



This field and also its adjoint anti-commute with the representation of (X. Thus, the 
observables constructed with their help 1 all commute with the elements of ir(Ol). 

The doubling of 7i (and in turn of Ci) is frequently done in the literature on 
"thermo field dynamics". Further discussion of the representation theory constructed 
on a thermal two-point function can be found in [Ha96j chapter V.1.4 , or |Fre99j . 



ai = *(e), 
ai = 




4>(g) = *[(0, 0, g+, = ^= (a a - b* 2 ) + g+ b\. 



The observable algebra is the algebra of even (i.e. containing even number of field operators) elements 

of a. 
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